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ABSTRACT 

We have simulated QCD using 2+1 flavors of domain wall quarks on a (2.74 fm)-^ volume with 
an inverse lattice scale of a^^ = 1.729(28) GeV. The up and down (light) quarks are degener- 
ate in our calculations and we have used four values for the ratio of light quark masses to the 
strange (heavy) quark mass in our simulations: 0.217, 0.350, 0.617 and 0.884. We have mea- 
sured pseudoscalar meson masses and decay constants, the kaon bag parameter Bk and vector 
meson couplings. We have used SU(2) chiral perturbation theory, which assumes only the up and 
down quark masses are small, and SU(3) chiral perturbation theory to extrapolate to the physical 
values for the light quark masses. While next-to-leading order formulae from both approaches 
fit our data for light quarks, we find the higher order corrections for SU(3) very large, mak- 
ing such fits unreliable. We also find that SU(3) does not fit our data when the quark masses 
are near the physical strange quark mass. Thus, we rely on SU(2) chiral perturbation theory 
for accurate results. We use the masses of the Q. baryon, and the n and K mesons to set the 
lattice scale and determine the quark masses. We then find /;;: = 124.1(3.6)stat(6.9)systMeV, 

/j^= 149.6(3.6)stat(6.3)systMeVand/j^//;r = 1.205(0.018)stat(0.062)syst. Using non-perturbative 
renormalization to relate lattice regularized quark masses to RI-MOM masses, and perturba- 
tion theory to relate these to MS we find m^(2GeV) = 3.72(0. 16)stat(0.33)ren(0.18)systMeV, 
mp(2GeV) = 107.3(4.4)stat(9.7)ren(4.9)systMeV, and m^j : m, = 1 : 28.8(0.4)stat(1.6)syst. For 
the kaon bag parameter, we find B^{2GeY) = 0.524(0.010)stat(0.013)ren(0.025)syst- Finally, for 
the ratios of the couplings of the vector mesons to the vector and tensor currents (fy and fy re- 
spectively) in the MS scheme at 2 GeV we obtain: /J//p = 0.687(27); /^*///f* = 0.712(12) and 
AV/^ =0.750(8). 



I. INTRODUCTION 

Numerical simulations of Quantum Chromodynamics (QCD) conventionally discretize four- 
dimensional space-time by introducing a lattice scale, a, which yields a well-defined path integral 
formulation appropriate for study via importance sampling techniques. Measurements of observ- 
ables at non-zero values of a can then be extrapolated to a = to produce continuum results and, 
if discretization errors are small, the continuum extrapolation is better controlled. An important 
aspect of simulating QCD is choosing a lattice discretization which reduces the a dependence of 
observables [1]. Here, one is helped by knowing that if an 0{a^) accurate gauge action is used, 
such as the Wilson action, and a discretization of QCD preserves the continuum chiral symmetries 
of massless QCD, then the lattice theory can only have errors quadratic in a, in combinations such 
as (aAqcD)^, and (ma)^, where m is a quark mass. These errors are much smaller than the 0{am) 
discretization errors which can occur if the lattice theory breaks chiral symmetry. Additionally, if 
chiral symmetry is broken at non-zero lattice spacing, renormalization of the simplest operators, 
while straightforward, is involved [2]. (For a recent review, see [3].) 

Preserving chiral symmetry at non-zero lattice spacing has a large impact on many of the more 
complicated observables, such as matrix elements of operators, that one wants to determine from 
lattice QCD. For observables such as hadron masses, which do not require any renormalization, 
controlling non-zero lattice spacing effects in the discretized theory is helpful. For observables 
requiring renormalization, such as weak matrix elements in hadronic states, the presence of chiral 
symmetry at non-zero lattice spacing can be vital to the renormalization and mixing of the relevant 
operators [4]. The chiral symmetries control the allowed operator mixings, so simplifications 
take place. Without this control, the number of operators and mixings amongst them can make 
renormalization very difficult, if not practically impossible. Chiral symmetry can now be preserved 
at non-zero lattice spacing with a variety of formulations [5-10]. (It is important to point out that 
some of the benefits of chiral symmetry can be achieved without its presence at non-zero lattice 
spacing. For example, twisted mass Wilson fermion discretizations [11] make judicious use of 
chiral transformations to calculate quantities without linear dependence on a and continuum-like 
renormalization properties, without having full chiral symmetry.) 

In addition to taking the limit a ^ 0, to achieve accurate physical results lattice simulations must 
reach the large volume {V -^ o°) limit and also the limit of physical light quarks. For discretizations 
which preserve chiral symmetry at non-zero a, simulations with arbitrarily light quarks can be 



done before taking the continuum limit - without chiral symmetry, lattice artifacts can alter the 
chiral limit at non-zero a and make the limits non-commuting. Current computer power does not 
allow simulations with physically light up and down quarks, so an extrapolation from masses used 
in the simulations to physical light quark masses must be done. Chiral perturbation theory [12, 
13] provides a theoretical framework for these extrapolations and, for lattice QCD discretizations 
which preserve chiral symmetry, the chiral perturbation theory is very similar to the continuum 
theory, since there are only a few lattice alterations to it. 

Of long recognized importance, the preservation of the global chiral symmetries in discretized 
QCD has been achieved by Kaplan's proposal [5] of four-dimensional fermions resulting from 
a defect in a five-dimensional theory. Further developments led to the domain wall fermions of 
Shamir and Furman [6, 7], that we use here, the overlap formulation of Neuberger and Narayanan 
[8, 9] and the perfect actions of Hasenfratz [10]. To date, the domain-wall formulation has proven 
to be the most numerically feasible. In this approach, one introduces a fifth dimension (which 
we label by the index s and which has extent Ls) and only achieves exact chiral symmetry in 
the Ls -^ °° limit. However, for finite Ls, chiral symmetry breaking effects can be made small 
enough to be easily controlled, as we will discuss in subsequent sections. The presence of the fifth 
dimension (we use L, = 16 in this work) increases the number of floating point operations required 
in a calculation by a factor of 0{Ls) over conventional QCD discretizations such as Wilson and 
staggered fermions, which do not preserve continuum chiral symmetries. 

The domain wall fermion formulation has been used extensively in numerical simulations for about 
a decade. The original works were primarily in the quenched approximation, although some early 
work did involve QCD with two light quark flavors [14]. More extensive 2 flavor simulations 
were done [15] and, with recent improvements in algorithms and computers, 2-1-1 flavor QCD 
simulations have been completed [16-20]. These previous calculations demonstrated that domain 
wall fermion QCD shows the expected consequences of having a controlled approximation to the 
full chiral symmetries of continuum QCD. In particular, the following important features were 
observed. 

1 . A mild dependence on a in the a ^ limit in the quenched approximation. 

2. For both quenched and unquenched QCD, the residual chiral symmetry breaking at non-zero 
a, as measured by the additive contribution to the quark mass m,,„ is readily made a small 
fraction of the input bare quark mass at practical values for L^. 



3. The operator mixing problem for domain wall fermion QCD is essentially the same as the 
continuum problem. 

These previous calculations were not able to fully exploit one of the main benefits of domain wall 
fermion QCD - the ability to probe the chiral limit of the theory at non-zero lattice spacing - 
and hence make accurate contact with the physical quark mass region. With recent advances in 
computers and algorithms, we have made considerable progress in simulating with light quarks 
and on large volumes. In this paper, we report on simulations of domain wall fermion QCD, with 
two light degenerate quarks and a single flavor heavier quark (a 2+1 flavor simulation) at a single 
inverse lattice spacing of a^^ = 1.729(28) GeV in a (2.74 fm)^ spatial volume. We used four 
different light dynamical quark masses in our simulations and the ratios of these light masses to 
the physical strange quark mass is 0.217, 0.350, 0.617 and 0.884. A single value for the heavy 
quark mass was used in all our simulations and its ratio to the physical strange quark mass is 1 . 150. 
(An accurate value for the physical strange quark mass was, of course, not known until after 
our simulation was complete.) For mesons made of degenerate light quarks, the corresponding 
pseuduscalar meson masses are 331 MeV, 419 MeV, 557 MeV and 672 MeV. We have also done 
measurements with a variety of valence quark masses, with a ratio of the smallest mass to the 
physical strange quark mass of 0.110, corresponding to a pseudoscalar meson with a mass of 
242 MeV. In a previous paper [16], we have given results from simulations with the same gauge 
coupling constant, but on a smaller volume, which gives us some understanding of finite volume 
effects. 

To extrapolate from our simulation quark masses to the physical values, we use chiral perturbation 
theory (ChPT), which is an expansion of low energy QCD observables in powers of the meson 
masses and momenta over the pseudoscalar decay constant. Within the general framework of 
ChPT one can consider only the pions to be light particles, yielding an SU(2)l x SU(2)r ChPT 
(which we will call SU(2) ChPT) or one can also consider the kaons as light, yielding an SU(3)l 
X SU(3)r ChPT (which we wiU caU SU(3) ChPT). In Section II, we discuss the DWF corrections 
to ChPT and develop SU(2) ChPT for the kaon sector, which we will later use to fit our data. 
In Section III we give details of our simulations, including the Rational Hybrid Monte Carlo 
(RHMC) algorithm that we use to generate our lattices. Section IV describes the sources and 
sinks we use for our pseudoscalar observables, our methods of determining desired quantities and 
results in lattice units. We will also use the mass of the Q. baryon, which has no corrections from 
chiral logarithms, as part of our scale setting and the details of the measurement of m^ are given 



in Section V. 

In Section VI we fit our lattice values for the masses and decay constants of pseudoscalars to 
partially quenched SU(2) ChPT at NLO. We find our data is well described by the theoretical 
formula from Section II, provided the pions have masses below about 420 MeV. We use the fits 
to SU(2) ChPT as the most accurate way to extrapolate our data to the chiral limit, since SU(2) 
ChPT does not require the kaon mass to be small, but only requires m-n <^ niK- Using the results 
for pseudoscalar masses from our SU(2) ChPT fits and our lattice values for the Q. baryon mass, 
we fix the lattice scale and bare quark masses using the known masses of the n, K and Q.. We find 
that our inverse lattice spacing is a^^ = 1.729(28) GeV. In a separate work [21], we have used 
non-perturbative renormalization to calculate the multiplicative renormalization factor needed to 
relate our bare lattice quark masses to continuum MS masses. We find 

m^(2GeV) = 3.72(0.16)stat(0.33)ren(0.18)systMeV, (1) 

mfS(2GeV) = 107.3(4.4)stat(9.7)ren(4.9)systMeV, (2) 

ihud ■■ ms = 1 : 28.8(0.4)stat(1.6)syst. (3) 

where (■ ■ ■ )stat5 (■ ■ ■ )ren and (• • ■ )syst show the statistical error, the error from renormalization and 
the systematic error. We assume the light quarks to be degenerate in this work. We now predict val- 
ues for /;, and /k and find ^ = 124.1(3.6)stat(6.9)systMeV and /k = 149.6(3.6)stat(6.3)systMeV. 
Our fits to SU(2) ChPT also determine the low energy constants (LECs) for pseudoscalar masses 
and decay constants in SU(2) ChPT. Furthermore, implications of our results to CKM matrix ele- 
ments are discussed. 

In Section VII we fit our light pseudoscalar data to SU(3) ChPT. Here we also find that our data 
is well represented by SU(3) ChPT at NLO, provided our pseudoscalars have masses below about 
420 MeV. The failure of NLO SU(3) ChPT to fit our data when pseudoscalar masses are near the 
physical kaon mass rules out using NLO SU(3) ChPT in this mass region. With light masses, 
we determine values for the SU(3) LECs which agree well with values determined by others. 
However, we find a small value for the decay constant in the SU(3) chiral limit, which we denote 
by /o (a complete description of our notation is given in the Appendix). Our fits give /o = 94 
MeV, with conventions such that the physical value is fj^ = 131 MeV, and this value is smaller 
than generally found phenomenologically, which we discuss further in Section VII. Along with 
this we find that the size of the NLO corrections to SU(3) ChPT, relative to the leading order term, 
are in the range of 50% or more. This makes the convergence of SU(3) ChPT for these quark 



masses unreliable. Thus, although it represents our data well and the parameters we find generally 
agree with others, we find the systematic errors in SU(2) ChPT substantially smaller and use it as 
our most accurate means of extrapolating our data to the chiral limit. 

In section VIII, we discuss our determination of the kaon bag parameter, Bk, which is needed to 
relate indirect CP violation in the standard model to experimental measurements. This section 
expands upon the data and analysis presented in [19]. Here we also find extrapolations to the 
physical quark masses to be under much better control with SU(2) than with SU(3) ChPT. We 
present our estimates of systematic errors, including finite size effects. We find 5^^(2GeV) = 

0.524(0.010)stat(0.013)ren(0.025)syst. 

In Section IX we present results for the couplings of light vector mesons to vector and tensor 
currents. The results for the ratios of the couplings of the vector mesons to the vector and tensor 
currents {fy and fy respectively) in the MS scheme at 2 GeV are: fp/fp = 0.687(27); fl*/fK* ~ 
0.712(12) and /;^//^= 0.750(8). 

II. CHIRAL PERTURBATION THEORY 

In this section, we discuss chiral perturbation theory for domain wall fermions in subsection II A. 
We develop SU(2) ChPT for kaons in II B. Our notation and an extensive list of the formulae from 
SU(2) and SU(3) ChPT that we use in this work are given in the Appendix. 

A. Chiral Perturbation Theory for Domain Wall Fermions 

Our simulations are done with domain wall fermions, which have explicit chiral symmetry break- 
ing effects at non-zero lattice spacing. These effects are controlled by the extent of the lattice in 
the fifth dimension, denoted by Ls, whose value is chosen to make such terms small, consistent 
with current computer power. The small chiral symmetry breaking that remains can be measured 
and its effects taken into account, as has been discussed extensively in the literature. A recent 
review of this topic is available in [22]. 

As previously discussed in [4, 23], for a theory with A'^ quarks these explicit symmetry breaking 
effects can be easily included by introducing an N x N matrix parameter Q. into the domain wall 
fermion action. This parameter connects four-dimensional planes at the mid-point of the fifth 



dimension and is included into the action by adding 

Sa = -E {^xX./2-l^L (^^ - l) ^xX./2 +^xX./2^i^ (^ " 1) ^xX./2-l } (4) 

to the conventional action for domain wall fermions [24] . Here, ^'x,s represents a five-dimensional 
fermion field with four spin components and suppressed flavor indices. We recover the conven- 
tional domain wall fermion action when we set Q = 1 . 
If we let Q. transform as 

a^VRavl (5) 

under SU(A^)l x 5t/(N)R, then the domain wall fermion Dirac operator possesses exact chiral 
symmetry. Thus we can use Q. to track the explicit chiral symmetry breaking from domain wall 
fermions, in both the Symanzik style effective Lagrangian for domain wall fermion QCD and in 
Green functions. However, we note that Q. itself is not a small quantity, since a general Green 
function involving the five-dimensional fermion fields does not have any approximate chiral sym- 
metry. However, for the effective action and low momentum limit of Green functions made from 
the four-dimensional fields at the boundaries of the fifth dimension, each power of Q. that enters 
should come with a suppression related to the ratio of amplitudes of the low-energy fermion modes 
between the boundaries and the midpoint in the fifth dimension. 

Consider a non-zero a effective Lagrangian description of QCD with domain wall fermions at finite 
Ls. The presence of the parameter Q. implies that the terms containing fermions, up to operators 
of dimension five, are 

Here nif is the dimensionless input bare quark mass in the domain wall fermion formulation, and 
C3 and C5 are dimensionless parameters that represent the mixing of left- and right-handed quarks 
between the five-dimensional boundaries. These parameters are of 0{e^^^') at weak coupling; for 
coarse lattices where there are localized dislocations in the gauge fields corresponding to changes 
in the topology, they are generically 0{aie^^^' + a2)/Ls, where a2 7^ is due to the density of 
localized topological dislocations [17]. 
Setting Q = 1 , we have 

ZmMf C3 

The combination Zminf + c^ is the total (dimensionless) quark mass and we choose L? to control 
the contribution of the second term, by changing the size of C3. Eq. (7) is identical to the result 



for Wilson fermions, except that the coefficients C3 and C5 are expected to be small, O(10~^), for 
realistic domain wall fermion simulations, compared to being 0(1) as for Wilson fermions. 
We can now discuss the application of chiral perturbation theory to our domain wall fermion 
simulations at a fixed lattice spacing. Our discussion will be for SU(3), but the results are easily 
generalized. We start from the conventional QCD SU(3) chiral Lagrangian in the continuum and 
make use of the presence of the Q. spurion field to add all additional terms to it. Initially we 
power count only in ^ and defer, for the moment, the additional question of power counting in 
a. We choose E(;c) = e2'<^W//o, where E transforms as: E -^ VlIIV^.Vl.Vr E SU(3). We define 
X = IBq diag(m„, m^, m.v), where Bq is one of the low energy constants (LECs) that enters in chiral 
perturbation theory. To 0{p^) the continuum Lagrangian is: 



^ = ^Tr[^^E^^Et] + ^Tr[2E+(2E) 



tl 



+ L^hv 



+ if^ |Tr[^^E(^^E)t 
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+ Lfhx 
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^^E(o)'^E)%rfe + (2E) 



Trfe+(2E)t 



+ Lf 



Tr (2E- (2E)t)] +4' ^Tr (2E2E+ (2E2E)^) . 



(8) 



For the domain wall fermion case, we can generate the new terms that arise by starting from 
the Lagrangian in Eq. (8). Since Q. transforms as 2 or E^ does, new terms can be created by 
substituting Q for 2 and E^ in Eq. (8), remembering that derivatives acting on Q. vanish. Since the 
dominant contribution of explicit chiral symmetry breaking in Eq. (7) is an additive contribution 
to the quark mass, we power count Q as 0{p^). Keeping terms of 0{p^) and using D,, D, for the 
LECs for these terms, we have additional contributions to Eq. (8) of 

:^DoTr[^E+ (^E)+] + :^DoTr[^2^ + ^^2] 

o o 



Z)4Tr 


'd^Ud^L)^' 


Tr[j^E+(J^E)''']+D5Tr 


(9^E(<9^E)+(^E+(r2E)'f) 


+De 


Trf^E+r^E)"^) 


Tr(2E+(2E)t)" 


+ ^6 


Tr(aE+(aE)'^) 


2 


+Z)7 


Tr('^E-(r 


iE)t); 


Tr(2E-(2E)t)" 


+ D7 


Tr[aE-(aE)'^) 


2 



+D8Tr I 2EJ^E + (2S^E)''' ) + DsTr ( J^EQE + (QEQE 



(9) 

Terms which involve two derivatives, a factor of 2 and a factor of Q will be 0{p^) and have been 
neglected. Note that we have kept the term involving Dq, even though it does not involve any E 
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fields. Such a tenii does play a role in determining the value for the chiral condensate through the 
variation of the partition function with respect to the quark mass. 

We would seem to have many new low energy constants to determine with domain wall fermions at 
non-zero a. However, the form of the Symanzik effective Lagrangian for DWF QCD shows that the 
leading order (in a) chiral symmetry breaking effect is a universal shift in the quark mass, i.e. c^ja 
multiplies the dimension three operator "y/i//. Thus, there is no difference in ChPT between nif and 
C3, and we can rewrite the terms in Eq. (9) in terms of the original LEC's of QCD, with a shifted 
quark mass, plus higher order corrections. Letting x = 25o[diag(m„, m^jm^) +C3diag(l, 1,1)] and 
looking at the D4 term as an example, we have 

D4Tr[<9^E(^'^E)-f]Tr(aE+(aE)'^) = L^ Tr [<9^E(<9^E)^] Tr (;tS+ (;tE)'^) 

+ 0{ac5)Lfhv [^^E(^^E)^] Tr (^E + E^) . (10) 

The last term is 0{a) x 0(05) x 0{p^). It is customary to power count 0{a) and 0{p^) terms as 
the same size for unimproved Wilson fermions, and this same term appears there [25], except that 
^Wilson {^Q(^\y While for Wilson fermions, this term must be kept at NLO in the chiral Lagrangian, 
for domain wall fermions, where C5 is very small, it can be neglected. 

Examining all the terms in Eq. (9), we see that the complete NLO chiral Lagrangian for domain 
wall fermions is given by Eq. (8), with X ^X^ where x is proportional to the sum of the input bare 
quark mass, nif, and the additive quark mass contribution which comes from C3. Since we will 
be working to NLO order of ChPT in this work, domain wall fermions at non-zero lattice spacing 
should be described by the chiral Lagrangian given in Eq. (8). When we fit our data for specific 
quantities to the chiral formula following from Eq. (8), the LEC's Li will differ at 0{a^) from their 
continuum values. Since we work at a single lattice spacing, we will not be able to correct for 
these deviations from continuum QCD. 

The size of the residual symmetry breaking terms represented by C3 and C5 in Eq. (6) is most easily 
studied by examining the five-dimentional current ^" which can be easily defined for domain wall 
fermions and is exactly conserved in the limit ruf ^ and L^ — > 00 [7]. The DWF equations of 
motion imply that this current obeys the divergence condition 

A^^;W = 2mfJlix) + 2r,^{x) , (11) 

where 7| is a pseudoscalar density made up of quark fields on the boundary of the fifth dimension 
and 7? is a pseudoscalar density containing quark fields at Ls/2 — 1 and Li/2. While the J^ term is 
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the result of the usual chiral non-invariance of the mass terai, the 7| operator is expected to have 
vanishing matrix elements at low energy as L^ -^ °°. It represents the effects of residual, finite-L^, 
chiral symmetry breaking. For low energy Green functions, the midpoint term in Eq. (11) can be 
expanded as 

7|^ ~ m^Jl - ^^^A^^; + c',01. (12) 

Here we have introduced the my^-independent parameter m,^„ related to the constant C3 in Eq. (6). 
We also have a new lattice operator, O^, similar to the axial transform of the C5 term in Eq. (6), 
which is carefully subtracted so that its matrix elements are of order a^ smaller than those of the 
operator 7| at long distances. Since the low energy matrix elements of the midpoint operator on 
the left-hand side of Eq. (12) will be suppressed by a factor exp(— aL^) we expect the quantities 
m,e„ (Z^ — 1) and Cg to all be of this order (at least in the perturbative regime). 
While an expansion such as that written in Eq. (12) is typically written in terms of the operators 
of the effective theory, the present form of this equation is useful because it can be combined with 
Eq. (11) to yield: 

Z^A^^;W = 2(m/ + m,j7f(.^) +40^ . (13) 

As is well known from the classic analysis of the case of Wilson fermions [2], the conservation of 
the vector current and the vector Ward identities imply that product my^7|(.)c) inEq. (13) approaches 
its continuum counterpart without multiplicative renormalization. This equation then implies that 
the product Z^£/"{x) will reproduce the standard continuum axial current when evaluated in low 
energy Greens functions. 

The appearance of the (Z^ — 1) term on the right-hand side of Eq. (12) and the consequent renor- 
malization of the five-dimensional DWF axial current is an effect that was not recognized in earlier 
RBC or RBC-UKQCD work. Such a possibility was raised by Steve Sharpe [22]. However, in this 
paper, Sharpe argues that (Z^ — 1) is expected to be of order m^^,, not the order m,^, scale suggested 
by the above argument. Sharpe's argument relies on the chiral character of the operators V^ ±A^, 
defined in terms of fields on the boundaries of the fifth dimension. The mixing between these 
operators implied by Z^ 7^ 1 can arise only if a fermion and an antifermion move from one wall 
from the other, a situation suppressed as m^^. We believe that this generally very useful argument 
does not apply in the present case because the right- and left-handed currents under consideration 
each span one-half of the entire five dimensions and are not localized on the left and right walls. 
For example, two fermions can propagate from the midpoint operator 75^ to the same wall with 
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only a suppression of [exp(— aL^/l)]^ = exp(— otLv). Note, this argument only applies to the per- 
turbative piece. For tunneling caused by near- zero modes of the four-dimension Wilson operator, 
two such modes are needed implying a suppression more like m^^^. 

Finally it is informative to consider the implications of this more complete analysis on the conven- 
tional calculation of the residual mass from a ratio of low energy pion to vacuum matrix elements. 
We measure a quantitiy commonly called m[^^{mf) and given by 

m (nif) = , , , , (14) 

(0| (m,J^ - (5^A;,< + c'gO") I k) 
= — ^^ C15) 

We can then use the lattice equations of motion to evaluate the A^.e/^ term so that Eq. (15) takes 

the form: 

1 c' / I \ (OlO"l;r^ 

<.K)='"» + (^-l)('"/ + '"...) + f(l + ^j^ (16) 

This equation governs the nif dependence of m[^^{mf). Since, by construction, the third term on 
the right-hand side is of order a^, we might be tempted to neglect it and use Eq. (16) to relate 
the nif dependence of m[^^{mf) to the difference (1/Z^ ~ !)• However, we can distinguish the 
dependence of m[^_,{mf) on the valence and sea masses, ni^?^ and m^.®^ (m^ and m/ in the general 
notation of the rest of this paper). We see that there is valence quark mass dependence in both the 
1/Z^ — 1 term and in the term containing Cg, while the sea quark mass dependence is only in the 
Cg term. As will be shown in Section IV B, m[^^ depends to a roughly equal degree on both m^^^ 
and my^. In particular, we find that d\n{m[^^ / dirC?^ ^ —4 and d\n{m[^^ / dmy^ ^ 6. This implies 
that the third term in Eq. (16), which can depend on both the valence and sea quark masses, must 
be of an equal size to the second. In fact, if we introduce physical dimensions, we expect the ratio 
(0|O5|;r)/(0|7||;r) to be of 0((aAQCD)^)- The effect of the quark masses on this ratio will be to 
replace one of the factors of Aqcd by the dimensioned factor m^?^ /a or w?f'^ ja. Thus, the quark 
mass dependent terms coming from the third term in Eq. (16) are suppressed by only by O(aAQCD) 
relative to the second term, which is apparently insufficient to permit them to be neglected, a fact 
alreadly pointed out by Sharpe[22]. 

The arguments above show that we expect Z,^ — 1 to be 0{m^^,) . If all of the unitary mass depen- 
dence of m[^^{mf) came from the 1/Z^ — 1 term, this would give Z^ — 1 = —0.003. Similarly, 
the valence mass dependence would give Za/ — 1 = 0.01. At present, we do not have sufficient 
data to measure Z^ — 1. without contamination from the c^ term. Thus, for the remainder of this 
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work, we take Z^ = 1 and we expect this to introduce an error of 1% or less in our axial current 

normalization. 

We also work with a local, four-dimensional axial current A^, which we renormalize so that 

ZaA^ = =c/^, so that 

ZAA^A';,{x) = 2mJl{x). (17) 

B. SU(2)l X SU(2)r Chiral Perturbation Theory for Kaon Physics 

1. Introduction 

An important goal of this paper is the comparison of the mass dependence of our lattice results 
with both SU(3)l x SU(3)r and SU(2)l x SU(2)r chiral perturbation theory (ChPT) at NLO {i.e. 
at one-loop order). For compactness of notation, throughout the remainder of this section we will 
simply refer to SU(3) and SU(2) ChPT. We will see that our data for pseudoscalar masses and 
decay constants only agrees with NLO SU(3) ChPT with good precision when the quark masses 
are small. This will be quantified in detail in Sec. VII, where we will see that NLO SU(3) ChPT 
provides a good description of our data when the average, dimensionless, input valence mass 
satisfies mavg < 0.01 (all our results will be obtained using data for valence u and d quark masses 
below this value). In order to attempt to extend the range of agreement to heavier masses (to 
include the strange quark for example), one possible approach is to use NNLO (or even higher 
order) ChPT. This introduces many additional low energy constants (LECs) and we find that we 
have insufficient data to determine all of these with sufficient precision. 

Instead of attempting to fit our data using NNLO (or even higher order) ChPT, we propose to use 
SU(2) ChPT at NLO. For pion physics, where it is possible to satisfy the condition, mavg < 0.01 
for the valence quarks, both SU(2) and SU(3) ChPT provide a good description of our data. 
Moreover, we will see in Sees. VI and VII that the SU(2) LECs obtained directly are consistent 
with those obtained by "converting" the SU(3) LECs to SU(2) [13, 26] (conversion formulae 
appear in Appendix B 2 c). For kaons on the other hand, the presence of the valence strange quark 
means that we do not satisfy the condition mavg < 0.01 and so we propose to use SU(2) ChPT at 
NLO. The effects of the strange quark mass are now absorbed into the LECs of the SU(2) ChPT. 
This eliminates the errors due to neglected higher powers of Zs/^csB prcsent when using SU(3) 
ChPT, but now one has to ensure that the u and d quark masses are sufficiently small to be able to 
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neglect higher order terms in Xud/Xs ■ Here Acsb ~ ^T^fit is the chiral symmetry breaking scale, 
while Xud and Xs are tree level masses of pseudoscalars made of ud and s quarks (see Appendix A). 
In this section we derive the NLO formulae for the behavior of the kaon's leptonic decay constant 
fK, its mass niK, and the bag parameter Bk as a function of the light-quark masses in SU(2) 
ChPT. At each stage, we start by presenting the arguments and calculations in the unitary theory, 
with valence and sea quark masses equal, and then proceed to the partially quenched theory. In 
the partially quenched theory in general the valence and sea strange quark masses, niy and m/, 
respectively, are also different and the SU(2) LECs depend on both these masses. For compactness 
of notation, we write these LECs with a single argument m/j, but it should be remembered that if 
niy 7^ m/j, then m;, should be replaced by the pair of variables niy, nih. We note that in this subsection 
quark masses written without a tilde have their usual interpretation as lagrangian mass parameters 
in continuum ChPT or PQChPT (a change from the notation defined in Appendix A and used 
elsewhere in this paper). 

We start by introducing our notation. In the unitary effective theory, the pion matrix, quark-mass 
matrix and the kaon fields are written in the form: 

(n'^I^Pi K+ \ (mi o\ (k+\ 

0= ,M=' \ ^ K=\ . (18) 

\ T^ -7tVV2j \0 mi) [k^ 

We work in the isospin limit in which the two light quarks are degenerate with mass m/. The pion 
matrices ^ and E are defined in the standard way: 

<^=exp(/0//) and 11 = ^^ (19) 

where the constant / is the pion decay constant (fjt) at lowest order in the chiral expansion (we 
use the convention in which fjt = 132MeV). Under global left- and right-handed transformations, 
L and R respectively, these quantities transform as follows: 

Z-^LZR^ ^^L^U"^ = U^R^ and K^UK, (20) 

where t/ is a function of L, R and the meson fields, but reduces to a global vector transformation 
when L = R. The pion Lagrangian at lowest order is 

4^i = 4 tr ^^Eo'^E^ + ^ tr(M^E + MEt) (21) 

o 4 

where / and B are the usual leading order LECs and, to this order, m| = 25m/. 
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The results for m| and /;^ at NLO in the chiral expansion are well known [12] and are presented 
in Eqs. (B34) and (B38) in the appendix. We now turn our attention to kaon physics using "Kaon 
ChPT" (KChPT). The corresponding chiral Lagrangian has already been introduced by Roessl [27] 
in order to study kK scattering close to threshold. At leading order the interaction of kaons with 
soft pions is described by the Lagrangian 

L^^l = D^K^D^K - M^K^K , (22) 

where the covariant derivative Z)^ is constructed using the vector field V^ so that D^j^K transforms 
like K under chiral transformations: 

D^K = d^K + V^K -^ UD^K . (23) 

Vn itself is constructed from the pion fields and transforms as 

V^ = ^(^^'^d^^ + ^d^^')^UV^U'^ + Ud^uK (24) 

We refer the reader to Eq. (1 1) of [27] for the higher order terms in the Kn Lagrangian (terms up 
to 0{mf) are explicitly listed). For completeness we also present the axial field: 

A, = ^ (r^^^ -^a^r) -. UA,U^. (25) 

In the following, when constructing Feynman diagrams from the Kn Lagrangian and the effective- 
theory local operators, we expand the vector and axial fields in terms of pion fields: 

V^ = [(l),d^(})]/2f + --- and A^ = -5^<^ // + ••-. (26) 

The use of Partially Quenched Chiral Perturbation Theory (PQChPT), in which we vary the va- 
lence quark masses independently of those of the sea quarks, gives us more scope to determine the 
low-energy constants of the unitary chiral theory and in this paper we will profit from this oppor- 
tunity. We consider the case with two degenerate light sea quarks, two degenerate light valence 
quarks and two 'ghost' degenerate light commuting quarks to cancel the loop effects of the valence 
quarks. Thus we are using the graded symmetry method [28, 29], based on Morel's ghost-quark 
trick [30]. See the lectures by Sharpe [31] for a recent review of PQChPT and references to the 
original literature (as well as many other applications of chiral perturbation theory to lattice QCD). 
We let the sea quarks have mass m/ while the valence and ghost quarks have mass nix. Thus we 
build an SU(4|2)l x SU(4|2)/; effective theory. We identify the quark flavors using the ordered list 
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where subscripts v and s are for valence and sea respectively, while the tilde denotes a ghost quark. 
We remind the reader that the strange quark mass does not appear explicitly here; rather, the LECs 
are functions of the strange quark mass. 
The matrix of Goldstone fields is written in block form 



*= . \ (28) 




where is a 4 x 4 block containing normal mesons, r/i = r/j is a 4 x 2 block of normal-ghost 
mesons and ^ is a 2 x 2 block of ghost-ghost mesons. The QCD pions are present in the central 
2x2 block of ^. As before, we work with the quantities ^ and Z 

<^=exp(/$//) and Z = ^\ (29) 

in terms of which the leading order PQ chiral Lagrangian is 

^p(2) = L str(D^E)'^D'^E + ^ str(^E^ + ILJl^) , (30) 

o 4 

where ^ is the mass matrix 

^ = diag(Mv,Mi,,Mv), with My = diag{mx,mjc) and M, = diag(m/,m/). (31) 
The quadratic terms in -S?pQ are 

-^paquad = \Md^^d^^) -5str(^$2) 



(32) 



5tr| ((^(^ + 771772) I " I -77277iM,-(^^M, 

The Lagrangian in Eq. (32) leads to a number of propagators which appear in the Feynman dia- 
grams in the following sections. We distinguish the propagators for the following mesons: 

• normal 'charged' (off-diagonal) mesons (~ qiqi)'- 

2^2 ^ m?2 = 5(mi+m2) = Ui+Z2)/2, (33) 

P -^n 

where mi 2 are valence or sea masses; 



17 



'charged' ghost mesons: (~ qiqi) 



p^ — m 



m?2 = 5(mi +m2) = Ul +X2)/2. 



(34) 
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Here mi, 2 are valence masses; 



'charged' quark-ghost mesons (~ qiq2 or ~ qiqi) 

2^2 ' '«?2=5('«l+'«2) = Ul+Z2)/2, 

jc^ — mj2 



(35) 



where mi 2 are a valence mass from the ghost and a valence or sea mass from the quark; 



'neutral' mesons from diagonal parts of and 0: 



i^^i^jj) = i 



. eA 



ij 



P^-Xi 



P^-Xi N{p^-Xi){p^-Xj)' 
for N sea quarks which are all degenerate with mass m/, with N — 2 here and 



(36) 



Si 



+ 1 valence or sea 



— 1 ghost 



(37) 



Note that the neutral propagators are determined after implementing the constraint str [4>] = 
0, which has the effect that neutral valence- valence meson propagators can have contribu- 
tions from neutral sea quark states (see [31] for example). 

We now add a kaon matter field 



K 



(K-] 




' SUv 


K^ 




sdv 


K+ 




SUs 


K^ 




sds 


K+ 




SUv 


i?0 , 




\Sdyj 



(38) 



The partially quenched Kn Lagrangian is the generalization of Eq. (22) written in terms of the 

partially quenched fields. 

In the following three subsections we derive the mass dependence of fx, m\ and Bk respectively. 
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2. Chiral Behavior of Jk 

The kaon decay constant, fx, is defined by the matrix element of the axial vector current: 

{Q\uy^y5s\K-{p)) = -I/kp^ . (39) 

We now need to match the QCD axial vector current to operators in KChPT. We start by doing this 
for left and right handed currents and then build the axial vector current. The left-handed QCD 
current for kaon decay is 

\qMi-r5)s (40) 

where q = u or d. It is convenient to promote c? to be a 2-component vector with components u 
and d and to introduce a constant 2-component spurion vector h in order to be able to project u or 
d as required; specifically we write the left-handed current as: 

2^%(l-r5)* (41) 

The current in Eq. (41) would be invariant under SU(2)l transformations if h transformed as 

h^Lh. 

We match to the effective theory by building quantities linear in h and with a single Lorentz index, 

which would also be invariant if h transformed as above. At lowest order, we identify two terms: 

(D^Ky^'^h, and K^A^^'^h. (42) 

By using equations of motion, operators with more covariant derivatives acting on the kaon field 

(which transform in the same way under chiral transformations) can be reduced to ones of higher 

order in the chiral expansion. 

For a right-handed current we simply take the /z-transformation to be /z ^ Rh and obtain the two 

operators 

{D^Ky^h, and K'^A^^h, (43) 

which transform in the same way as qhffiil + 75)5 . 

Under a parity transformation, the quantities in the effective theory transform as: 

K^-K, ^^^\ A^^-A^. (44) 
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Noting that parity transforms the left-handed current into the right-handed current and vice- versa, 
we deduce that the currents are of the form: 

Jfl = -Lai {D^Ky^h + ilAiK^A^^'^h (45) 

jI^Lai {D^K)'^^h + iLA2KU^^h, (46) 

where Lai and La2 are low-energy constants. To the extent that the kaon is regarded as heavy, the 
La2 term is subleading (and in any case does not contribute to fx at tree or one-loop level). The 
axial vector current J^ — 7^ is therefore 

jl=LAi{D^,K)\^+^'^)h + iLA2K'^A^{^-^'f)h. (47) 

We have introduced a factor of / in the La2 term to make the low energy constants real. This can 
be seen by considering charge conjugation ^. For the quark current (taking h to be real), we have 

J^^i = qhy^Yss ^ sy^Ysh^q = (-Z^)^- (48) 

In the effective theory, the charge conjugation transformations are 

^^^\ A^--Aj, y^--y; (49) 

together with 

Ka -^ Kl {D^K)a = d^Ka + {V^K)a ^ d^K^ - {K%)a = {D^K)l, (50) 

where a is a flavor label. Using these transformations, the current in the effective theory is trans- 
formed into its Hermitian conjugate under charge conjugation provided Lai and La2 are both real. 
The leading contribution to fx in the chiral expansion can be readily deduced. The first term in 
the axial current has a dp^K factor with no pion fields, so that at tree level /k is fixed by Lai, 
specifically 

/(^)K)=2Lai, (51) 

where nih is the mass of the strange quark (we are using the notation defined in Appendix A) . 
The NLO contribution to fx is obtained from the tadpole diagram in Fig. 1(a). From the Knn 
vertex in the Lai term in the axial current we obtain the contribution 

—iCf f d^k 1 —iCf 2 I '^n I 

2LAi{-ipp)^:^ J ^27t)H^-ml = ^^Ai{-iPn)j^^m^ I log ^ + constant 1 , (52) 
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where Cf = 3/4 is the eigenvalue of the quadratic Casimir operator for SU(2) in the fundamental 
representation and A;^ is an arbitrary renormalization scale. There is no contribution at one-loop 
order from the term in the axial current proportional to Lai and so we have shown that the chiral 
behavior of fx is of the form: 

A = /«K){l+cK)0-^£|j^log^}, (53) 

where c{mh) is a low-energy constant and m| = Xi ■ 

Before proceeding to discuss the chiral behavior of fx in partially quenched SU(2) ChPT, we 
compare the above calculation with that of fg in Heavy Meson ChPT [32]. In that case the NLO 
contribution has the corresponding contribution to Eq. (53), but in addition it has a second contri- 
bution from the self-energy diagram in Fig. 1(b). To discuss this for both /g and fx simultaneously, 
let P and P* be pseudoscalar and vector mesons containing a light quark (u or d) and a heavier 
antiquark (e.g. s or b). Let the masses of the P and P* mesons be M and M* respectively. The in- 
teraction of pions with pseudoscalar and vector mesons, P and P*, takes the form Mgd^nPP^/f, 
where the 1// arises because there is one pion and the Mg is put in for compatibility with heavy 
meson chiral perturbation theory (where the fields are usually normalized with an implicit factor 
of v^) so that the coupling constant g is dimensionless. The diagram in Fig. 1(b), for a P meson 
with momentum p, is proportional to {Mg/ f)^! where / is the Feynman integral: 

, f dH lH'[g^^-{p-l)^{p-l)y/Ml] 

We now compare the behavior of this integral with m| in the two cases: i) the heavy quark limit in 
which M ^ oo with M^ —M^ = 0{Aqqj^) which is an approximation for 5-physics and ii) M and 
M* both of O(Aqcd) and not degenerate, a situation which is appropriate for kaon physics. 
In the heavy quark limit, we write p = Mv + k, where v is the meson's four- velocity (v^ = 1) and 
the leading term in / in the l/M expansion is: 

2MJ (2;r)4(/2_m2)(v(/t-/)-A)' ^ ' 

where A = {Ml —M^)/2M which is of 0(Aq^j~)/M) and hence negligible compared to nij,; in the 
heavy quark limit. We stress that in this analysis we take the heavy-quark limit before the chiral 
limit, i.e. we keep Aqcd/^ ^ mn/AqcD ■ By power counting we see that the component of / 
contributing to the wavefunction renormalization can have a term proportional to 



21 

and so can lead to terms proportional to (^^//^)m|ln(m|/A^) in the behavior of /g. An explicit 
evaluation of the diagram confirms that such terms are indeed present [32]. 
In contrast, there is no contribution of the form m|logm| in KChPT, i.e. when we take the chiral 
limit of small m| while keeping M and M^ fixed, non-degenerate and neither especially large or 
small. To see this we combine the denominators in / using Feynman parametrization and write 

^ /■! f dH N{a) 

1= da -—-j—^-^, (57) 

Jo J (2;r)4D2(a)' 

where a is the Feynman parameter, 

Ma)-(/ + arf-'" + "'"-"„."'""'"' (58) 

Mi 

and 

D{a) = l^ + a{l- a)p^ -{I -a) ml- aMl . (59) 

The question we are addressing is whether there is a contribution to / of the form m|logm|. If 
there is such a term then we can isolate it by differentiating with respect to m|, setting m| -^ 
and searching for a logarithmic (infrared) divergence in 

dl ^ f\ r. ^ f d^l N{a) 

—^ =-2 da(l-a) -—-j—^, — V — ^- (60) 

dml ^2^^ Jo ^ 'J {2nYD\a,ml = Q) 

We use dimensional regularization in D = 4 + 2e dimensions. The contribution from the first term 
inA^inEq. (58)is 

dh 



dml 



r , . f d^l /^+«v 

Jo ^^ ^'J {2K)D[f + a{\-a)p^-aMl]^ 



/^ Ja (1 - a)a^[M2 - (1-a)/]^ I (2+£)r(-e) - ^^^-^^} 



(4;r)2+e 

(61) 

The contribution to the self-energy is the coefficient of p^ at p^ = M^, which we can isolate by 

differentiating with respect to p^ and setting p^ = M^ . By inspection we can readily verify that 

there are no infrared singular terms and hence no term proportional to m^ logm^ in / . Note that 

M^ > M and hence M^ — (1 — (x)p^ > for p^ = M^ throughout the integration region in a. A 

parallel argument shows that there is also no contribution to / of the form m|logm| from the 

second term in //(a) on the right-hand side of Eq. (58) . 

Note also that / does contain a term proportional to m^logm|. In this case we seek an infra-red 

divergence after differentiating / twice with respect to m| and setting m^ ^ . The power of 
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l^ + a{l — a)p^ — OlMI in the integrand is —4 as compared to —3 in Eq. (61). After performing 

the /-integration, we find that for small a the integrand ~ a^+^, which diverges in 4 dimensions 

and is a signature of the presence of a logm| term {i.e. a m^ logm| term in /). 

There are also no chiral logarithms from the tadpole diagram in Fig. 1(c). The KKkk term in L)^f. 

in Eq. (22) is proportional to K{dp^K) [(^,d^(^] and so the integrand is an odd function of the pion's 

momentum and the integral vanishes. 

Finally in this section we derive the behavior of fK in the partially quenched case, i.e. in 

PQKChPT. The axial vector current is of the form of that in Eq. (47), but uses the PQ K and ^ 

fields. As before, the spurion h is chosen to pull out the appropriate flavor. The tree and quadratic 

terms from this current are 

jl- = 2LAid^Kl(l-^ +■■■ (62) 

V ^J J ba 

where we have dropped terms with a single derivative acting on a Goldstone field. 
To calculate the decay constant for a kaon containing a valence light quark, we take the flavor 
index a = 1 or 2. Choosing, for illustration, a = 1 we find 

1 

There are also terms containing <?^^2' ^^^ these have Goldstone field pairings which cannot be 
contracted to give a one-loop tadpole contribution and thus are not listed above. The one-loop 
tadpole contribution to /^+ is 

_A 1/^!L/ 2/ / 2/ / 1 i^-xi \ .... 

The 5 terms in the integrand of Eq. (64) have the following origin: 

1. The first term comes from the quark-ghost mesons in the tadpole loop. These have valence- 
valence quark masses and the minus sign arises from the closed loop of anticommuting 
fields. 

2. The second term has the 0i202i propagator with both quarks having valence masses. 

3. The third term has the ^13^31 + ^14^41 propagators with one quark having the valence-quark 
mass and the other the sea-quark mass. 

4&5. The final two terms have a neutral propagator in the tadpole diagram from the 0^j term on 
the rhs of Eq. (63). Both quarks have valence-quark masses. 



^LAid^K^ 



1 - ;T72 (^l'11^2,ll + ^1,12^2,21 + ^H + ^12^21 + ^13^31 + 014^41; 



(63) 
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We observe, as expected, that the tadpole contribution with the valence-ghost propagator cancels 
that with the valence- valence one. Extracting the chiral logarithms from the loop integrals we 
arrive at the expression: 



fxh = 2Lai < 1 - 



1 



{AnfY 



Xt ' Xx ^ Xi ' Xx . Xi ^Xx 1 Xx 



+ • 



(65) 



The analytic terms can be proportional to Xi and Xx and so we obtain the final result for the mass 
behavior of a meson with a light valence quark x and a heavier (strange) valence quark h: 



fxh = f{mh)l\ + 



4-2 ^^ ~^ fl '^^ 



Xx~T~ Xi 1 Xxi Xi . Xi ^Xx 1 Xx 



(66) 



where A3 ,4 are LEC's and we remind the reader that if the valence and sea strange quark masses 
are different then the LECs p'^' and A3,4 depend on both these masses. Eq. (66) agrees with the 
corresponding calculation of the chiral behavior of fg [32], when the terms proportional to the 
square of the BB*k coupling are neglected. The right-hand side of Eq. (66) reduces to Eq. (53) in 
the unitary limit Xx = Xi ■ 



3. Chiral Behavior of rri^ 

We now consider the chiral behavior of mj., starting in the unitary theory in which the valence and 
sea masses are equal. In principle the chiral logarithms could come from the tadpole diagram in 
Fig. 1(c). However, just as for the wavefunction renormalization, there is no such contribution at 
one-loop order. The KKkk term in L^,.^ in Eq. (22) is proportional to K{d^K) [0, f?'^^] and so the 
integrand is an odd function of the pion's momentum and the integral vanishes. This is also the 
case in the partially quenched theory and so, at NLO in the chiral expansion, the mass-dependence 
of mj^ comes from the analytical terms coming from the higher order terms in the Kn Lagrangian. 
For the partially quenched theory at NLO 



m 



xh 



= 5(^)(m;,)mJl + 



r2 ^' ~'~ f2 ^^ 



(67) 



In the unitary theory we have the natural simplification of only a single low-energy constant. 
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4. Chiral Behavior ofBfc 

The non-perturbative strong interaction effects in neutral kaon mixing are contained in the matrix 
elements of the QCD operator 

SLYndLSLf'dL, (68) 

between K^ and K'^ states. This operator is part of a multiplet transforming under SU(2)l on the 
two down quarks and is symmetric under the interchange of the two dis. Hence it transforms as 
a triplet under SU(2)l and is a singlet under SU(2)r (the analogous situation for neutral 5-meson 
mixing is that the corresponding operator is part of a 6 representation of SU(3)l [33]). 
In our effective kaon theory the combination ^K transforms as ^K ^ L^K under SU(2)l, and we 
can use it to build the leading A5 = 2 operator, 

^ab = 2li{^K){,{^K),}, (69) 

which is symmetrised in the flavor indices a and b and has a LEG /3 (the analogous construction 
for neutral 5-meson mixing in SU(3)l x SU(3)r heavy meson chiral perturbation theory leads to 
an operator 4-j5Q{^P^^'^)[a{^P^^')},j [33], where P^^' and P^^' destroy mesons containing heavy 
Q and Q quarks respectively). As in the discussion of the axial current, operators with one or both 
of the kaon fields on the right hand side of Eq. (69) replaced by D"K, where D is the covariant 
derivative (the Lorentz indices have been suppressed) transform in the the same way as ^ab- Again, 
by using the equations of motion, the leading component of these operators can be reduced to i^ab ■ 
In order to evaluate Bk we will need the matrix element {K'^\&22\K^) and thus need to use the 
K^K^ piece of ffji- At tree level we find 

4j8 = \mlflBK. (70) 

In order to evaluate the one-loop contributions to the matrix element (see Fig. 2) we expand ^22 
up to second order in the pion fields: 

--l^KaKt-'^U^K)a{Wb+\[{<^^K)aKj,+Ka{fK)h\\+---. (71) 

We observe, by comparing to the expansion of the axial vector current in Eq. (47), that the two 
terms containing 0^ lead to relative one-loop corrections each of which is the same as the relative 
one-loop correction for f/^^. Since we are calculating a matrix element proportional to f^Bx, these 
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two corrections will not affect Bk and we need only calculate the one-loop correction from the 
{^K)a{^K)i, term. Thus the relevant component of ^22 is the K^K^ term in 

-^^i<^K)2mh = -'^'^K'K' + - . ill) 

Contracting the two pions into a loop leads to a contribution 

-J^J Jl^^W^ml = -7^(4^-.logX| + -^lyticterms. (73) 



11 9. m? 



Hence we find that 

This agrees with the result for 5g in Eq. (3.8) of Sharpe and Zhang [32] when the BB*7C coupling 

We now promote the above discussion to the partially quenched case. The effective theory operator 
is still of the form of Eq. (69), but now the flavor labels take 6 values and ^ is expanded in terms 
of 4> in Eq. (28). This echoes the discussion for partially-quenched neutral 5-meson mixing in 
Sharpe and Zhang [32], starting from the QCD discussion in [33]. 

To obtain the mixing matrix element for neutral kaons with valence down quarks, we look at the 
K2K2 piece of the ^22 operator, expanded up to second order in the Goldstone fields: 

^22 = 215 K2K2 - ^ I (*i^)2(*i^)2 + ^ [{^^K)2K2+K2{^^K)2] I + ■ ■ ■ . (75) 

As in the unquenched calculation above, the (4>^^)2^2 terms lead to chiral logarithms which are 
cancelled by those from /^ when extracting Bk from the A5 = 2 matrix element. Hence the pieces 
of the partially-quenched operator we need are: 

2l5(l-jj^i2]K2K2. (76) 

The ^22 field is a valence- valence neutral meson, so we find a one-loop correction 

_4^ fJ^S^ 1 l^-Xi \_ 2/3 X. . .™. 

f^ J {2nr\l^-X. 2{f-x.f] (4;r/)2^'°^A|+"" ^"^ 

For Bjc the result is: 
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This agrees with Eq. (3.9) in Sharpe and Zhang [32] with ^ ^ and Nf = 2. It also reduces to the 
QCD result above, Eq. (74), when Xx = Xi = ^|- Noting that there are analytic terms proportional 
to both Xi and Xx the final result for the mass dependence of Bk is that in Eq. (B48) of the appendix 

D dW^ \i^ , M^^ ^biimh) Xl ,_ Xx \ .„3. 

B^h = 5^s ("^h) < 1 + -J^Xi + -j^Xx - y2n^ ^ A^ f ' ^ ^ 

where bi^2 are LECs. 

5. Comments 

We conclude this section with the observation that for all the physical quantities considered in this 
section the chiral logarithms in the SU(2) theory can be simply deduced from those in the SU(3) 
theory. In all these cases the chiral behavior is of the generic form: 

O = Olo ( 1 + chiral logarithms + analytic terms) , (79) 

where O is the physical quantity (pseudoscalar mass-squared, decay constant or bag parameter) 
and there is a single LEG at lowest order (Olo)- The generic form in Eq. (79) holds in both the 
SU(2) and SU(3) theories. As an example consider the pseudoscalar decay constant /^/,, with a 
degenerate heavy (strange) valence and sea quark with mass m/, and partially quenched up and 
down quarks. In the partially quenched SU(3) theory at NLO we have from Eq. (B18): 
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Xrj Xx ^ Xx Xt] Xh Xr\^ . 
"Converting" to the SU(2) theory, we can expand in m^/nifi and nii/nih obtaining: 

+ analytic and higher order terms 



(81) 
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The analytic terms now include contributions proportional to Xh which do not vanish in the SU(2) 
chiral limit. They are absorbed into the lowest order SU(2) LEC /, which now depends on m/,. 
Whether or not the relation between / and /o is well approximated by one-loop SU(3) ChPT or 
whether higher order {Xh ^^'^ higher powers) terms must be included depends on m;, and on the 
convergence of the series. In any case, \i f^h is to satisfy the generic form dictated by SU(2) ChPT 
given in Eq. (79) the chiral logs are fixed from Eq. (81): 



r Y* j 1 A*-^ ~^ f^l 1 f^^ ' A*/ I /^ ,, ,, \ 1 A--^ 



+ analytic and higher order terms 



(82) 



which agrees with the result from the direct evaluation in the SU(2) theory (see Eq. (66)). The 
same is true for the other quantities being studied here. Of course, in general there may be more 
than one operator in the effective theory with the same quantum numbers as the QCD operator 
whose matrix element is to be evaluated. In such cases there are more than one LEC at leading 
order, Eq. (79) does not apply and the simple arguments presented here have to be generalized. 

III. SIMULATION DETAILS AND ENSEMBLE PROPERTIES 

Following the work in [16, 34] we have used the Iwasaki gauge action and the domain wall fermion 
action. By producing smoother gauge fields at the lattice scale, for a fixed lattice spacing in 
physical units, the Iwasaki action removes some of the gauge field dislocations that contribute to 
the residual chiral symmetry breaking for domain wall fermions at finite Ls. While suppressing 
such dislocations improves residual chiral symmetry breaking, it also suppresses topology change 
in the evolution [17]. Since we want our ensembles to sample topological sectors of the theory as 
well as possible, we have found the Iwasaki gauge action to provide a reasonable balance between 
these two, contradictory goals. 

We generate ensembles with two degenerate light quarks, whose bare input mass is given by m/, 
and one heavy strange quark of mass nih. We use the exact Rational Hybrid Monte Carlo (RHMC) 
algorithm to generate ensembles. During the course of this work, we made improvements to the 
RHMC algorithm, as detailed in [34], yielding three variants of the RHMC algorithm, 0, I and 
II. The original RHMC algorithm [35, 36] was used for the 2+1 flavor simulations in [16] and 
the RHMC I was used for most of the simulations in [34]. We also compared the RHMC I and 
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II algorithms in [34]. We find the RHMC II algorithm to be the best version to date. It uses a 
single stochastic noise source to estimate ratios of determinants, which reduces the size of the 
forces in the molecular dynamics integration [15]. It produces ensembles that change topology 
more rapidly than RHMC I, likely due to this decrease in the fermionic force. It uses a multiple 
time scale Omelyan integrator and light quark preconditioning [37, 38], which has the effect of 
making the time spent solving the light quark Dirac equation less than the time spent solving for 
heavier quarks. 

We can briefly define the RHMC II algorithm. Letting ^(m,) = Dp^p(M5,m,)DDWF(^5,w/), 
where Ddwf(^57 "^i) is the domain wall fermion Dirac operator, M5 is the domain wall height and 
m, is the mass of the quark we wish to simulate, we write the fermionic contribution to the path 
integral as 

3 

(83) 
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(Recall that for domain wall fermions, a regulator must be added to remove the bulk infinity that 
would arise as L^ -^ °°. The ^(1) terms in Eq. (83) are this regulator.) Each ratio of determinants 
on the right-hand side of Eq. (83) is represented by a single stochastic estimator, i.e. a single 
estimator is used for the ratio involving m/ and m/, and three estimators are used for the three, frac- 
tional power determinants involving m/, and 1 . The light quark ratio is integrated on the coarsest 
time scale, using a conventional Hybrid Monte Carlo algorithm. The three 1/2 powers of ratios 
involving the heavy dynamical quark are integrated on a 2 times finer time scale using the RHMC 
algorithm. For each integration step of the RHMC and for each field, a single solution of the Dirac 
equation is required with mass nih and two solutions with the regulator mass of 1. (The solutions 
are found with the conjugate gradient algorithm.) The regulator mass solutions take roughly half 
as many conjugate gradient iterations in the RHMC as the m^ solutions, but there are twice as 
many required. Thus the time spent solving the Dirac equation in the two cases is comparable. 
Because of the three different stochastic estimators needed in the RHMC, and the coarser time 
step for the HMC, the number of conjugate gradient iterations in the RHMC are larger than in 
the HMC part of the algorithm for the quark masses currently used. This has made simulating at 
lighter quark masses much less expensive. The gauge field is integrated on an even finer time scale 
than the RHMC. For a further review of the RHMC, see [39]. 

As mentioned, we had already begun evolving configurations with the RHMC algorithm as 
improved versions were being developed. Table I gives the molecular dynamics time when we 
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changed from one algorithm to the next, for each of our ensembles. Since, the RHMC II gives more 
decorrelated lattices than earlier versions, all of the observables in this paper have been measured 
only on configurations generated by the RHMC II algorithm. The configurations generated with 
the earlier algorithms are being used only for equilibration. 

As one indicator of how well the RHMC II algorithm is decorrelating our ensembles, we have mea- 
sured the global topological charge. Of course, for infinite volumes, global topological charge is 
not relevant to local physics. However, in finite volume knowing that global topology is changing 
gives us evidence that there are local topological fluctuations. To calculate the topological charge 
the configurations are first cooled by applying 30 updates with a quenched, 5 loop improved gauge 
action with zero coupling strength, /3 = oo. After cooling, the topological charge is measured using 
a 5 loop improved gluonic topological charge operator, the 5Li method of [40] and the results are 
shown in Fig. 3. This figure shows that the RHMC algorithm is sampling different topological 
sectors for each of our ensembles and the histograms indicate, by their symmetry and shape, that 
the topological landscape has been sampled reasonably well. 

Figures 4 and 5 show the integrated autocorrelation time for the pion correlator at a separation of 
12 lattice spacings for the m/ = 0.005 ensemble. Both show an integrated autocorrelation time 
for this quantity of 10 to 15 molecular dynamics time units. Figure 6 shows the evolution of the 
plaquette for the m/ = 0.02 ensemble. For this short-distance observable, equilibration took a few 
hundred molecular dynamics time units. 

IV. LATTICE RESULTS FOR PSEUDOSCALAR MASSES AND DECAY CONSTANTS 

A. Measurements 

Since we are most interested in the light quark limit where pion masses are close to the physical 
values, we put significant computational effort into the measurements on the two lightest ensem- 
bles with nil = 0.005 and m/ = 0.01. On these two ensembles, three separate measurements have 
been done, which we shall refer to as Full-Partially-Quenched (FPQ), Degenerate (DEG) and 
Unitary (UNI). In the FPQ measurement we computed hadron two-point correlators for all the 
combinations with valence quark masses m^, niy E {0.001, 0.005, 0.01, 0.02, 0.03, 0.04}. The 
DEG dataset consists of hadron correlators with degenerate quarks with the same list of masses. 
In the UNI dataset only quark propagators with masses equal to the light and strange sea quark 
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masses were calculated, and the light-light, light-strange and strange-strange meson correlators 
were then constructed. On the m/ = 0.02 and m/ = 0.03 ensembles, only the DEG and UNI cal- 
culations were performed. Details of the measurements, including the gauge configurations used, 
separation of each measurement and the total number of measurements, can be found in Table II. 
For the m/ = 0.005 and m/ = 0.01 ensembles, we blocked the data so that each block contains 
measurements from every 80 molecular-dynamics time units, while for the m; = 0.02 and 0.03 
ensembles the block size was chosen to be 40 time units to yield a reasonable number of jackknife 
blocks to perform the analysis. 

In each of the FPQ measurements, we used Coulomb gauge fixed wall (W) sources of size 24^, 
which were placed at t = 5 and t = 59. For each source, two quark propagators were calculated, 
one with the periodic, and the other anti-periodic, boundary condition in the temporal direction. 
The sum of these two quark propagators (as a single quark propagator) was then used to construct 
the meson correlators. The resulting cancellation of backward propagating states has the effect of 
doubling the temporal extent of the lattice, so over much of the lattice volume, there is no excited 
state contribution to the hadron propagator. We find this works well with our Coulomb gauge fixed 
wall sources, which generally have small statistical errors but are not tuned to remove excited state 
contaminations. The long plateaus and small statistical errors allow us to work far enough from 
the source that excited states are not a worry. 

The DEG measurements used a Coulomb gauge fixed box source (B) of size 16-' which we have 
found to give the optimal early onset of the plateaus for pseudoscalar mesons. The sources were 
placed at two timeslices, t = and 32. 

In the UNI measurements, the propagators were calculated from Coulomb gauge fixed hydrogen 
S-wavef unction (H) sources [41], with radius r = 3.5, in lattice units, or gauge invariant gaussian 
(G) sources and sinks with radius r = 4. Again the sources were placed at multiple timeslices for 
a better sampling over the gauge fields. 

We also used two different interpolating operators for the pseudoscalar state, namely, the pseu- 
doscalar operator P"{x) = q{x)Y5{T" /2)q{x), which we refer to using a short-hand notation P, and 
the axialvector operator A" (;c) = q{x)Y^Y5{'^"/^)li^)^ which we refer to as A. Here t" is a flavor 
symmetry generator. Unless otherwise specified, we only consider the time component of A^{x), 
i.e., A'^{x). Together with the different source/sink smearings described above, this allows us to 
construct several different pseudoscalar meson correlators, which are tabulated in Table III. (Ta- 
ble IV has similar information for vector and tensor correlators, which are discussed in detail in 
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Section IX.) The notation in Table 111 follows Ref.[16], where quark propagators are specified 
by the smearings applied in the source and sink. For example, WL means a quark propagator 
calculated with a wall (W) source and a local (L) sink. Meson correlators are then denoted by the 
the type of quark propagators used in the contraction, e.g., WL-WL refers to a meson correlator 
with two WL quark propagators. In this case we may also use WL to refer to the meson correlator 
unless ambiguities arise. 

In this paper we only consider meson correlators with zero momentum projection. Our meson 
states are normalized such that the time dependence of the correlators in the large t limit can be 

expressed as 

(0\n'M7r\(7r\ni^\0\ r .^ ..i 

(84) 






where the superscripts specify the smearings for the quark propagators, and the subscripts specify 

the interpolating operators used, nixy is the ground-state mass of a pseudoscalar meson composed 

of two valence quarks with masses nix and niy. For convenience the amplitude of the correlator 

will be denoted as 

_ (0|a;'|;r)(;r|Q-^|0) 

^0.02 = ^;;;^ • (85) 

B. Residual Chiral Symmetry Breaking 

As discussed in Section II A, the finite size in the fifth dimension of the domain wall fermion 
formulation in our current simulations permits the mixing between the two light quark states bound 
to the boundary walls, resulting in a residual chiral symmetry breaking [7, 14, 17]. Close to the 
continuum limit, this effect can be quantified as a residual mass term (denoted as m„J, which, up 
to 0{a^), behaves as a regular quark mass, making the total quark mass effectively the sum of the 
bare input quark mass and the residual mass. We determine the residual mass by computing the 
ratio [14, 16]: 

. <^^^^(^')^'°.°'). (86, 

where 7| is the "mid-point" operator which mixes the quark states from the left and right 
walls [7, 14], and P" is the pseudoscalar operator defined earlier. For t large enough that only 
pions contribute to the correlators in Eq. (86), R{t) -^ m[^^{mf). From the DEG data sets, we ob- 
tained the residual mass for each pair of valence quark masses by averaging the plateaus oi R{t) 
from t = \Otot = 32. Figure 7 shows the residual mass as a function of the valence quark mass. 
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We can see that the residual mass has a linear dependence on the input quark masses, and the 
magnitude of the slope is larger in the valence sector than the sea sector. The slopes have opposite 
signs, as was also observed in our previous simulations on smaller volumes [16, 34]. In particular, 
d\n{m[^) / din'?^ = —4.1 for the ensemble with m/ = 0.005 and —4.5 for the m/ = 0.03 ensemble. 
We also find d\n{m[^^ / drn^r^ = 6 and, for the dynamical quark mass dependence in the unitary 
case, the result is d\n{m[^^) / dnii = 1.1. In terms of the discussion at the end of Section II A, 
it appears that the Z^ — 1 and dimension five terms are comparable in magnitude. We define a 
mass -independent residual mass for our simulations by evaluating m^,.,{mi) at the m/ = limit. The 
straight line in Fig. 7 shows a linear fit to the four unitary points with nix — my — mi. And the star 
is the resulting residual mass at m/ = 0, which we determined to be 

m„, = 0.00315(2). (87) 

It is worth noting that this result is in decent agreement with that obtained from earlier simulations 
at the same gauge coupling but with a smaller lattice volume [34], indicating that the finite volume 
effect for the residual mass is small. 

C. Pseudoscalar Meson Masses and Decay Constants 

1. Results from the nii = 0.005 and 0.01 ensembles 

To obtain the masses and decay constants for the pseudoscalar mesons from the m/ = 0.005 and 
nil = 0.01 ensembles, we focus on the FPQ data sets, since the long plateaus of the FPQ data sets 
lessen the uncertainty in choosing the fit ranges. In fact, including the UNI and DEG measure- 
ments in the simultaneous fits (described below) does not change significantly either the central 
values or statistical errors of the final results. Thus here we only quote results from the FPQ data 
sets. 

For each mass combination in the FPQ data sets, we performed a simultaneous fit to the five 
correlators as specified in the second column of Table III, to determine a common mass and a 
different amplitude for each correlator. This approach has the advantage of averaging out the 
systematic uncertainties due to different characteristics of the interpolating operators (P or A), and 
reducing the uncertainty of fit range choices. One drawback is that, as the number of data points 
included in the fit (~ 220) is greater than the number of independent measurements we have in 
hand, the covariance matrix is poorly resolved, and we lean on uncorrelated fits to obtain results. 
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The resulting pseudoscalar meson masses for all the mass combinations of the m/ = 0.005 and 0.01 
ensembles are given in Table V. Note that the quoted ;t^^/dof 's only give a qualitative indication of 
how different fits compare to each other. They are not meaningful quantitative indications of the 
goodness of the fits. 

Besides a common mass for each valence mass combination, we also obtained an amplitude for 
each correlator included in the simultaneous fit, namely 

^^,^p7,^^^,^pr and y^/^''. (88) 

Combinations of these amplitudes can be used to determine the pseudoscalar decay constant, 

which we will describe in the following. 

The pseudoscalar decay constant can be determined by 

u^Mami, (89) 

where Za is the axial vector current renormalization constant which relates the (partially) con- 
served axial vector current in the original five-dimensional theory of domain wall fermions to the 
local four-dimensional axial vector current [4]. Alternatively, the axial Ward identity in Eq. (11) 
connects the divergence of the partially conserved axial vector current to the pseudoscalar density, 
allowing us to determine the pseudosclar decay constant through 

_ {m, + my)\{n\P\0)\ 

Jxy — 2 ' ^^^^ 

where m^^^ = m^_y -|- m„,. 

We determined Za from the DEG data sets using the improved ratio [4] of (^(?)i'(0)) to 
(A4(?)P(0)), where £/^{x) is the partially conserved axial vector current for domain wall 
fermions [4, 7]. Similar to the residual mass, we computed the value of Z^ at each unitary quark 
mass, and extrapolated it to the chiral limit at m/ = — m,^„ finding 

Za = 0.7161(1), (91) 

which agrees very well with the result of our smaller- volume simulations [34]. 
There are in principle five methods [42] to determine the matrix element necessary for the fxy 
calculation in Eq. (89) or (90). These use jackknife ratios of different amplitudes in Eq. (88). 
We found that different methods give statistically consistent results for all quark masses except 
for the lightest or heaviest quark masses, where systematic deviations for results obtained using 
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Eq. (89) and Eq. (90) are observed. The methods where the pseudoscalar density is used make 
use of the equation of motion for the conserved axial current and are more indirect, since we 
have not systematically controlled higher order renormalization effects here. As was discussed 
in Section II A, the renormalization of the conserved axial current is known to 0{m,^,), and we 
know the renormalization between the local and conserved current to high accuracy. Therefore, 
we have chosen as our preferred result a determination of f^^y using axial current matrix elements 
and postpone further investigation of the differences we are seeing for future work. We also find 
that the statisical errors are the smallest when only the axial vector current matrix elements are 
used in the analysis. Thus, we determine fxy by 



2 ^LW^ 

The final results for the m/ = 0.005 and 0.01 ensembles are shown in Table VI. 
Comparing the results from the m/ = 0.01 ensemble with those obtained from a smaller volume 
reported in Ref.[34], we find that there is a statistically significant 4% difference at the unitary 
point of nij^ = niy = nii = 0.01, as shown in Fig. 8. As we shall see in Section VIA, such a 
difference is about twice as large as predicted by finite volume chiral perturbation theory. While it 
is possible that ChPT still underestimates the finite volume effect for fxy, we want to point out that 
the 16^ X 32 gauge configurations were generated using RHMC I which could intrinsically induce 
long-range autocorrelations. Another difference is that there we determined f^y using correlators 
with much shorter plateaus compared to the 24-^ x 64 ensembles discussed here, which means the 
smaller- volume results may suffer from larger systematic errors. Although we tried our best to 
estimate uncertainties on the physical quantities of the 16^ x 32 simulations, the errors might still 
be underestimated. With these factors taken into account, the discrepancy between the lattice data 
and the predictions of the finite volume ChPT may not be as significant as it appears. 

2. Results from the mi = 0.02 and 0.03 ensembles 

For the two heavy ensembles with m/ = 0.02 and 0.03 we performed uncorrected simultaneous fits 
to four correlators from the DEC data sets: ^if (0, ^pp (0. ^af (0 and ^//(O- The time slices 
included in the fits are from ? = 10 to 32 for each correlator. The pseudoscalar meson masses 
obtained this way are presented in Table VII. The corresponding amplitudes of these correlators 
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can be used to detennine the pseudoscalar meson decay constants by the following relation: 
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^ = ^-^— ^^^W^' (93) 
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the results of which are given in Table VIII. 



V. RESULTS FOR THE OMEGA BARYON MASS 

We have chosen to set the lattice scale using the Q. baryon mass because it is made of three 

strange valence quarks and so does not have non-analytic light quark mass terms at NLO in chiral 

perturbation theory [43, 44]. Thus, it is justified to linearly extrapolate in m/ to obtain the physical 

value. 

The interpolating field for the decuplet baryons is 

Zdec = ^abc\q^aCyyiqb\qc, (94) 

which couples to both spin 3/2 and 1/2 states, the former being the ground state [45]. For q = s, 
the strange quark, the operator destroys an ^ baryon. 

The correlation function is constructed using the BL quark propagator described in Section IV 
computed from sources set at time slices and 32 for every twentieth trajectory, starting from a 
thermalized one. Results are averaged over these sources and blocked in bins of size 80 trajectories 
to reduce auto-correlations. The total number of measurements made varies between about 100 
and 200, depending on the light sea quark mass (see Table IX). To avail ourselves of all possible 
statistics, we average over forward propagating particle states, and backward propagating anti- 
particle states, as well as all three spatial directions for the interpolating operator. 
The average correlation function is fit using a jackknife procedure where the covariance matrix is 
estimated for each jackknife block. The minimum time slice used in the fit was varied in the range 
4 < tmin < 10; the maximum was fixed to 14. Though the mass values do not depend sensitively 
on the choice of tmim we chose to take as central values the masses obtained from tmin = 8 to 
minimize excited state contamination without losing control of the statistical error and because x^ 
is acceptable for this tmin value. The baryon mass values for valence quark masses m^ = 0.03 and 
0.04 are summarized in Table IX and plotted in Fig. 9. 

Also shown in the figure is a linear extrapolation, for each degenerate valence mass, to the physical 
value of the light sea quark mass (solid lines). The linear fit to the nix = 0.03 data has a ;t^/d.o.f. 



36 

of 1.2 and the m^ = 0.04 data has a j^/d.o.f. of 1.4. In the next section we use these values to 
interpolate to the physical strange quark mass in order to set the lattice scale using the Q. baryon 
mass. 

A systematic error in the Q. mass, and hence in the lattice spacing, results from the fixed unphys- 
ical value of the strange sea quark mass in the simulation. This effect can be estimated from the 
slope of the baryon mass with respect to the light sea quark mass (see Fig. 9). The effect, which 
is similar for both valence quark masses, is about 2-3 % for a (required) shift in sea quark mass 
of about 0.008. Since this results from varying both light quarks together, it should be divided by 
two for the single strange quark in the simulation. This is about a one percent downward shift in 
the mass, which is smaller than the statistical errors, as shown in the next section, so we ignore it. 

VI. LATTICE SCALE, QUARK MASSES AND DECAY CONSTANTS USING SU(2) CHPT 

A. SU(2) X SU(2) Chiral Fits 

In this subsection we present a brief overview of the techniques we use to perform the chiral 
extrapolations for pseudoscalar masses and decay constants. For our simulations with 2+1 flavors 
of dynamical fermions, the traditional approach would be to use the SU(3) chiral perturbation 
theory to guide the extrapolations. However, as we shall see in Section VII A, NLO SU(3) ChPT 
does not describe our data for the pseudoscalar mesons which contain a quark as heavy as the 
strange quark. In order to extrapolate our results reliably to the physical quark masses (and in 
particular to the physical strange quark mass) we therefore use a different approach and impose 
only SU(2) chiral symmetry in the light quark sector, as has been discussed in Section II. We 
divide our discussion into two parts. In Section VI A 1 we discuss the chiral extrapolations in the 
pion sector and in Section VI A 2 we extend the discussion to kaons. We present physical results 
for the quark masses and pseudoscalar decay constants using the SU(2) chiral fits in Section VI B 
below. 

1. The Pion Sector 

In the pion sector, the chiral dynamics of SU(2)l x SU(2)ij symmetry is sufficient to describe the 
physics we are interested in. The effects of the dynamical strange quark are fully contained in the 
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low-energy constants of SU(2)/^ x SU(2)^ ChPT [12]. There is a caveat. A posteriori we learn 
that the dynamical strange quark mass used in our simulations is a little larger than the physical 
one. Ideally, we would need to perform the simulations with several different dynamical strange 
quark masses and interpolate the results to the physical strange quark mass. In the absence of such 
simulations, the dynamical strange quark mass is fixed to m/, = 0.04, and we need to be aware that 
the corresponding LECs will have a corresponding small systematic error (see Sec. VIC 3). 
The next-to-leading order partially quenched SU(2) formulae for the squared masses, mL, and de- 
cay constants, fxy, of the light pseudoscalar mesons composed of two valence quarks with masses 
nix and niy can be derived from [46], as has been done, for instance, in [47]. For reference, we 
summarize these formulae in Appendix B 2. As discussed in Section II, for domain wall fermions 
the effective quark mass is the sum of the input quark parameters m^jih and the residual mass, 
which we denote as in^j^i^h = '^x.yj,h + "^res (here / and h denote the light and heavy sea quarks 
respectively). To next-to-leading order, this is the only correction to the continuum chiral formu- 
lae. Since fxy and m^ share two unknown low energy constants B and /, we perform combined 
fits to both of them, using the two ensembles with dynamical light quark masses of m/ = 0.005 
and 0.01. Such combined fits with valence quark masses m^, niy = 0.001,0.005, and 0.01 (which 
therefore satisfy the cut mavg = {nix + niy)/! < 0.01) are shown in Fig. 10. The values of the fit 
parameters are given in Table X, where we quote the scale-dependent low energy constants at two 
commonly used renormalization scales A;^ = 770 MeV and 1 .0 GeV. These partially quenched 
LECs (of SU(4|2) PQChPT) can be related to those in the physical unitary theory (SU(2)), l^ 
and l^, by Eqs. (B35) and (B39). In Table X we also quote the values of It, and I4, which are 
conventionally defined at the scale of the pion mass, see Eq. (B44). 

Using partially quenched NLO SU(2) chiral perturbation theory we are able to describe the data 
with the mass cut mavg = {mx + %)/2 < nicui = 0.01. It is not possible however, to extend the 
range of masses significantly. For example, the quality of the fit degrades significantly if instead 
we impose a mass cut of mavg < 0.02: the ;^^/d.o.f increases from 0.3 {nicut = 0.01) to 1.2 (mcut = 
0.015) and even 3.6 for nicut = 0.02. (We only report the uncorrelated ;^^/d.o.f. since we could 
not reliably determine the correlations, so these numbers most likely underestimate the "true" 
deviation from the fit.) This suggests that higher orders of chiral perturbation theory (NNLO or 
even more) are needed in this mass range. (See also the discussion in Sec. VIC 1 regarding the 
inclusion of analytic NNLO terms.) 
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2. The Kaon Sector 



In order to incorporate the valence strange quark mass in the chiral extrapolations, we need to 
extend the SU(2) ChPT to the kaon sector. As has been discussed in Section II, this can be done 
by considering the kaon as a heavy meson under the assumption that m/ <^ m^, in analogy to the 
heavy-light meson (or heavy baryon) chiral perturbation theory [32, 48]. In this framework the 
chiral expansions are ordered in powers of m//m/, or equivalently m^/mj., as well as m|/Ay. We 
present the NLO chiral formulae for m^ and fxy, the masses and decay constants of pseudoscalar 
mesons composed of a light valence quark with mass nij^ and a heavy valence quark with mass niy 
in Eqs. (66) and (67) (see also App. B 3 for the complete set of formulae at NLO). 
For kaons a new set of low-energy constants f^^', B^^> and Ai, 2,3,4 are introduced. These LECs 
depend on the valence and dynamical strange quark masses and we have made this explicit by 
assigning a functional dependence, writing for instance, /(^-'(m/,). Here m^ represents both the 
masses of the heavy (strange) dynamical and valence quarks. The parameters B and / are the usual 
LECs in the (pionic) SU(2) theory, and their values are determined from the quantities in the pion 
sector, i.e., their values are fixed to the ones given in Table X. 

The strange quark mass dependence in the formulae for the meson masses and decay constants, 
Eqs. (66) and (67) are contained in the low energy constants. Since it is not possible to tune the 
dynamical heavy quark mass m/, perfectly to the mass m, of the physical strange quark, in principle 
we would wish to perform the calculations at several values of m/, and interpolate the results to 
m.s. Since we only have results at a single value of nih {nih = 0.04) we rely as far as possible on 
using the behavour of our results with the valence strange quark masses, niy, to extrapolate to the 
physical kaon. 

Finally, we briefly summarise our fitting procedure. We start by fitting the mass dependence of 
m^ and fj^, using Eqs. (66) and (67) for nix ^ [0.001 , 0.01] and niy = 0.04, with B and / fixed to the 
values in Table X. Since these two formulae do not have any unknown parameters in common, the 
fits can be done independently. We then extrapolate to the physical value of the light quark mass to 
obtain mK{my = 0.04) and fximy = 0.04) . (The physical light quark mass niud is determined from 
the SU(2) pion sector fits, see VI B 1 for a detailed description.) This step is shown in Fig. 11, 
where the circles give the extrapolated values mfc{my = 0.04) (right) and fximy = 0.04) (left). 
Then similar fits are performed with niy = 0.03, giving values for niximy = 0.03) and fK{tny ~ 
0.03). Eventually, we take the results at niy = 0.03 and 0.04 and interpolate linearly between them. 
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As discussed in Sec. VI B, we define the physical strange quark mass from the interpolation of mjy 
to the physical value of mj.. Having determined the physical value of the strange quark mass, we 
extract the kaon decay constant from the interpolation offj^y to this mass. The linear fits are shown 
in Fig. 12. In the same figure we also include the points obtained with niy = 0.02. Although at such 
a low mass the use of SU(2) chiral perturbation theory is questionable, we see that nevertheless the 
points lie very close to the line obtained from the interpolation between niy = 0.04 and niy = 0.03, 
suggesting that the higher-order corrections of 0{mi/fhh)^ may be small. 



B. Lattice Scale, Physical Quark Masses and Decay Constants 

Having established the procedure for the chiral extrapolations, we now turn to obtaining the phys- 
ical results. In this subsection we start by explaining how the lattice scale and physical bare quark 
masses were determined. We then use the renormalization constant for the quark masses which 
was determined non-perturbatively in [21] to obtain the physical quark masses in the MS scheme at 
/i = 2GeV. Finally we present the results for /;;:, fx and the CKM matrix element \Vus\ determined 
from our measured ratio fx/ f% together with experimental inputs. 

1. Determination ofniud, tn^ and a^^ 

In order to determine the quark masses m„^, nis and the lattice spacing a, we need to compare our 
lattice results for three quantities to their physical values. We take the pseudoscalar masses m^ 
and nix as two of these. Natural choices for the third would be the p mass or the Sommer scale. 
However, since the p meson is a resonance with a finite width and there is an uncertainty of order 
10% in the value of the Sommer scale, we choose instead to use the mass of the Q. baryon. This is 
a state composed of three valence strange quarks. One advantage of using this baryon mass, is that 
up to NLO in ;^PT it is free of logarithms containing the light quark masses [43, 44]. Therefore the 
extrapolation of the measured masses to the light physical mass can be readily performed using a 
linear ansatz without an uncertainty due to chiral logarithms. 

The physical quark masses were obtained from the SU(2) x SU(2) fits described in Section VIA. 
For the average light quark mass, m„^ = (m„ + m^)/1 , we solved for a pion mass of niji = 
135.0MeV, corresponding to the physical neutral pion mass, while for the strange quark mass 
nis the fit to the kaon mass was solved at niK = 495.7 MeV, which is the quadratically averaged 
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neutral and charged kaon mass. 

The determination of the three parameters is a coupled problem. The lattice scale is needed to 
convert masses between lattice units and physical ones, whereas the quark masses are needed for 
the extrapolation and interpolation in the light and strange quark masses to obtain the mass of 
the Q. baryon. We performed the determination iteratively, starting with an initial "guess" for the 
quark masses, fixing the lattice spacing by requiring that m^^ takes its physical value, then using 
this value of a to adjust the quark masses by requiring that niji and nix take their physical values 
and so on until no further changes in the parameters were observed. The final values for I /a, a, 
niud, ms can be found in Table XI, with only the statistical errors shown. (We will discuss the 
systematic errors in Sec. VIC.) We also find the quark mass ratio to be m„j : fhs = 1 : 28.8(4). 
An independent check of the lattice scale has been done from the spectrum of heavy-heavy and 
heavy-strange quark states using the non-perturbatively determined relativistic heavy quark ac- 
tion [49, 50]. The preliminary analysis [51] gave l/a= 1.749(14) GeV, where the error is statisti- 
cal only. The agreement of this result with that obtained from the niQ^- shown in Table XI suggests 
that the associated systematic errors are small. 

2. Quark Masses in MS Scheme 

In this section we present the results for the quark masses in the MS renormalization scheme. The 
renormalization factor Z^ = I/Z5 (where S represents the scalar density) needed to convert the 
bare quark masses to the commonly used MS scheme at a scale of 2 GeV has been calculated with 
the same action on a 16^ x 32 x 16 lattice [21]. (For details about this ensemble of configurations, 
see [34].) We first matched the bare lattice operators to the RI-MOM scheme using the non- 
perturbative Rome-Southampton technique [52], and then performed a perturbative matching to 
the MS scheme. The relation between the renormalized and bare quark masses is given by 

mf^(2 GeV) = Z^(2GeV) a'^ (amx) , (95) 

where X = ud, s and for clarity we made explicit the factors of the lattice spacing a. In Ref. [21] it 
was found that Z^^ (2 GeV) = 1.656(48) (150), where the first error is statistical, while the second 
is the combined systematic uncertainty due to residual chiral symmetry breaking, the use of a 
three-loop matching factor between the RI-MOM and MS schemes and mass effects. With this 
renormalization constant, the results for the average light quark mass and the strange quark mass 
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are determined to be 

m^(2GeV) = 3.72(16) (33) MeV, (96) 

mp(2GeV) = 107.3(4.4)(9.7)MeV, (97) 

where the first error is the combined statistical error from the bare lattice quark masses, the lattice 
spacing and Z^^ and the second error is the systematic error in Z^^. 

Related to the renormalization of the quark masses is the renormalization of the lowest order LEG 
B, which is proportional to the quark condensate and therefore is renormalized by Zs = l/Z^, 
ensuring that the product of B and a quark mass does not depend on the renormalization scheme 
or scale. We have 

5^(2GeV) = fz^{2GcY)y^ ■ (l/a) ■ aB (98) 

= 2.52(11) (23) GeV, (99) 

where the first error is statistical and the second is the systematic uncertainty from renormalization. 
The renormalized value for the chiral condensate, E = f^B/2, reads 

E^(2GeV) = (^255(8) (8) Mev) . (100) 

•?■ Ik, Ik and \Vus\ 

The chiral fits described in Section VI A also allow us to determine the decay constants /;j and 
/k- Extrapolation to the physical light quark mass niud using the SU(2)xSU(2) fits in the pion 
sector gives /;,; = 124. 1(3.6) MeV, while interpolation to the physical strange quark mass m^, as 
described in Section VIA, gives fx = 149.6(3.6) MeV. In both cases only the statistical errors 
are given. Systematic errors will be discussed in Section VI G. These results can be compared to 
the physical values [53], fn =130.7(0. 1)(0.36) and /k =159.8(1.4)(0.44) MeV. For the ratio of the 
decay constants we find fx / fn = 1.205(18) compared to the experimental value of 1.223(12). See 
[102] for possible implications of an updated result for the measured value of the GKM matrix 
element y„d [54]. 

An important application of the result for /k/Ik is the determination of the GKM matrix ele- 
ment \Vus\, as has been pointed out in [55]. Using the experimentally determined branching ratios 
r(^ -^ IJ.v{y)) and r(;r -^ ij.v{y)), together with the radiative electroweak corrections from [53], 
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we obtain: 

\Vus\/\Vud\ = 0.2292(034),tat(118)syst(005)other, (101) 

where the first two errors are statistical and systematic (as discussed in Sec. VIC below), while 
the third error is the combined uncertainty from the measurement of the branching ratios and the 
radiative electroweak corrections [53]. Taking \Vud\ = 0.97377(27) from super-allowed nuclear 
j8 -decays [53], we obtain 

|K,| = 0.2232(033)stat(115)syst(005)other (102) 

(The third error now also contains the uncertainty in \Vud\)- This implies 

\Vus\^+\Vud\^ = 0.9980(15)stat(51)syst(06)other = 0.9980(54)total- (103) 

Since \Vub\ is negligible, our result implies that the unitarity relation is satisfied within the uncer- 
tainties, constraining the possible breaking of quark-lepton universality in models of new physics. 
For completeness we now present the corresponding results obtained using the recently published 
value \Vud\ = 0.97418(26) [54] (instead of the PDG value of \Vud\ = 0.97377(27) which was used 
in obtaining the results in Eq. (102) and Eq. (103) above): 

|K,| = 0.2232(033)stat(115)syst(005)other, (104) 

\Vus\^ + \Vud\^ = 0.9989(15)stat(51)syst(06)other = 0.9989(54)total. (105) 

We have also determined Vus (with a smaller error) from the experimentally measured rates of 
semileptonic K£2 decays, by computing the form factor f^{q^ = 0), where q is the momentum 
transfer. The results (in particular IV^^I = 0.2249(14)) and an outline of the calculation are pre- 
sented in [18]. 

C. Systematic Errors 

In this section we will discuss the systematic errors resulting from various sources: 1. chiral 
extrapolations; 2. finite volume effects; 3. the unphysical dynamical strange quark mass and 4. 
scaling errors. 
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1. Errors arising from the chiral extrapolations 

As explained above, we obtain our results by using NLO chiral perturbation theory to fit the 
measured values of the meson masses and decay constants in the range of valence quark masses 
satisfying mavg < 0.01 and with m/ = 0.005 and 0.01. We do not have enough data to extend the 
analysis fully to NNLO where the chiral logarithms depend on the Gasser-Leutwyler coefficients 
and there are too many parameters to make the fits meaningful. In order to estimate the effects of 
the (neglected) higher order terms in the chiral expansion we have therefore extended the range of 
quark masses to mavg < 0.02 but included only the analytic NNLO-terms in the fit functions. In 
the pion sector, the analytic terms at NNLO are quadratic in the ;i^'s and symmetric under Xx'^ Xy\ 
there are four such terms: 

Ux + Z>■)^ Ux-Z>■)^ Ux + Z.y)Z/ and %}■ (106) 

We therefore have four new parameters for the behavour of the mass of the pseudoscalar mesons 
and four for the decay constants. Unfortunately, with only two values for the dynamical sea quark 
mass (which is proportional to ;^/), we were not able to resolve the complete sea quark mass 
dependence up to NNLO in our fits. We therefore dropped the term proportional to j^ and were 
then able to obtain stable fits with a good ;^^/d.o.f. (about 0.1) . We stress that the NNLO analytic 
terms are necessary in order to obtain good fits; NLO chiral perturbation theory does not represent 
our data in the extended range of masses mavg < 0.02. With the NNLO analytic terms included, 
not only do we obtain good fits to our data but the behavour of the chiral expansion is sensible, 
i.e. the relative importance of successive terms in the series is as expected. We illustrate this point 
in Fig. 14, where, in the top two graphs, we plot the size of LO, NLO, and NNLO contributions 
(normalized such that the LO contribution equals one) for the masses-squared and decay constants 
of mesons with nix = niy = mi {i.e. for mesons in the unitary theory). Up to a quark mass of 
0.01 (the 3rd vertical dashed line in the plots, counting from the left) the NLO-contributions are 
much larger than the NNLO terms, which are negligible in this region. As expected, the NNLO 
corrections become more important for larger quark masses and in particular, for masses close to 
0.02 (rightmost vertical dashed line) they represent a significant fraction of the decay constant. 
This study is reassuring, but, since we have neglected the chiral logarithms at NNLO terms, it is 
only a phenomenological estimate of the likely contributions from higher order terms in the chiral 
expansion. 
As our final estimate of the errors due to the chiral extrapolation we take twice the difference 
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between the results from our standard NLO SU(2)-ChPT fit in the mass range m^^g < 0.01 and 
those from including the analytic NNLO terms in the range mgyg < 0.02. The factor of 2 was 
included, rather cautiously, because of the phenomenological nature of the NNLO analysis, and in 
particular because we were unable to study the sea-quark mass dependence. The systematic error 
is reported in Table XII. 

The systematic error due to the extrapolation in the mass of the light quark in the kaon quantities, 
nis and fx, was also estimated as twice the difference between our standard results and those 
obtained using kaon SU(2) PQChPT in the mass range m_^ < 0.02 with the values of / and B 
obtained from the NNLO fit described above. The errors are also tabulated in Table XII. For 
this case we get a good description of our data up to m^ = 0.02 without the necessity of adding 
NNLO analytic terms (x^, xf, and XxXi)- Therefore one expects those terms to have a negligible 
effect even in the mass range nix < 0.02. The lower panels in Fig. 14 show the LO and NLO 
contributions, normalized so that the LO contribution equals one, to mjf^ and //;,. The horizontal 
axis is Xh which is proportional to the total light sea quark mass. We see that the size of the NLO 
contributions are at most 35% for the range of masses where we have data. 

2. Finite volume effects 

In this section we estimate the errors due to the fact that the simulations are performed on a finite 
volume, (aLy ^ (2.74fm)-^. We do this by following the procedure proposed by Gasser and 
Leutwyler [56-58] in which one compares results obtained in ChPT using meson propagators in 
finite volume, G^^'{x), to those obtained with the infinite- volume propagators, G{x). The finite- 
volume and infinite- volume propagators are related by 

G^^\x) = G{x) + Y,Gi^ + Ln), (107) 

with the 3-vector n El? . This allows one to calculate the corrections in the PQChPT formulae 
for the meson masses and decay constants directly. We list the expressions in Appendix C, where 
details of the numerical implementation can also be found. (See also refs. [59-61] for the anal- 
ogous discussion with staggered PQChPT.) We now repeat the chiral fits using PQChPT with the 
finite-volume propagators and obtain a new set of fit parameters. For each physical quantity, we 
assign the difference between the results obtained using the finite-volume and infinite-volume for- 
mulae as the systematic error and tabulate these finite- volume errors in Table XII for both the pion 
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and the kaon sectors. In the latter case we find that the finite-volume corrections in m^ and fK are 
actually negligible. In Fig. 13 we plot the correction factor for the decay constants and squared 
masses, ( 1 — A^/ ) and ( 1 — A^'" ) (cf . App. C for their definition), respectively, for our two values 
of the dynamical light quark mass (m/ = 0.01 and 0.005). Since the volume is reasonably large, 
these corrections are found to be below 1 percent, except for the very light valence masses. 
We now compare our results for the finite- volume corrections obtained using ChPT (as described 
above) to the theoretical predictions by Colangelo et al. [62], who use a resummed Luscher for- 
mula. The results in ref. [62] are presented for the unitary theory {m^ = my = nii in our notation) 
and so we restrict the comparison to this case. We take the results in Tabs. 3 and 4 of [62] which 
we interpolate to the volume and pion masses used in our simulation. In Table XIII we present the 
finite-volume corrections for the pion mass and decay constant estimated using 3 methods, SU(2) 
and SU(3) ChPT and that of Ref. [62]. The quantities Rm and Rf are defined by 

mu{L)-mu{oo) 1 ,„2 ML) - fuH tf ,,p,g. 

Rm = ^—. = 7:^11 7 Rf = T~r~\ ^^n ' ^^^^^ 

Within the uncertainties, we see a reasonable agreement between the estimates obtained with the 
different methods. 

In an earlier publication [34], we presented results for the masses and decay constants obtained 
with the same action and coupling as in this paper, but on a smaller volume, about (1.83fm)^ 
compared to (2.74fm)^ used here. (In [34] a different method was used to extract the lattice 
scale. Since the gauge action and coupling are identical, for consistency we quote here the volume 
obtained with the value for a^^ used throughout this work.). 

Specifically the lattice in [34] has 16 points in each spatial direction as compared to 24 in this 
paper. We are therefore able to compare the theoretical estimates of the finite- volume corrections 
described above, with the difference of our results on the two lattices with different volumes. We 
compare the measured pion mass and decay constant at the common unitary point at m^ = niy ~ 
nil = 0.01, nifi = 0.04 for which the pion mass is ^ 427 MeV (this is the smallest mass for which 
we have results on both lattices). On the smaller volume, both finite-volume corrected fits (SU(2) 
and SU(3)) give results for the correction factors comparable to those interpolated from Colangelo 
et al.. To obtain these correction factors, we used the parameters obtained from the fits to the 24^ 
data and evaluated the finite-volume corrections for L = 16 (instead of the original L = 24). The 
SU(2) fits resulted in 7?1|^''^ = 0.42(04)% and R^^^"^ = 1.66(14)%, where the superscript is meant 
to indicate that this is the finite- volume correction for a spatial volume of L^ = 16^. In Table XIV 
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we list the measured values of the meson masses and decay constants on the 24^ and 16^ lattices 
and compare their ratios to the prediction from finite volume ChPT and the resummed Liischer 
formula (interpolated from [62]). We obtain this prediction from the ratio of the finite-volume 
correction factor for the 16^ spatial volume and that for the 24-^ spatial volume, 

(109) 

The superscripts in Eq.(109) indicate the spatial volumes on which the results for niji and f-^ were 
obtained or for which the finite-volume corrections Rm and Rf were evaluated. The question 
mark reminds us that we are checking whether the measured values of the ratios are equal to 
the theoretical predictions for the finite- volume effects. Whereas (within the errors) the observed 
(2.0 ±1.2)% effect for the pion mass is somewhat better reproduced by the FV ChPT predictions 
than the (3.7 ± 1.1)% for the decay constant, in both cases all three method tend to underestimate 
the observed FV effects. We make the following two comments. First, the precision of FV ChPT 
at NLO and the resummed Liischer predictions for the finite-volume corrections are expected to 
improve with increasing volume. A possible reason for underestimating the difference between 
the results from the 16^ and 24^ lattices may be that the smaller volume is already borderline for 
the methods used here. Whereas the values for nijiL would appear to be sufficient in both cases; 
on the 16^ lattice mjiL ■^ 4.0 compared to ?^ 5.8 on the 24^ lattice, it may be that the smaller lattice 
does satisfy the relation L ^ (Vlfji)^^ ~ Ifm sufficiently. The second comment is that, as has 
already been pointed out in Sec. IV, the errors on the measured quantities on the smaller volume in 
general, and on the decay constants in particular, may be underestimated due to shorter plateaus. 
Furthermore, in this present paper we have obtained the decay constants using an improved ansatz 
(see also Sec. IV). For these reasons, we believe that the finite volume effects on the larger lattice 
are well described by our approach and give a reliable estimate for the systematic error. 

3. Effects of the unphysical dynamical strange quark mass 

In this subsection we estimate the systematic error due to the fact that the mass of the dynamical 
heavy quark (m;,) turns out, a posteriori, to be about 15% larger that of the physical strange quark 
{nis). To perform this estimate we exploit the fact that the LECs of SU(3) ChPT are independent 
of the quark masses, and use the NLO conversion formulae relating SU(3) and SU(2) LECs [13] 
given in Eqs. (B40-B43), together with the results for the SU(3) LECs presented in Sec. VII A. In 
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this way we obtain the shifts in the SU(2) LECs: 
5K)-5K) = 5o|^(2Lf -Lf))(z.-Z.) -:^(zaog^-Z/Jog^) 



Jo ^ ^ '72;r Jq ^ -^^..;|, -^^■■;f' 

(110) 

/K)-/K) ^ /o{|4^^(^-^.) -^(zaog^-z.iog^)}, (111) 

r,{m,)-r,{m,) = ^log^, (112) 

l\{ms)-l\{mh) = TT^log — • (113) 

64;r^ Xs 

(Note: here one has to use Bq to determine Xx-) We use the shifted values of the SU(2) LECs 
to determine m„^ and f^i and again use the shift in those two quantities as the estimate of the 
corresponding systematic uncertainties. All these uncertainties are presented in Table XII and are 

(2) 

small compared to the total errors. (We did not quantify the systematic error on the separate L^- , 
just on the /s 4 since these are the phenomenologically relevant observables.) 
For the kaon sector the shift due to changing the dynamical heavy quark mass to its physical value 
can be obtained from Eqs. (B49-B54). We perform the conversion for fixed niy at niy = 0.03 and 
0.04 separately and then in the same way as before linearly interpolate the heavy valence quark 
mass niy to m^ to determine the shift in the kaon mass and the kaon decay constant. A shifted 
value for the physical strange quark mass is obtained by solving the following system. Find a 
value of m.v such that the interpolation in niy to niy = nis between m^^ Jnis, niy) between niy = 0.03 
and 0.04 gives the squared physical kaon mass. Here m^^ (nis.niy) was obtained from the unitary 
extrapolation in the light quark mass to ni^d using shifted LECs. The shifted LECs are found by 
starting from their values at the dynamical heavy quark mass used in the simulation and changing 
to their values at the new value of m^ (at fixed niy). By this method, we observe a 2% shift in m^ 
and less than 0.5% shift in /k- Therefore, conservatively we will quote the systematic error from 
tnh 7^ ms to be 2% and 1% for fhs and fx, respectively. 



4. Scaling errors 

At present we only have data at a single value for the lattice spacing (a = 0. 1 1 fm) and so can- 
not perform a scaling study to extrapolate our results to the continuum limit. (We are currently 
generating an ensemble of configurations on a finer lattice and so will soon be able to study the 
discretization errors directly.) We therefore take our central values and assign to them a systematic 
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uncertainty of 4% due to the missing continuum extrapolation. The 4% is simply an estimate of 
(aAQco)^ with our value of a^^ = 1 .73 GeV, but we stress that it is only when we will have results 
at more than one value of the lattice spacing that we will be able to quantify this error reliably. 
Table XII contains the 4% estimate for scaling error for the measured quantities. Again, we did 

(2) 

not calculate this uncertainty for the (phenomenologically uninteresting) LECs L] , but only for 
/3^4. Here one has to keep in mind, that a relative error affects the universal low energy scales 

A3 4 = A;^-exp(^^/^4 ) = (139MeV)-exp(/3 4/2), (114) 

V 73,4 ' / 

meaning that a 4% uncertainty translates into an absolute error of 0.08 for /3,4. 

D. Comparison of our results with other recent determinations. 

In this subsection we compare our results for the SU(2) LECs to those obtained in the continuum 
[63, 64] and in other lattice simulations with either Nf = 2+\ [65, 66] or Nf = 2 [67, 68] dynam- 
ical fermions. We also compare the values we obtain for the quark masses with those from other 
recent simulations [66, 69-73]. 

Our value for the decay constant in the SU(2) chiral limit, / = 114.8(4.1)stat(8.1)syst MeV, is 
consistent, within our uncertainties, with the phenomenological estimates of 122.3(0.5) [63] or 
121.9(0.7) [74]. For the ratio fjt/f we obtain 1.080(08) (with only the statistical error quoted, 
since the numerator and denominator are likely to be highly correlated), which agrees very well 
with the phenomenological results 1.069(04) [63] and 1.072(07) [74]. Previous lattice simulations 
give, e.g., 1.052(2)(+^) [66] or 1.075 [68] for the ratio and for / the resuks read / = 124.2MeV 
[66] and 121.6MeV [68], respectively, where in Ref. [68] no estimate of the error was provided. 
From MILC's results for / and the chiral condensate, (E^^(2GeV)) '3 — 
278(l)stat(5)ren(l3)systMeV [66] (compare to 255(8)stat(8)ren(13)systMeV obtained in this 
paper), one can via E = f^B/2 deduce their resuh for the LO-LEC B^^{2GeY) = 2.79 GeV (no 
estimate for the error). Comparing this number with our value of 2.52(. 1 l)stat(-23)ren(- 12)syst GeV 
shows approximately the same relative deviation as for the decay constant in the chiral limit, /. 
In Table XV we compare the values for the NLO LECs li,4 (we also include the values obtained 
by converting the SU(3) LECs, see Sec. VII C). Within the quoted uncertainties all these numbers 
agree, except the result for l^ as quoted in [75], which used the MILC SU(3)-LECs as input 
obtained from a "NLO plus analytic NNLO and beyond" fit. A recent update of the MILC results 
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finds a more consistent value for h obtained from a pure NLO fit. The two results in the table for 
/4 from the CERN Collaboration, where one was obtained by including NNLO analytic terms, not 
only agree with our result for the central value, but also suggest a comparable magnitude for the 
systematic error (as estimated from the difference between the NLO and partial NNLO fits). 
Quark masses are computed in lattice simulations using a variety of actions. We end this section 
by making the observation that results obtained using non-perturbative renormalization (such as 
those in this paper) appear to be generally higher than those obtained by renormalizing the mass 
perturbatively (mostly using 2-loops). Whether or not this is significant or merely a coincidence 
is still to be investigated. We illustrate this point by tabulating recent results obtained using non- 
perturbative renormalization [69-71] and perturbative renormalization [66, 72, 73] in Table XVI. 
Note that for the simulations in this table, both different fermion actions and different numbers of 
dynamical fermions were used. The ratios nis/niud are also tabulated and agree better across the 
different computations. 

VII. FITTING MASSES AND DECAY CONSTANTS TO NLO SU(3) CHPT 

In the previous section, we have fit our lattice data to NLO SU(2) ChPT formulae and found good 
agreement. Additionally, for quark masses where we have data, we find that for m^^, the NLO 
terms are at most a 10% correction, while for /// they are less than 30%, as shown in Fig. 14. Our 
estimates for NNLO effects in this range are appropriately small, < 5%, leading to the conclusion 
that SU(2) ChPT converges reasonably well here. This implies that our range of light quark masses 
is small enough that NLO SU(2) ChPT provides more than just a useful, smooth phenomenological 
function to fit our data to, rather it represents the theoretically correct description of our data. 
We now turn to fitting our data to SU(3) ChPT, with the following two points in mind: 1) we want 
to determine how well our data is fit by the NLO SU(3) ChPT formulae and 2) if the fits agree 
with our data, how convincing is the convergence of the (now known) SU(3) ChPT series for the 
quark masses where we have data. We will see that the answer to the first point is yes, if we only 
include observables involving light quarks, and is no for observables including the strange quark. 
For the second point, even in the light quark case, we find the convergence of the series to be poor. 
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A. SU(3) X SU(3) Chiral Fits 

At NLO, the 2+1 flavor, partially quenched chiral formulae [46] for the squared masses, m^ , and 
decay constants, f^y, of the light pseudoscalar mesons composed of two non-degenerate quarks 
with masses m^ and niy involve six unknown low energy constants (LECs), which we denote as 
Bq, /o and ^4 5 g g- (The complete formulae used in our chiral fits are summarized in Appendix 
B 1.) As fxy and m^ share the unknown low energy constants Bq and /o, we performed combined 
fits to both of them, using the two ensembles with dynamical light quark masses of m/ = 0.005 and 
0.01 . A reasonable ;^^/d.o.f. of 0.7 (uncorrelated) could only be achieved by imposing a cut in the 
average valence quark mass of mavg = {m^ + my)/! < 0.01, corresponding to partially quenched 
pion masses in the range of about 250 MeV to 420 MeV. The chiral fits with such a cut are shown 
in Fig. 15, and the fit parameters are given in Table XVII, where, for convenience, we quote the 
scale-dependent LECs at two commonly used chiral scales of A;^ = 770 MeV and 1 GeV Only 
statistical errors are quoted in this table, as discussed in more detail below. In these fits, the valence 
quark masses are all in a region where NLO SU(3) ChPT might be expected to be reasonably 
reliable (corresponding to masses below 420 MeV for pseudoscalar mesons made of such valence 
quarks), but we point out that the dynamical heavy quark mass, which is approximately 15% 
higher than the physical strange quark mass, lies outside of our fit region (mavg < 0.01). Even the 
combination of our lightest valence quark mass of 0.001 with the dynamical heavy quark mass 
leads to a pseudoscalar mass of '^ 554 MeV. 

Extending the fit range to valence quarks satisfying mavg = {mx + yny)/! < 0.03, or a partially 
quenched valence pseudoscalar mass of up to 660 MeV, we find that the fit curves miss almost all 
of the data points, resulting in a poor ;^^/d.o.f of 5.7. This is shown in Fig. 16. Similar behavior 
is seen for a lower mass cut of mavg = {tn^ + niy)/! < 0.02, which allows the largest valence 
pseudoscalar masses to be just above the kaon mass, and gives an uncorrelated ;^^/d.o.f of about 
3. From these fits, we conclude that we cannot use SU(3) ChPT at NLO to obtain physical results 
for kaon observables, such as /a: or the physical strange quark mass from niK, with well-controlled 
interpolation or extrapolation errors. For these important quantities, we have used kaon SU(2) 
ChPT as described in previous sections. 
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B. Systematic Errors 

For the SU(3) PQChPT results, we have to consider the same sources of systematic uncertainties 
(chiral extrapolation errors, finite volume effects and scaling errors) as in the SU(2) case, except 
that, in principle, the SU(3) approach allows one to directly extrapolate the dynamical heavy quark 
mass to the value of the physical strange quark mass. 

For the chiral extrapolations, we first investigate the systematic error from the convergence of NLO 
ChPT. We have already seen that we cannot use NLO SU(3) ChPT for valence quarks near the 
physical strange quark mass, since the fits do not agree well with our data. Focusing on the results 
from fits to the light quark region (mavg = (m^ + m^ ) /2 < 0.0 1 ), where our data is well represented 
by the NLO formula, in Fig. 17 we plot the LO and NLO contributions to the pseudoscalar decay 
constant /// in the unitary case. The values of /// are plotted as functions of the light quark mass 
parameter (Xi °^ ^/) ^^^ the heavy quark mass parameter (Xh °^ f^h) using our results from the 
SU(3) fit with the mass cut mavg < 0.01. We see that the corrections are generally large, even for 
relatively small light-quark masses. For example, the top left-hand panel shows that for m/, = 0.04 
and for pseudoscalar masses in the 250-400MeV range, i.e. in the region where we have data, the 
NLO correction is as large as 70% of the LO term. 

In the preceding paragraph we found that, for the unitary case with nij,; ^ 400 MeV, the SU(3) 
ChPT expansion for the decay constant is approximately of the form /o(l +0.7 + 0{p^)), which 
leads one to expect the O(p^) term to be about (0.7)^ ^ 0.5. This would be a 30% correction to 
the sum of LO plus NLO terms. However, the NLO fits agree with our data quite well, certainly 
excluding corrections beyond the few percent level. From this we conclude that the agreement 
between our data and the SU(3) ChPT formula is not indicative of NLO SU(3) ChPT working well 
for these quark masses, since this would imply anomalously small 0{p'^) corrections. It seems 
that the SU(3) ChPT formula are serving as smooth fitting functions with sufficient flexibility 
in the parameters to absorb the effects of higher order terms and to match our data. We would 
need more data to investigate this hypothesis fully and, in particular, this data should come from 
ensembles with more than the single m/, we currently have. This will be done in the future and 
if our hypothesis is correct then, as discussed in the following section, part of the reason for the 
small value of /o is due to the fact that it absorbs some of the NNLO corrections. We now briefly 
discuss our limited attempts to use our current data to obtain additional information but stress that 
a complete analysis of the range of validity and precision of SU(3) ChPT will have to wait until 
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we have data at more values of quark masses. 

One can also try to extend the range of validity of ;^PT by going from NLO to NNLO. The 
complete continuum NNLO formulae are available in the literature [76] and we have done some 
preliminary fits of our data, augmented with results from an ensemble with m/ = 0.02 [77], to these 

(3) 

formulae. Many more LECs are needed at NNLO (L- for / = to 9 and 12 linear combinations 
of Ki) and it is not currently clear how stable such fits will be. We defer a discussion of these fits 
pending completion of the ongoing analysis. 

To further probe the convergence of the series, we have dropped the additional logarithm terms 
that appear at NNLO and fitted to the analytic terms, as we did for our SU(2) fits in Section VI A. 
The analytic terms for the meson masses and decay constants for SU(3) are 

ixx + Xyf , {Xx - Xyf , {Xx + Xy)x , z^ z^ = {^xf + X^)/^ ■ (115) 

With our present data, limited to a single value of Xh and to a small set of ;f/'s, we were not 
able to include the last two terms in our NNLO phenomenological fits. The fits led, e.g., to an 
approximate 10% increase in the value of /o (for this particular case a fit range of mavg < 0.01 has 
been used) and to even more significant changes in some of the NLO LECs. This supports our 
conjecture that SU(3) ChPT shows a slow convergence and NNLO terms are indeed important. 
On the other hand, since NNLO-terms are not negligible, taking into account only (some of) the 
analytic NNLO terms and neglecting the logarithmic terms is not sufficient to determine the chiral 
behaviour of observables quantitatively. It is for this reason that we choose to use SU(2) ChPT to 
determine our physical results and at this time we do not quote any estimates for systematic errors 
for quantities from our fits with SU(3) ChPT. 

In spite of our reservations about the convergence of SU(3) ChPT, we have used the corresponding 
formulae in finite volumes to estimate finite volume effects (see Appendix C), as we had done 
previously for the SU(2) case. We find similar results for the correction factors. Table XIII also 
contains the corrections obtained in the SU(3) case for the unitary pions; for more details see the 
discussion in Sec. VI C 2. 

C. Comparing SU(3) x SU(3) and SU(2) x SU(2) Chiral Fits 

In this section, we compare the results of our SU(3) ChPT fits with other results and also compare 
them with the results of our SU(2) fits. This serves to further probe the behavior of ChPT in the 
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region of quark masses we are studying. We stress however, that we believe that the convergence 
of the SU(3) series is relatively poor and therefore at this stage any quantitative conclusions will 
be limited. 

To compare the SU(3) fit results with the previous ones obtained in SU(2) PQChPT (Sec. VIA), 
we use Eqs. (B40)-(B43) (cf. also [13, 26]) to match the three flavor ChPT to the two flavor 
case at NLO. The results for Bq, /q and the low energy constants ^3^4 (for a definition of the lat- 
ter see Eq. (B44) and [13]) are shown in Table XVIII. Indeed, having a fixed dynamical heavy 
quark mass in the SU(3) x SU(3) theory is equivalent to the SU(2) x SU(2) theory, up to terms 
of 0((m//mv)^). The remarkable agreement between the LECs obtained directly from the fit to 
the chiral behaviour using SU(2) ChPT and those obtained by converting the SU(3) LECs using 
Eqs. (B40)-(B43) (see Table XVIII) is evidence that these terms of 0{{mi/ms)^) are small. Al- 
though reassuring, we stress that this in itself is insufficient to demonstrate fully the validity of 
the SU(3) ChPT at NLO. For example, if the m^ term in the expansion for the decay constant 
were large, this would appear in our NLO fits as a shift in /o but the behaviour with m/ would be 
equivalent to that in SU(2) ChPT (up to corrections of 0((m//ms)^) of course). We suspect that 
this is the case and use SU(2) ChPT to obtain physical results nevertheless, for the remainder of 
this section, we take our SU(3) results at face value and compare them to previous determinations 
of the LECs. 

We first consider the ratio of the decay constants in the two and three flavor chiral limit. From 
our fits to SU(2) and SU(3) ChPT we obtain f/fo = 1.229(59) (here we quote only the statistical 
error), showing the influence of the strange quark loops. An important observation is that the 
value of /o, the pion decay constant in the SU(3) chiral limit, is much smaller than the measured 
pseudoscalar decay constant /// at, say, m/ = 0.01 and also much smaller than /;;:. This is due 
to the large NLO correction shown in Fig. 17. The small value of /o may be a contributor to 
poor convergence, since the chiral logarithms come in with a factor of I//0 • Of course, one 
can rearrange the series and expand in m^,//| consistently to NLO and possibly improve the 
convergence. We are exploring these options, but if such a rearrangement, which only effects the 
series at NNLO, markedly changes the convergence properties, one is still lead to the conclusion 
that the series is not well controlled. 

If we compare our result for /o to phenomenological estimates by Bijnens [78], our value of 
/o = 93.5(7.3) MeV (statistical error only) turns out to be substantially lower than the preferred 
value from Bijnens, which is 124.0 MeV (NNLO, alternative fits also published there give values of 
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1 14.7 (NLO), 99.6 (NNLO), and 1 13.7 (NNLO) MeV). His ratio f^/fo is also different from ours 
(in contrast to the SU(2) case, where for the ratio /;r// a better agreement has been achieved): we 
obtain a value of 1.33(07) (statistical error only), whereas Bijnens' preferred fit suggests fn/ fo ~ 
1.05. Interestingly, the Nf = 2 + \ lattice simulation by MILC [65, 66] also observed a higher 
ratio of f^/ h = l-21(5)(l3'^), which translates into a central value for /o = 106.0MeV (using 
their value for /;;: = 128.3(0.5) (^'^jg), we do not quote an error for this number, since we do not 
know the correlation between /;;: and fn/fo-) 

For the ratio B/Bq we obtain B/Bq = 1.03(05). The very small deviation from 1 is perhaps sur- 
prising, and indicates good agreement with the predictions of the large Nc approximation and only 
small Zweig rule violations in this case. A comparison of the LO-LEC Bq is not directly possi- 
ble, since this number depends on the renormalization scheme (see Sec. VI B 2). Instead, Bijnens 
quotes a number for the renormalization scheme independent and dimensionless ratio 25ow„rf/m| 
of 0.736 (from the preferred NNLO-fit; alternative fits give: 0.991 (NLO), 1.129 (NNLO), 0.958 
(NNLO)). From our SU(3) fit we obtain 0.995(41) (statistical error only), which agrees better with 
the alternative phenomenological fits than the preferred one. 

The ratio of E = f^B/2, the chiral condensate in the two-flavor theory, and Eq = /q Bq/I in the 
three-flavor theory, can also be compared directly. We obtain E/Eq = 1.55(21), whereas MILC 
quotes a value of 1.52(17)(^jg). It should be noted however, that we obtain slightly different 
values for f/fo and for B/Bq than the MILC collaboration. 

In Table XIX we compare the NLO-LECs to phenomenological NLO and NNLO fits [78] and 
the results of MILC's Nf = 2+ I dynamical lattice simulations [65, 66] (all LECs in that table 

(3) (3) 

are quoted at the scale A;^ = 770MeV). In the fits by Bijnens, L\ ' and Lg ^ were set to zero (at 
A;^ = 770 MeV ^ mp) as motivated by Zweig-rule and large A^^ arguments. Whereas for the latter 

(3) 

our result agrees with this assumption within the statistical uncertainty, for L^ we observe some 

(3) 

discrepancy (without taking systematic errors into account). Interestingly, our NLO results for Lg ^ 

f3) 

and Lg ' agree very well with Bijnens' numbers from the NNLO fit, but not for his NLO fit (for 
which no systematic error is given however). Within the reported uncertainties, our NLO LECs 
agree nicely with the set of values quoted by the MILC Collaboration. 

In Fig. 18 we summarize some of the results presented in this section. Here we plot our mea- 
sured values for the decay constant, converted to physical units, versus the degenerate valence 
pseudoscalar mass squared, along with the results of fits to SU(2) and SU(3) ChPT. The data for 
degenerate quarks is denoted by filled symbols and for non-degenerate quarks, open symbols are 
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used. The graph shows that we see small effects for non-degenerate quarks. The results of SU(2) 
and SU(3) partially quenched ChPT fits to our data are shown, and both fits agree well with the 
data. From the (now determined) fit functions, we plot the unitary SU(2) chiral extrapolation. Two 
of our data points lie on this curve, but one sees that our measured decay constant for a pseu- 
doscalar mass of 420 MeV is about 30% above the chiral limit value, /. We also plot the SU(3) 
chiral limit curve. For this curve, the horizontal axis is the unitary meson mass in the theory where 
all quarks are degenerate. Therefore, none of our measured values must lie on this line. Fig- 
ure 18 shows graphically the large difference between /o and our measurements. This is graphical 
evidence of the poor convergence of NLO SU(3) ChPT, with LECs as determined from our data. 
In summary, we have found that we can fit our data well using NLO SU(3) ChPT for average 
valence quark masses < 0.01, corresponding to pseudoscalar masses below 420 MeV. The SU(3) 
LECs we determine are generally in reasonable agreement with continuum phenomenology and 
other lattice results (note however, the small value of /o which we find). We see that the LECs 
from our SU(3) and SU(2) fits agree well when a conversion is performed from SU(3) to SU(2). 
However, since we find large corrections at NLO we would expect significant ones also at NNLO. 
We are therefore not at all confident that the systematic errors in the SU(3) LECs are currently 
under control. In addition, the fits do not agree with our data when we extend the NLO SU(3) 
ChPT theory to the kaon mass scale. It is possible that with more data it will become feasible to 
perform NNLO SU(3) fits and to control the systematic uncertainties with sufficient precision. We 
will investigate this in the future but we stress that whether or not this proves to be the case, we 
can happily use SU(2) PQChPT to obtain predictions for physical results (including quantities in 
kaon physics as explained in section II B ). This is what we do in this paper. 

VIII, Bk 

A. Pseudoscalar Bag Parameter on the Lattice 

The kaon bag parameter is defined as the ratio of the neutral kaon mixing matrix element and its 
expectation value from the vacuum saturation approach, 

Bk = ^ ' ^^ 2 ' (116) 

j,JK^K 
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where mk is the neutral kaon mass, /k the kaon decay constant, and 



9A5=2 



(^s{\-r5)r^d)(^s{i-r5)rnd) (iiv) 



is the local, effective four quark operator, which couples to the left-handed quarks and induces a 
change in strangeness by AS = 2. In our simulations, we define the corresponding pseudoscalar 
bag parameter for a meson made from either valence or sea quarks with masses m^ and niy, which 
we measure by fitting the ratio 






BM = g .J'"!' r/w ^ (118) 



to a constant Bj^y over some range in time t. Here, t denotes the time coordinate at which we insert 
the four quark operator, tsic and ?snk are the temporal coordinates at which we insert Coulomb gauge 
fixed (spatial) walls of anti-kaon and kaon interpolating operators, respectively. The correlators 
using wall source and local sinks, which we construct from the WL -propagators obtained with 
quark masses nix and niy, are defined as 

^ yev 
^ yev 

"^ApitJsnk) = TT £(^(j,0^^(3^,0^w(?snk)/''?w(?snk)), (121) 

^ yev 
where qw{t) and q{yj) denote a Coulomb gauge fixed (spatial) wall smeared quark field and 
a local quark field, respectively. A summation over all spatial points for the local operators is 
performed, and a balancing volume factor W = l} included giving a spatial volume average that 
is statistically efficient. In addition, we average over propagators obtained with periodic and anti- 
periodic boundary conditions in the time direction, which results in a doubled time extent available 
for the plateau fit. 

The values measured on the 24^ x 64, L^ = 16 lattices are given in Table XX. There, the fit to 
the plateau was performed over the range t E [12,52] using an uncorrelated fit, since, due to the 
large time extent, correlated fits became unstable. Previously, our collaboration also obtained 
the pseudoscalar bag parameter on smaller 16^ x 32, Ly = 16 lattices with dynamical light quark 
masses m/ G {0.01,0.02,0.03} and m., = 0.04 at the same gauge coupling [19, 34, 79, 80]. For 
comparison we quote those values in Table XXI. In this case the plateau was fitted over a range 
?G [12,22]. 
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We will have to extrapolate our measured values to the physical light and strange quark masses 
for which PQChPT in the SU(2) x SU(2) formulation will be used. These fits are discussed in 
the next part of this section; the result obtained for Bk was already published in [19], but here 
we provide some more details on the chiral fitting procedure. In the remainder of this section, we 
employ SU(3) x SU(3) PQChPT to extract the bag parameter in the chiral limit and briefly discuss 
the data from the 16^ lattices. 

B. SU(2) X SU(2) Chiral Fitting Procedure 

For the chiral extrapolation of Bk, we use SU(2) x SU(2) ChPT for the kaon sector, as introduced 
in Sec. II and applied to the kaon masses and decay constants in Sec. VIA. Specifically we use 
Eq. (B48) to extrapolate to the physical average light quark mass m„j determined in Sec. VIB 1. 
Here we consider the data from the 24-^ lattices at values of niy = 0.03 and 0.04 for the heavier 
quark mass approximating the strange quark mass. Of these two values only the latter describes a 
truly unquenched quark and we account for possible effects of a partially quenched strange quark 
treatment in our discussion of systematic errors. The fits using SU(2) x SU(2) ChPT for the kaon 
sector are shown in Fig. 19, whereas the fitted parameters are given in Table XXII. The parameters 
B and / have been fixed to their values obtained in the combined SU(2) x SU(2) ChPT fits for 
the pseudoscalar masses and decay constants, cf. Sec. VI. We applied a cut in the light valence 
quark mass of nix < 0.01 . The low energy constant bi (niy) multiplying the dynamical light quark 
mass in Eq. (B48) shows almost no dependence on the heavy quark mass, whereas b2{iny) and 

(K) 

Bpg{my) vary by 12% and 4%, respectively, when going from niy = 0.03 to 0.04. Also the value 
of the kaon bag parameter extrapolated to the physical average light quark mass, B^dy, increases 
by approximately 4% due to this change in niy. 

The interpolation to the physical value of the strange quark mass (obtained from the kaon mass, 
cf. Sec. VIB 1) is done linearly as shown in Fig. 20. Included in the plot (and also in Table XXII) 
is a point at niy = 0.02 for information only. It was not included in the linear fit and interpolation 
because this supposedly heavy quark mass is likely too light to be sufficiently separated from 
the simulated light quark masses to enable convergent effective treatment. However, the linear 
interpolation does not seem to deviate much even at this point, giving us confidence that non- 
linear effects in the strange quark mass region at least do not show up at the level of the numerical 
precision achieved. Finally, we quote for the unrenormalized physical value of the kaon bag 
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parameter B^^ = 0.565(10). We will now discuss the renormalization of i^ff^^ and estimate all 
significant contributions to the systematic error. 

1. Renormalization 

The renormalization of ^ff^^ has been treated in [21] using the RI-MOM non-perturbative renor- 
malization technique. Couplings to wrong chirality four quark operators are in principle admitted 
by the presence of a non-zero residual mass. It has been demonstrated in that paper that they are 
highly suppressed and small, using non-exceptional momenta to remove the spontaneous chiral 
symmetry breaking effects that might obscure this fact at the modest, accessible lattice momenta. 
Thus we only need to consider a simple, multiplicative renormalization. Here we quote the renor- 
malization factors for the regularization independent (RI) and the modified minimal subtraction 
(MS) schemes both at /i = 2GeV as well as the renormalization group invariant (RGI) result: 

zg(2GeV) = 0.910(05)stat(13)syst, 5|'(2GeV) = 0.514(10)stat(07)ren, (122) 

Z|f (2GeV) = 0.928(05)stat(23)syst, B^{2GeY) = 0.524(10)stat(13)re„, (123) 

Zff = 1.275(10)stat(25)syst, Bk = 0.720(13)stat(14)re„, (124) 

where the first error is the (combined) statistical error and the second error the systematic error 
from the renormalization. 

2. Systematic Errors 

a. Finite Volume Ejfects The spatial volume in our simulation is approximately (2.7 fm)^, 
therefore, from [81] it follows that the difference between the B^y measured on the finite lattice 
volume and the result in infinite volume is negligible for all points except the lightest ones (both 
valence quarks at a mass of 0.001), where the difference may be as large as 2%. Excluding the 
point nil = 0.01, mx = niy = 0.001 from the fit, the final result for B^^^ remains almost unchanged. 
Comparing points with m/ = 0.01 from the smaller volume simulation (16^ lattice, ~ (1.8fm)^ 
spatial volume) from Table XXI with the corresponding ones from the 24-^ simulation (Table XX), 
statistically marginal differences of up to 1% are observed, see Sec. VIII D for a more detailed 
discussion. Conservatively, we adopt this as an estimate for finite volume effects affecting our 
final number for Bk- (Note, that for the lightest point in the 16-^ ensemble (nix = niy = nii) we 
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have m^L ^ 4.0, whereas for our lightest valence pion in both 24^ ensembles {m^c = niy = 0.001, 
nil = 0.005 or 0.01) we still have nijiL ^ 3.4, so it is reasonable to assume the finite volume effects 
in those points to be comparable to the observed effect from comparing 16^ and 24^ lattices.) 

b. Scaling Effects Having only data at a single value for the lattice spacing does not allow to 
quantify the scaling effects affecting our data. The CP-PACS Collaboration calculated Bk in a 
quenched calculation using the Iwasaki gauge action with Domain Wall valence quarks at two 
different lattice scales, namely 1/a = 1.81(4) and 2.81 (6) GeV (determined from the p-meson 
mass) keeping the physical volume approximately fixed [82]. Extrapolating their observed scaling 
violation to our coarser lattice spacing, we would expect a 3.5% scaling effect. We assume a 4% 
systematic error, which is in agreement with the scaling violations discussed in conjunction with 
the meson masses and decay constants in Sec. VIC 4. 

c. Interpolation to the Physical Strange Quark Mass Currently only ensembles with a fixed 
value for the dynamical heavy quark mass are available. We estimate the effect of this 15% too 
high sea quark mass from the measurements done at the 16^ lattices. There dynamical light quark 
masses of m/ = 0.02 and 0.03 have been simulated, which are closer to the physical value of the 
strange quark mass than the dynamical light quark masses available from the 24-^ data. Comparing 
Bxy for nix = 0.01, niy = 0.04 (lightest valence quark mass and dynamical strange quark mass) at 
the two aforementioned light sea quark masses, an increase of 3% is observed (Table XXI). In that 
case, two dynamical quark flavors were changed by Am = 0.01, whereas the single strange quark 
mass only has to be changed by 0.0057. Accordingly, the systematic error has to be scaled down 
by 2 (from number of flavors) times ^1.8 (from scaling Am), resulting in an 1% effect. (This is 
an exclusive sea quark effect and should not be confused with the 4% difference discussed above, 
when changing the valence niy from 0.03 to 0.04 in the kaon itself.) 

d. Extrapolation in the Light Quark Mass Including terms up to NLO in our chiral fit functions, 
we have to assign a systematic error to our extrapolation resulting from neglecting NNLO and 
higher order terms. The linear fit to the 16^ data gives a 6% higher value for Bk, see Sec. VIII D. 
Assigning this difference to the (here included) NLO terms, the size of NNLO contribution can 
be estimated to be 2% by scaling the observed 6% difference by fhi/fhs ~ 0.4 taken at the lightest 
value for quark mass in the linear fit, m; = 0.01, and the strange quark mass at the physical point 
m.. 
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3. Final Result 



Combining the 1% finite volume, 4% scaling, 1% heavy quark mass interpolation, and 2% ChPT 
extrapolation systematic errors with the one from the NPR (except for 5^^ of course), our final 
result in the different renormalization schemes considered reads: 

Bf = 0.565(10)stat(27)syst, (125) 

5|I(2GeV) = 0.514(10)stat(07)ren(24),yst = 0.514(10)stat(25)comb, (126) 

5f^(2GeV) = 0.524(10)stat(13)ren(25)syst = 0.524(10)stat(28)comb, (127) 

Bk = 0.720(13)stat(14)ren(34)syst = 0.720(13)stat(37)eomb , (128) 

where "stat" denotes the statistical error, "syst" the systematic error as discussed above and "ren" 
the error due to renormalization, "comb" is the combined systematic error from the latter two. 

C. Fits to SU(3) X SU(3) PQChPT 

Here we apply Eq. (B22) from SU(3) x SU(3) PQChPT to fit our data. That means, we now also 
try to describe the dependence on the heavier (valence) quark mass by (PQ)ChPT. As discussed 
in Sec. II B 5, if Eq. (B22) were applied for a fixed value of the heavier valence quark mass (also 
fixing the heavy sea quark mass), it would naturally revert to the kaon SU(2) x SU(2) form. If we 
were able to describe the valence mass dependence with the SU(3) x SU(3) PQChPT form up to 
the strange quark mass, this form would be applicable to a determination of 5^^ However, it is 
also possible that this form can be used to describe only light valence masses m^, niy to obtain the 
SU(3) X SU(3) low energy constant 5pg representing the pseudoscalar bag parameter in the three 
flavor theory in the limit of all three masses being zero, which is of phenomenological interest 
[83]. The only remnant terms involving a large mass in such a fit of the SU(3) x SU(3) PQChPT 
arise from the fixed dynamical heavy quark mass m/,, and we will estimate a systematic error from 
this. 

We observe the same behaviour as for the meson masses and decay constants: a reliable fit includ- 
ing quark masses up to the strange quark mass is not possible with only terms of NLO included 
in the fit formula. In Fig. 21 fits to the 24^ data are shown with two different ranges for the mass 
cut mavg. While the fit with mavg < 0.01 describes the data inside the fit range but badly fails at 
the heavier points, going to mavg < 0.02 already fails to describe the data at the lowest masses. 



61 

The fitted parameters are given in Table XXlll. Here we fixed Bq to its value obtained from the 
SU(3) X SU(3) fit for the meson masses and decay constants (mavg < 0.01), cf. Sec. VII. Due to 
the failure to describe the data with one quark mass as heavy as the strange quark, a determination 
of 5)^' is not meaningful using this ansatz. However, limiting the fit range to small masses, we will 
estimate the SU(3) x SU(3) LEC 5pg. In Fig. 22 the dependence of 5pg on the applied mass cut is 
shown. From that we conclude that at least within the statistical error the result is stable, therefore 
we quote (from a fit using mavg < 0.01): 

Bff = 0.266(26)stat, (129) 

5g'(2GeV) = 0.242(24)stat(03),en, (130) 

5^P(2GeV) = 0.247 (24)stat(06),en, (131) 

5^s = 0.339(33)stat(07),en, (132) 

using the same renormalization factors as in Eqs. (122-124) (first error statistical, second error 
systematics from renormalization). One has to keep in mind that this result was obtained by 
extrapolating m/j from its value used in the simulation to zero (i.e., extrapolating m/, -^ —m^es)- 
Considering Eq. (B22) which was used to fit the data, one sees that — except for the analytic terms 
in /disc where it enters via Xri — the heavy quark mass only enters via x, which multiplies the LEC 
d. By only using one dynamical heavy quark mass value, this LEC is exclusively determined by the 
dependence of Bxy on the light dynamical quark mass, a very mild dependence as can be seen by 
comparing the columns for m/ = 0.005 and 0.01 in Table XX. Two problems might arise. Firstly, 
the slope with respect to variations of m/, may differ from the slope with respect to variations of m/ 
due to higher order corrections. Secondly, with our fit procedure we might not be able to reliably 
determine this parameter. We assign the full size of this NLO estimate of the contribution. 






PS 'X(A-rrf^\2 



H^nfoY 






(133) 



to be the systematic error of our chiral extrapolation in the heavy quark mass, which is a 10% effect 
for the chosen cut-off, mcut = 0.01. (Going to higher values in mcut the contribution increases to as 
much as 28% for mcut = 0.02.) The NLO-contribution of the extrapolation in the light quark mass, 
which should be well under control by our fit procedure, we get from the difference of linearly 
extrapolating the measured B^y values at the dynamical points with nii = nix = my = 0.01 and 
0.005 to nil = nix = niy = — mres ('•£. fhi = mx = thy = 0) and the chiral extrapolation to this point 
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(m/, = 0.04 is fixed in this procedure). By multiplying the resulting difference of 

A = Bf^^^imi = 0,mh = 0.04) - 5j^^(m, = 0, m,, = 0.04) ^ 0.154 (134) 

by Xl/i'^^fo)^^ we determine the uncertainty due to neglected NNLO contribution to be 8%. 
For the systematic errors due to finite volume and scaling, the same applies as for the SU(2) x 
SU(2) analysis of Bk, so we assume them to be 1% and 4%, respectively. Eventually, with the 
10% strange and 8% NNLO light quark mass extrapolation, the 1% finite volume, 4% scaling and 
NPR systematic error, we obtain (first error statistical, second combined systematic error): 

5g' = 0.266(26)stat(36)eomb, (135) 

Bff{2GcY) = 0.242(24)stat(03)ren(33)syst = 0.242(24)stat(33)eomb, (136) 

5^^'(2GeV) = 0.247(24)stat(06)ren(33)syst = 0.247(24)stat(34)eomb, (137) 

5^s = 0.339(33)stat(07)re„(46)syst = 0.339(33)stat(47)comb. (138) 

Our final number for 5pg is in agreement with the phenomenological estimates from (i) the large 
Nc approximation, B^^ = 0.38 ±0.15 [83] (B^ in their notation), and (ii) B^^ = 0.39 ±0.10 [84] 
using the QCD-hadronic duality approach (which is close to the chiral limit value, see remark in 
[83]). 

D. Results from Simulations on a Smaller Volume 

As mentioned in the beginning of this section, our collaboration first performed a study of the 
pseudoscalar bag parameter at a smaller volume of (1.8fm)^ with higher dynamical masses 
nil £ {0.01,0.02,0.03}, see Table XXI. It turned out, that these masses were not light enough 
to probe the regime of ChPT, as was revealed by comparing a simple linear fit in the light quark 
mass with a fit to (PQ)ChPT. The extrapolations from these fits gave similar results within the 
numerical uncertainty, suggesting that the chiral logarithms occurring in NLO were not repre- 
sented correctly, or — in other words — the influence of the linear terms in the PQ;^PT formula 
was overestimated. Here, we just quote the result of extrapolating the three unitary points 
(nix = nil E {0.01,0.02,0.03}, niy = nih = 0.04) by a linear fit to the light quark mass at the 
physical point m„j 

^idV'"'"'('«/. = 0-04) = 0.611(08) , (139) 
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which gives the value of the kaon bag parameter at a slightly too high strange quark mass. This has 
to be compared with B^^y = 0.5789(97) at niy = 0.04 from Table XXII, yielding a 6% difference 
between these two approaches. (More details on the fits to the 16^ data have already been published 
in [80].) 

Also this data enables us to check for possible finite volume effects. The spatial volume increases 
by a factor 1.5^ = 3.375 when going to the 24^ lattices. Fig. 23 shows a comparison of the 
measured Bxy values at m/ = 0.01 for the two different volumes, revealing only small changes 
in the range of valence quark masses 0.01 < m^^my < 0.04, where data is available from both 
simulations, see also Table XXIV. The maximal deviation of 0.8% is observed for the lightest 
points (m^, niy) = (0.01 , 0.02) and (0.02, 0.02) and decreases down to 0.4% for the heaviest point 
(0.04,0.04). (Accidently, (0.01,0.01) shows almost no, /.e., < 0.1% deviation in the central value.) 

IX. VECTOR MESON COUPLINGS 

In this section we discuss the couplings of the light vector mesons V to vector and tensor currents. 
These couplings fy and fy are defined through the matrix elements: 

{0\q2{0)fqm\V{p;X)) = fynivE^ (140) 

(0|^2(0)a^V(0)|V(;^;A)) = if^{^) (e^p'-elp^) , (141) 

where p and A are the momentum and polarization state of the vector meson V{p;X) and £x is 
the corresponding polarization vector. The tensor bilinear operator qiO^'^qi (and hence fy{iJ,)) 
depends on the renormalization scheme and scale ji. The final results will be quoted in the MS 
scheme at jl — 2 GeV We have presented preliminary results for the ratios fy /fy in Ref. [85]. 

A. Experimental Determination of fy 

The decay constants fy can be determined experimentally. For the charged p and K* mesons, one 
can use T decays to deduce fp and f^* as illustrated by the diagram in Fig. 24, where the curly line 
represents the VF -boson. From the measured branching ratios one obtains the following values for 
the decay constants [86]: 

Br(T~ -^ p-Vr) = (25.0±0.3)% => fp- ~ 208MeV (142) 

Br(T"^i^*-v^) = ( 1.29 ±0.03)% => /^*-~217MeV. (143) 
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One can also determine /pO from the width of the decay of the p^ into e+e~ which gives fpO ~ 
216(5) MeV. Similarly from the width of the decay ^ -^ e+e^ one deduces f(j, ~ 233 MeV . 
The couplings fy are not known directly from experiment but are used as inputs in sum-rule cal- 
culations (see, for example, refs. [87, 88]) and other phenomenological applications to 5-decays 
(see, for example, refs. [89-92]). Previous lattice results for the vector meson couplings are re- 
called below; determinations obtained using QCD sum-rules are nicely reviewed in [87]. We now 
present our calculation and results for fy / fy, which can then be combined with the experimental 
values of fy to obtain fy . For the ^ we neglect the Zweig suppressed disconnected contribution. 



B. Lattice Calculation of fy /fy 

In order to determine fy/fy it is sufficient to calculate the following zero-momentum correlation 
functions for large values of the Euclidean time t: 

aj^it) ^ Y.{^\V'^^{t,x)V'^\0)\0) = 3f'y'f'y^mye-"^^^/^coMmy{T/l-t)) (144) 

x,i 

Crfit) ^ Y.^0\Tll{t,x)V;\0)\0) = 3f^''f^^mye-'^^^l^smh{my{T/l-t)), (145) 

x,i 

where V,- and T^i represent the vector and tensor currents and / = 1 , 2, 3 is a spatial index, si and S2 
label the smearing at the sink and at the source, respectively. From the ratio 

^(0 = ^ = T- tanh(mK(r/2 -0) (146) 

Cyy JV 

we readily obtain the ratio of (bare) couplings. 

C. Results 

In Table XXV we summarize our results for the vector meson masses from fits to (144) on the 
DEG data set (cf. Table II) and from fits to (144) and (145) on the UNI data set. In the latter case 
we average over various choices of the source smearing function (cf . Table IV) while always using 
a point sink. We restrict our study to the unitary case in which the sea and valence quarks have 
the same masses. The bare strange quark mass is always fixed at 0.04. On the UNI data set, again 
averaging over the same choices for the source and the sink, we also evaluate the ratios fy/fy. In 
each case Eq. (146) exhibits well pronounced plateaus which we fit to a constant. 
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In Table XXVI we present the bare values of fy / fy. It can be seen that the measured results are 

obtained with excellent precision. We have also compared our results with those obtained on a 

16-^ X 32 X 16 lattice for m; =0.01, 0.02 and 0.03 [85] (the properties of the ensembles on the 16^ 

lattice have been presented in ref. [34]). No significant finite volume effects were found. 

^From Fig. 25 it can be seen that the dependence of the bare fy/fv on the masses of the light 

quarks is very mild and so we restrict our chiral extrapolation to linear and quadratic functions in 

the quark mass as shown in the figure. For the ratio of bare couplings in the chiral limit we obtain: 

f /-r fl 

^=0.619(15)(18); ^ = 0.6498(62) (60); ^ =0.6838(32)(22) , (147) 

fp JK* T(p 

where the central value corresponds to the linear extrapolation and the second error is the differ- 
ence between the results from the linear and quadratic extrapolations. 

The bare results in Eq. (147) were obtained with the notional strange quark mass of m/, = 0.04 
rather than the physical value of m^ = 0.0343 (see Table XI). The values of the ratios in Eq. (147) 
are very similar for the p , K* and (^ mesons and we correct for the change in m^ by linear inter- 
polation in the valence quark mass (m/j is fixed at 0.04). Thus, for example, for the K* meson we 
interpolate between the K* and the p: 

^(m, = 0.0343) = ^(m, = 0.04) + — — ^ ^ (0.0343-0.04) , (148) 

JK* JK* (0.04 + mresJ 

where A = f^*/fK* {fn^ = 0.04) — fp /fp. After carrying out a similar extrapolation for fl / f^ the 
corrected bare values are then 

f /-r fl 

^=0.619(15)(18); ^ = 0.6457(62) (60); ^ =0.6753(32)(22) . (149) 

fp fK* fcj) 

We determine the renormalization constants non-perturbatively using the Rome-Southampton 
method and run the results to 2 GeV. The results for the renormalization constants of the ten- 
sor and (local) vector currents were presented individually in ref. [21] and for the ratio we find 
Zj/Za = Zj/Zy = 1.1101(92) ~ 1.11(1), where Za is the renormalization constant of the local 
axial current. The relation between the ratios of bare and renormalized matrix elements is then: 
/f (2GeV) _ Zr(2GeVa) /f'^^^(a) _ ^ ^w^ /v'"^(^) .^^q. 

fv Zy /bare " ^ ' /bare 

In the MS scheme with /i = 2 GeV we finally obtain: 

fT fT fl 

^=0.687(27); ^=0.717(12); ^ = 0.750(8). (151) 

Jp JK* J(p 
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These results can be compared with previous quenched lattice results which we summarize in Table 

'P 



XXVII. The QCDSF/UKQCD collaboration has also presented the resuh /J = 168(3) MeV using 



an Nf = 2 0{a) improved clover action with a range of lattice spacings (0.07 < a < 0. 1 1 fm) [93]. 
Combining our result for the ratio from Eq. (151) together with the experimental value for fp we 

P 



obtain a smaller value fl = 143(6) MeV . 



X. CONCLUSIONS 

This paper reports the results of simulations of 2+1 flavor QCD with domain wall fermions at a 
single lattice spacing, with a larger volume than has been achieved previously. The ensembles 
we have generated have a length of ^ 4, 000 molecular dynamics time units, and we have seen 
a good distribution of global topological charge, indicative of our algorithm's ability to sample 
phase space. With our (2.74 fm)^ spatial volume, and improvements in the RHMC algorithm, our 
lightest pion made of dynamical quarks has niji = 331 MeV and our lightest valence pion has a 
mass of niji = 242 MeV. For these two cases, we have nijiL = 4.60 and 3.36, respectively. 
Because of the good chiral properties of domain wall fermions at non-zero lattice spacing, the only 
correction to continuum chiral perturbation theory at NLO order is the inclusion of the residual 
mass, m,,„ in the quark mass appearing in the ChPT formula. This means we have the same number 
of ChPT LECs to fit as for continuum physics, although the LECs we determine can differ from the 
continuum values by 0{a^) terms. The small number of free parameters in our fits is a powerful 
advantage of working with a fermion formulation with good chiral and flavor symmetries. To have 
good control over our chiral extrapolations, we have developed SU(2) ChPT for kaon physics, 
which does not assume the kaon mass is small. We have fit our data to both SU(2) and SU(3) 
ChPT, at NLO, and find both approaches give good fits to our data. However, for SU(3) ChPT, the 
NLO corrections are very large - as much as 50% of the size of the LO term for the pseudoscalar 
decay constants in the range of quark masses where we have data. This poor convergence of the 
series does not seem to be due to a single term, such as having the dynamical heavy quark mass 
close to the physical strange quark mass. Our fits have used pseudoscalar particles with mass 
< 400 MeV, where SU(3) ChPT might be expected to be more accurate than we have found it to 
be. We are investigating this further, making use of the full continuum NNLO expressions of [76]. 
For our most reliable extrapolations to the physical values for the light quark masses, we turn to 
SU(2) ChPT for kaons. Here we have found good fits, with the NLO corrections no larger than 
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35% of the LO terms for the light quark mass values we used in our simulations. This is still a 
somewhat large NLO contribution, but makes the dropping of NNLO terms more sensible. We 
have tried to estimate our systematic errors by including possible analytic NNLO terms in the fits. 
For our central results, we use m^, m^z and mk to set the scale and determine m„j and m^. For m^j, 
there are no light quark chiral logarithms and for m^d and m^ we use the SU(2) ChPT formula to 
perform the extrapolation. We can then predict the decay constants, f-n, fx and their ratio. The fits 
also give values for the LECs of SU(2) ChPT. 

We have estimated our systematics through several methods. For the systematic errors from our 
chiral fits, we have tried varying the fit ranges and also including analytic NNLO terms, to see the 
effects of these changes on our results. Finite volume effects are also estimated through ChPT, 
as are the errors due to the heavy quark in our simulations being 15% heavier than the physical 
value. For finite lattice spacing effects, we expect them to be 0((3Aqcd)^ or about 4%. Scaling 
errors no larger than this have been seen in the quenched case [94], and, since dynamical quark 
loop effects change the scaling at a higher order in as, similar scaling is expected for the full QCD 
case studied here. The size of this scaling error also agrees well with the deviation of our central 
value of fn from the experimental results, so we have adopted a uniform 4% 0{a^) error for all 
of our results. (Simulations are underway at a smaller lattice spacing, so that we will soon have 
better control over the a ^ extrapolation.) 

YovBk, which has already been published [19], we have given more details about our use of SU(2) 
ChPT to extrapolate to the light quark limit. We have also measured the vector meson couplings 
for light vector mesons. 
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Our major results come from using SU(2) ChPT fits and are summarized below: 

/ = 114.8(4.1),tat(8.1)systMeV, (152) 

5^(2GeV) = 2.52(0.11)stat(0.23)ren(0.12)systGeV, (153) 

E^(2GeV) = (255(8)stat(8)re„(13)systMev) , (154) 

h = 3.13(0.33)stat(0.24)syst, (155) 

k = 4.43(0.14)stat(0.77)syst, (156) 

A3 = 666(1 10)stat(80)systMeV, (157) 

A4 = l,274(92)stat(490)systMeV, (158) 

m^(2GeV) = 3.72(0. 16)stat(0.33)ren(0.18)systMeV, (159) 

mf^(2GeV) = 107.3(4.4)stat(9.7)ren(4.9)systMeV, (160) 

mud : m, = 1 : 28.8(0.4)stat(1.6)syst, (161) 

fn = 124.1(3.6)stat(6.9)systMeV, (162) 

fK = 149.6(3.6)stat(6.3)systMeV, (163) 

fK/h = 1.205(0.018)stat(0.062)syst, (164) 

5f^(2GeV) = 0.524(0.010)stat(0.013)ren(0.025)syst. (165) 

These results have total errors for the decay constants of about 6% and for the MS quark masses 
of about 10%. The systematic error due to 0{a^) errors will be reduced by the lattices we are 
generating at a smaller lattice spacing (1/a ^ 2.3 GeV) and lighter quark masses on a similar 
physical volume. The errors on the quark masses are largely due to renormalization and should 
improve as better renormalization conditions are chosen [21]. 

With the expected increases in computer power over the next few years, we should be able to 
push to much lighter quark masses and minimize our reliance on chiral perturbation theory in the 
extrapolation to physical results. This is possible for domain wall fermions, since the chiral limit 
is decoupled from the continuum limit. With lighter quark masses, we should be able to check the 
convergence of the chiral perturbation theory expansion and achieve better control over the low 
energy constants. Finally, the ensembles that we are generating in these basic studies of the low 
energy properties of QCD are useful for a wide variety of other measurements, including heavy 
quark systems and hadronic weak interaction matrix elements. 
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APPENDIX A: NOTATION, CONVENTIONS 

1. Quark Masses 

Here we discuss our notation for quark masses. In Section II B, where we discuss Kaon ChPT, 
continuum notation is used and m represents the total quark mass. Throughout the rest of the 
paper, we need to distinguish between the input bare quark mass in the domain wall fermion 
formulation and the total quark mass, which includes the additive contribution to the quark mass 
from finite L^. 

Input (sometimes also called bare) quark masses for domain wall fermions will be denoted by the 
symbol nix, where the index X E {/, h, x, y, ud, s} distinguishes between different types of quarks: 

• m/: dynamical light input quark mass (two degenerate flavors), 

• nih'. dynamical heavy input quark mass (one flavor), 

• mx, niy-. valence input quark masses, 

• niud'- average light input quark mass at the physical point, i.e., a meson made of two quarks 
of mass niud acquires the experimentally measured mass of the neutral pion (this is related 
to the input quark masses of the up- and down-quarks via: niud = {mu + md)/2), 

• nis'. heavy input quark mass at the physical point (strange quark mass), i.e., a meson made 
of one quark of mass m„j and one of mass m^ acquires the quadratically averaged experi- 
mentally measured mass of the neutral kaons. 

The total quark mass, which has the residual mass mres added, will be denoted by in, 

mx = nix + Wres- (Al) 

The total renormalized quark masses at the physical point carry a superscript indicating the renor- 
malization scheme used. For example, physical masses in the MS scheme renormalized at a scale 
ji =2GeVread 

m^(2GeV) = Z^(2GeV) ■ (m^^ + mres) = Z^(2GeV)m„d, (A2) 

mf^(2GeV) = Z^(2GeV) ■ (m, + mres) = Z^(2GeV)m,, (A3) 

where Z^^ (2 GeV) is the mass renormalization constant. 
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2. Pseudoscalar Quantities 

The pseudoscalar quantities considered in this work (masses, decay constants, and bag parameters) 
will carry two indices, denoting the quark content. For example the mass of a pseudoscalar made 
from two valence quarks with mass nix and niy will be labeled nixy, the decay constant of a meson 
with one valence quark and one dynamical heavy quark is f^h and the pseudoscalar bag parameter 
of a meson with one quark at the physical average light quark mass and one valence quark is B^dy 
Also we will use the obvious symbols for quantities at the physical point: 

mTt = niudud, inx = in^ds, 

In Judud 7 JK Juds •) 

Bk = Buds- 
3. Low energy constants in ChPT 



The following table lists the LECs we use in the PQChPT formulae: 





LO 


NLO 


PQ SU(3) X SU(3) 


Bq, fo, 

dX 
Bps 


.(3) .(3) .(3) .(3) 

b, c,d 


PQ SU(2) X SU(2) 


B,f 


(2) (2) (2) (2) 


PQ SU(2) X SU(2) kaon 


5WK),/WK), 


Ai(m/,), hiinh), h{inh), ^4{inh), 
bi{mh),b2{mh) 



The renormalization scale dependent NLO LECs are defined at a scale A;^ . For the decay constants, 
we use the normalization such that the experimentally measured value for the pion decay constant 
is fn ^ 132 MeV ^ v^-93MeV. 
Further, we will use the following abbreviations for (tree-level) masses in the formulae: 



Xx = 25omx, SU(3)xSU(3) 
Xx = 25 mx, SU(2)xSU(2), 



(A4) 

(A5) 



where it should be clear from the context, whether the SU(3) x SU(3) or SU(2) x SU(2) expres- 
sion is to be used. In addition, we define the average dynamical mass to be 



m = -{2mi+mh) 



(A6) 



and the dynamical r/-mass 



f^ri = -{mi + 2mh) 
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(A7) 



(Note: this is not the tree-level r^-mass at the physical point, which would be {niud + 2ms)/3.) 

APPENDIX B: PQCHPT FIT FUNCTIONS 

1. SU(3) X SU(3) 

The following formulae hold for Nf = 2 + I sea quark masses and two valence quarks (SU(5|3)). 
The formulae for masses and decay constants were derived from [46], while the ones for the bag 
parameter 5ps can be found in [98]. 

a. Squared Pseudoscalar Mass 
a. non-degenerate (nix 7^ ^yj 



m^y^m^y^ niy-. 

Xx~T- Xy 



m 



'xy 



+ 



1 



\Xx ~ Xl) \Xx ~ Xh) 1 Xx [Xy ~ Xl) [Xy ~ Xh) 

[Xx ~ Xy) [Xx ~ Xrj) Ay [Xy ~ Xx) [Xy ~ Xrj) ^v 



24;r2/2 

[Xri —Xl)[Xri ~- Xh) ■. Xr]_ 
[Xi] ~ Xx) [Xi] ~ Xy) Ay _ 



Xy\og^' 



(Bl) 



nix^niT^ y^ niy-. 

Xx' Xy 



m 



xy 



1 



/o JO 

[Xy — Xl)[Xy — Xh) , Xy 
{Xy-XxY ^'^^\ 



[Xx Xy) [Xx Xh)Xx~T~[Xx Xl)['^Xx Xh) [Xx Xl)[Xx Xh)Xx 



[Xx Xy , 



1 AjX 



[Xx Xl)[Xx Xh) 
Xx ~ Xy 



(B2) 



b. degenerate (m^ = niy) 
• nix = my ^ nin: 



1 + 1^(243, _^(3,,.^16p^(3,_^,3,)^_ 



, -^ f ^Xx Xl Xh _ [Xx Xl)\Xx Xh) \ , %x 

24n^fi[\ Xx-Xn {Xx-Xn? ^^' ^A| 



, iXx — XuiXx — Xh) f . ,1 Xx\ , \Xi] ~ Xi)\Xi] ~ Xh) ^^ 1 Xr\ 
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(B3) 



mx = My -^ m 



m^ 



/o /o 



24;r2/o' 



XlXh , ^Xx JXI ■-'Xh , /r, N, Xx 

2^ + 2 + {^Xx-Xi-Xh)^og^ 

A/ ■ 



(B4) 



c. Pion and Kaon masses m\ and m\ From Eq. (B3) we obtain for the degenerate meson mass 
in unquenched QCD {mx — mi) 



.2 -.. >. ,^(243,_^(3,).^16,„,3,_,,3,, 
JO JO 



4 = Z/ ■ 1 + 1^(24'^ -Lr)X + -^{2q^-Lf)xi 



Xl 1 
^, log— --^^ log 



•A| 2' 



247t^fi 
and from Eq. (B 1) for the kaon mass (m^ = mi,my = mh): 



Xr\_ 



(B5) 



^ih = 



Xi^Xh 



8 



JO JO 



1 + -^(2L ^^ -L^^OZ + T7(24^ -4 0(Z/ + Z/, 



+ 



■Zrilog 



Z77 



24;r2/o2^" ^Al f 



(B6) 



So that for mx = mi= mud and my = mh = ms we get m| and m|^, respectively. These expressions 

agree with [13]. 

For a "kaon" made from light sea quark and valence quark {m^ 7^ m/), one obtains: 



nij^^mn: 
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m 



/o Jo 



lAK^fl 



Xx Xh 



Xx Xr\ Xh 



Xt] 



Xx Xrj ^Y Xr\ Xx ^y _ 



(B7) 



nix^m^. 



m 



u - 2^■^ + 5(24"-4'')^ + 4(24"-4")to+xo 

/o h 



Xx 



lAK^fl 



Xx~Xh + {1Xx-Xh)'^og^ 



(B8) 



Whereas for a "kaon" made from a valence quark {m^ 7^ m;,) and a heavy sea quark, we have: 



m^^mr^: 



m 



/o /o 



1 



Zx Z/ 



Xx 



- , .^o Zxlog -y H ^ Xn log -y 

24;r2/o LZx-Zri A2 Z77-ZX ' Ai 



(B9) 



nir^mr 



f^xh = 



Xx'\' Xh 
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.(3) ,(3)^ 



■(3) ,(3)^ 



24;r2/2 [ 



1 + -^{ILl'' -L\'')x + -iilL',' -L's')iXx + Xh) 
h 

{2Xx-Xi)^og^+Xx-Xi 



Xx 



(BIO) 



b. Pseudoscalar Decay Constant 



a. non-degenerate (nix 7^ tny) 



mx^nir^^ niy-. 
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/. 



24 



(3), 



■(3), 



xy 



Jo h 



Xx'T~ Xl -, Xx'T Xl I Xy "T Xl 1 Xy "T Xl 
^^ log — -^ + -^-^ log - 



8^Vo' 



2A| 



2A| 



Xx I X/i 1 „ „ Zx + Z/i I X>' + a/i 1 Xy ' Xh 



■log- 



2A| 



+ ■ 



8 



2A| 



96;r2/o' 



Xx 



iXl - Xx) (Xh - Xx) , ^, (Z/ - Xy) iXh - Xy)_ ^ . _ y )2^ ^°^ ^ 



Zrj Ax Zrj A,y 



"T/tT] \/Cy Ax 



\Xi] ~ Xx) \Xi] ~^ Xy) ^Xx — Xr\ Xy ~ Xr\ 
{Xl-Xx){Xh-Xx) {Xl-Xy){Xh-Xy. 



Xrj Xx 



Xt] Xy 



Xx Xy 



(BU) 



nij 



lUr^ 7^ niy-. 



24 



(3), 



■(3), 



/xy = /O ■ <; 1 + -:2^rX + ^2^ (Xx + Xy) 
JO JO 

Xx'T Xl -, Xx'T Xl Xy ~T Xl 1 Xy ~T Xl 
—, log^T7^ + — ^^ log- 



^n^fo 



2A| 



2A| 



, Xx + Xh 1 Xx + Xh , Xy "T Xh 1 Xy "T Xh 
+ ^ log ^,9 +^ log 



2A| 



8 ° 2A2 ' 8 

1 \Xl + xl^Xi + Za - 5z.y) + Zx (5Z>-(Z/ + Xh) - ^XiXh - 2Z; 



96;r2/o' 



Zx (3Zy + Z/Zfe - 2Zj(Z/ + Xh)) + Z>- (2Z/Za - XyiXi + Z/.; 

IZx ~ Xy) 



log — 

Xy ^ 



XyXlXh 



(B12) 



Zj. degenerate (nix = my) 



24. 



(3), 



8 



■(3), 



/xx = /O ■ <; 1 + -12^4'^ + 1^^5 Xx 
JO JO 



l6K^f^ [ 



U, ,, \ 1 Xx'T Xl Xx'T Xh 1 Xx 'T Xh 
x + Z/)log ^,^ + — ^ — log 



2A| 



2A2 



(B13) 
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c. Pion and Kaon decay constants fjt and /k From Eq. (B13) we obtain for the pion decay 
constant in unquenched QCD (m^ = mi) 



1 



24-(3) 8 (3) 

/o 
Xi + Xh, Xi+Xh 



/// = /o ■ <; 1 + j2Lrx + y2Lrxi 

Jo /o 



16;r2/o' 



■log- 



Xi 



+ 2Z/log 



, -^i->.^2 



and from Eq. (Bl 1) for the kaon decay constant (m^ = nii, my = m/,): 
//;. = /o ■ <! 1 + ^4'^^ + ^Lf\xi + Xh) 



64;r2/o' 



o Z/ + Za 1 Xl^Xh , ^, , Xl , ^, ■, Xrj 



(B14) 



(B15) 



Again, for mx = mi= niud and niy =mh = nis we get /^ and /a:, respectively, and these expressions 

agree with [13]. 

For a "kaon" made from a light sea quark and valence quark {nix ¥" ^i), one obtains: 



nix^m^. 



fix = 



/o-<; i + ^Lfz + 44'^(z/+Zx) 



+ 



1 



8^Vo' 

Z/ + Z/J 



i*'"^! 



Xx~T Xl ■, Xx' Xi 

■log- 



4 



2A| 



8 



1 Z/ ~r Xh Xx + ZA 1 Xx'T Xh 
log-r-^ r log- 



2A| 



2A| 



(Zx Z/J \Xr\ Xh 



Xr\ 



12Ur,-Zx)2U,-Z/)^'^''^A| 



, Xx\Xh~ Xl~ Xr\) ~Xr\Xl\Xh~ ^Xx) , Zx 

mxrj-xxr ^^A2 

IZ/ AxJ iXfe Xx) 

12Ui7-Zx) 



(B16) 



m^—^m^. 
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fix 



Jo JO 

•J ., 1 Xi Xx~T~ Xi -, Xx' Xi 

Xl~^ Xh , Xl ' Xh Xx'T Xh ■, Xxi Xh 

■log-— ^ — log- 



8 



2A| 



8 



2A| 



Xi iXh - ^Xx) + XxiXh + Xx) Xi (Xh -Xi) . x^ 
24Xx +12U/-Zx) °^A| 



(B17) 



Whereas for a "kaon" made from a valence quark (m^^ ^ m^) and a heavy sea quark, we have: 



m^^^nin 



my-^^nir 



L 



xh 



fo-i l + ^L^^h + ^Lf\xx + Xh) 
Jo Jo 



1 



8;r2/o' 

Xx I Xh 



Xx'T Xi ■, Xx' Xi I XhiXii Xh I Xi 
—, — log-— ^ + — - — log- 



2A| 



2A| 



1 Xx ~T~ Xh , Xh 1 _ Xh 



1 



96;r2/2 



Xx-Xh-Xx- — jlog- 
Zr/ — Xx ^ Xh 



+Zr] \Xh — Xx 



Xr\ Xl 



Xv 



Xv 



'osf logf 



Xri — Xx^ Xx Xf] Xh Xi] 

\Xl ~ Xx) [Xh ~ Xx 



Xr\ Xx 



(B18) 



/. 



xh 



1 



24.(3) 

'4 
Xx + Xl 



fo-i l + ^L)^'x + -:24\Xx + Xh] 
Jo Jo 



8^Vo' L 

Xx ~r Xh 



, Xx + Xi ^Xh + Xi, Xh + Xi 
log ^77^ + ^ — log- 



2A| 



2A| 



1 Xx + Xh , Xh 1 Xh 

8 ^°^^Ar + T^°^A|, 



1 



96;r2/2 



Xx + Zx (Z/ - 4Zfe) + ^XxXh - Xlxi 
■^XxyXx ^ Xh) 



Xx{Xh-Xi)+Xh{Xi-Xh) ^^ X. 

IZx^Z/jj Xh 



(B19) 
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c. Scale-dependence of the T^LO-LECs 



The NLO-LECs L^. (A;^) at a renormalization scale A;^ can be converted to a different scale Ay 
using [13] 



■(3)/A' 



(3) 



A. 



with 



Lr(A;,)=Lr(A,) + ^log-L 



1 ^ 3 ^ 11 ^ 5 

r4-g, r5--, r^-— , Tg--. 



(B20) 



(B21) 



d. Pseudoscalar Bag Parameter 



Here we follow the notation of [98] for the LEC's relevant for the pseudoscalar bag parameter and 
denote them by b, c, and d. Also here we only work in (PQ) SU(3) chiral perturbation theory. 
a. non-degenerate (nix 7^ tny) 



1 



5x>. -5ps<;i+ 2 



^conn + ^disc , Zx + Z>', , IZx Xy) , 3 
Xx'T' Xy ^ Xx~T~ Xy ^ . 



(B22) 



3 ,/ 

/conn = 2(^-^ + ^.y) (^-1-log 



Xx + Xy 



-2i^Xx + Xy)Xx'^Og^ - -{Xx + 3Xy)Xy^Og^, 



-/disc = L + Iy+I; 



y-r^ri, 



(B23) 

(B24) 



{^Xx + Xy){Xl-Xx){Xh-Xx) 



1-log 



Xx \ 

7q) 



'^[Xf] ~ Xx) 
OXx + Xy){Xi-Xx){Xh-Xx) , {^Xx + Xy){Xi-Xx){Xh-Xx) 



+ ■ 



'^[Xf] ~ Xx) [Xy ~ Xx) \Xr\ — Xx) 

2(Z/ - Xx) iXh - Xx) - OXx + Xy) (Xh - Xx) - (3Zx + Xy) iXi - Xx) 



^{Xii — Xx) 



*. 1 A'-t 

■^-^log-j, 



ly = Ix{x^y), 



\Xx Xy) \Xx~\~ Xy~^^Xr\)\Xl Xr])\Xh Xrj) , Xrj 



^{Xx Xt]) [Xy Xrj) 



Xri^Og 



^¥ 



(B25) 
(B26) 
(B27) 
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For nij^^ m-q, nix^ niy, one obtains from the disconnected diagrams 



^disc Xl 



Xl + Xl OXh - 16z.v) + Z7;Z>-( 16Zfe - 9Z>0 + XhXy 

^XniXn-Xy) 
VzJ - 9XhXy + ^Xy + XnXyi-'^^Xh + Vz>-) + X^ iXh + 18z>0 



4(Zr]-Zj; 



1 Xri 
■log-y 



24 - 8^J]Z>- - Z>^ [xhXy + Xl {Xy - ^Xh) 

+XnXy ixi OXh - 9Xy) + 3;t:y(Zy - ^Xh) 

+Xq ixi {Xh - 2Xy) + Zy ( 1 5Xy -^Xh)) 



1 



2{XrT Xy)^^""^ 4 



Xy[X^i2Xh-^Xy)+XTlXyi9Xh " H^y) +Z> to + 2Z>- 

+Z/ ( - Zj (Zft - 2Z>0 + XyiXy - ^Xh) + Z77Zj(9z>- - 7Za)) 
fe degenerate (nix = fny) In this case /disc = and therefore 

1 



Bxx — ^PS^ 1 + 



PS 



24;r2/o' 






(B28) 



(B29) 



c. Kaon bag parameter Bk From Eq. (B22) one obtains the full (unquenched) QCD result 

(nix = nii,my = nih) 



Bih — Bnc < 1 



1 



PS 



16;r2/o2 

5Xi + ^Xh 

2{Xi + Xhr''^^^Al 



U/+Zft)(l+log- 



Xi + Xh \ 3xi + Xh , Z^ 
2A| y'^2U/ + Z;.)^' ^Ai 



^ 1 Xv ,Xi + Xh 2{xi-Xhf ,- 



6 ^3U/ + Z/^) 

so that at the physical point (m/ = m„^, m/, = m^) we have 



(B30) 



Bk 



^PS 



(AKhY 



Mi(l+log^) + 



4M2 ^A2 



+ 



IMl-M: 



M2 



AMI "^Al 6^^^ 3 



2c (Mi ~ Mi 



■2\2 



M2 



--(M^ + 2Mi 



(B31) 
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Here we expressed everything in terms of meson masses: M^ = Xud, M^ = j{Xud + Xs), Ml ~ 
Xn\mi=mud,mh=ms- (This agrees for the non-analytic terms with Eq. (12) of [81], analytic terms 
were not considered in that reference.) 

2. SU(2) X SU(2) 

The following formulae hold for Nf = 2 (degenerate) sea quark masses (mass m;) and two valence 
quarks. (SU(4|2)). The formulae for masses and decay constants were derived from [46], also cf. 
[99]. 

a. Squared Pseudoscalar Mass 
a. non-degenerate 



m 



i - ^ ■ ] 1 + §(2Lf -Lf )Z, + ^i2l^-l^)iX.^Xy) 



\6n^f 



^2v 4 .^. f2 

^^Zx log ^ + ^3^Z>' log ^ 

Xx Xy J^x '^y ^^ X 



(B32) 



b. degenerate 



m^ 



P(2Lf)-Lfto + 15(24^) _if; 



+ 



l6Tt^f 



Xx-Xi + {2Xx-Xi)^og^ 

^x 



(B33) 



c. Pion mass m^ From Eq. (B33) one obtains the degenerate meson mass in full (unquenched) 
QCD (m^ = mi): 



ml = Xi-h+^^{{2Lf^-Lf^)+2{2L, 



(2) ,(2). 



Xl + 



1 



16;r2/- 



-;i:ziog 



Xi 



(B34) 



which gives for mx = mi = mud the pion mass m| and agrees with [12]. Furthermore, we can relate 
the LECs 4^] g g from PQChPT to /^ in ChPT: 



4((242)-Lf) + 2(2Lf-L 



(2) ,(2). 
■4 ' 



I!. 



(B35) 



b. Pseudoscalar Decay Constant 



a. non-degenerate 



U = f-y + ^Lf^Xi + ^Lf\xx + Xy) 



1 



32;r2/2 
1 



iX. + Xi) log ^1^ + {Xy + Xi) log ^^ 



6AK^f 



Xx + Xy-^Xl + 



■^XxXy XiyXx'T' Xy) 1 Xx 
^ log — 



Xy Xx 



Xy . 



b. degenerate 



J XX — J \ ^ ^ r2^4 XI ^ r2^5 Xx 16;j;2/-2 & 



2A| 



c. Pion decay constant f^z 
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(B36) 



(B37) 



8 ^Or(2) : ^(2)^ 



/„ = /.{i + ^(2Lr+Lr)*-3,,^,.^.^. 



^' log^}. 



(B38) 



Again, for mx = mi= niud this reproduces the pion decay constant f-^ and agrees with [12]. Fur- 



(2) 



thermore, we can relate the LECs L\ ^ from PQChPT to l\ in ChPT: 



4(2Lf+Lf 






(B39) 



c. Conversion from SU(3) x SU(3) to SU(2) x SU(2) 

The conversion from the SU(3) LECs to the SU(2) ones can be done (inch terms up to NLO) 
following [13] (see [26] for NNLO, note the latter organizes the expansion in a slightly different 
way compared to the former). 



/ = /o<^ 1 



^Og^.9 + ^7^-4 Xs 



32n^fi "2A2-/2 4 



B = BnU 



Xs 



log 



72;r2/2 ^3A2 



||.|(2.f-.f')-} 



r. = AiiL 



■(3) ,(3) 



+ 8(24^)-Lf 



57^0+^°^^) 



II = SLf^+4L^^^ 



64^0+^°^^) 



(B40) 

(B41) 
(B42) 

(B43) 
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Where one has to use Bq to evaluate the Xx on the right hand side. The relation of the SU(2) LECs 
in PQChPT, if'^ to the unquenched SU(2) LECs /^ 4 is given by Eqs. (B35, B39). To define them 
at the scale of the pion mass m^i, one has to use 

i^^/f-log§ (B44) 

73 = -2 '74 = 2. (B45) 



3. SU(2) X SU(2) for the kaon sector 

Here, we choose to denote by nix always a light and by niy a heavier valence quark mass. 
The additional low-energy constants appearing here (LO: B^ \mh), /(*■-' (m/,), B^^inih), NLO: 
^\,2,3,'x{fny), b\^2{tn}i)) are in general dependent on the dynamical heavy quark mass m/, and the 
valence heavy quark mass niy. To simplify the notation the argument m;, has to be viewed as 
placeholder for both of those. For more details see Sec. II. 



a. Light-Heavy Squared Pseudoscalar Mass 



ml = bC^) (m.) «, { 1 + '^x, + ^Z, } (B46) 



b. Light-Heavy Pseudoscalar Decay Constant 



„, _/mK)^i + M^;f, + M^ 



Xx'T' Xl 1 Xx 'T' Xi , Xi ■^Xx 1 Xx 

log^T7^ + -. log 



2 ^ 2A2 4 ^A2 



(B47) 



c. Light-Heavy Pseudoscalar Bag Parameter 



D dW/ \) ^ , ^^^^h) ^, , himh) Xl 1 Xx , ,-r,.g, 

Bxy = B'ps'imh)^ 1 + -j^Xi + -J^^--y2^^^7J ( ^ ^ 
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d. Conversion from SU(3) x SU(3) to kaon SU(2) x SU(2) 

For the LECs appearing in the kaon mass and decay constant formulae (Eqs. (B46, B47)), we 
provide the relation to the SU(3)-LECs as well: 



iK) 



(3), 



^(3), 



/^^K,m,) = Ul + -L\'>Xh + -piL\'^Xy + 



f^ 



/, 







32;r2/o2 



Xy Xh 
.^Xy-2Xh 



Xy 



,y , XhXy \qo.^^ 



Xh 1 Xh ,, 1 Xy , 

■ylog2Xj-»l"g2XJ+(3z,-2z.)^'-A| 



X 



XhXy 



rlog 



^Xh Xh "T Xy 1 Xh "T Xy 



-log 



2A| 



{^Xy-^XhY ^3A2 : 

5WK,m,) = 5o|l + §16(2Lf)-Lf) + |8(243)-Lf) 



(B49) 



^ ^Zylog-y + 2: 



Z, ^^g2z. 



11%^ fl [ 3z, - 2z/, •'■ ^ A| 3z, - 2;^;. ^ 3A2 



(B50) 



AiK,m,) = 32(2Lf -Lf ) + :^ 



^KXy-Xh? Xy 
{^Xy-2Xh? ^A| 



+ 2 

Xy 



„ 1 X/i I Xh 



(3Zy-2Z/.)2^'"^3A2 ^ 3^,-2^;. 



(B51) 



A2K,m,) = ^ + 8(2Lf ) -Lf ) + 16^(Lf - (24^) -L^ 

Xy Xy ^ 



+ 



1 



4;r2 



Xy Xh , Xy Xh , Xh . Xh\jXy Xh) , -^Xh 

2{3Xy-2Xh) °^a|"4^°^2A| 9^,(3^,-2^;.) °^3A| 



(B52) 



(3) -^ 

X'i{mh,my) = \6Ly + —^ 



1 Zy 4^ft-ll^/^^y + 8Zy 
2^°S2A|- 



0Xy-2Xh)^ 

Xh) 
{3Xy-2Xhy 



4(3z, - IXhf ^' °^ 3A2 + i3Xy - 2Xh? ^^ A| j 



(B53) 



f3) 1 

X4{mh,my) = 4Lg + 



16;r2 



Xv 



-log 



2Xh 



4°^2A| 3z, - 2Xh °^ a| " 2(3z, - 2Xh) "^ 3A| 



(B54) 



APPENDIX C: FINITE VOLUME CORRECTION FOR PSEUDOSCALAR MASSES AND 

DECAY CONSTANTS 



In the following we will provide the finite volume corrections for the meson decay constants and 
squared meson masses obtained in PQChPT for Nf = 2+1 and Nf = 2 sea quarks (cf. [56-61]). 
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The corrections are labeled A^/ and A^f , respectively, for decay constants and squared masses 
of mesons made of quarks with masses m^ and niy in a finite (spatial) volume L?. Labeling decay 
constants in finite volume f^ and those in infinite volume f^°° we have the following relations: 

4"°° = 4(1 -A^/) (CI) 

= /o(i + zPT4), (C2) 

4 =/o(i+ZPT4)(i+A^/) (C3) 

= /o(i + ZPT4 + A^/) (C4) 

and similar for squared meson masses. Here we summarized the NLO contribution in infinite 
volume (PQ)ChPT (as given in the previous section) by X^^xy (These equalities hold up to terms 
of NLO, higher order terms are neglected. Note further, that in case of SU(2) x SU(2) the /o has 
to be replaced by /. 

The Bessel functions of imaginary argument (modified Bessel functions of the 2"*^ kind) Kn{x) 
enter the expressions for the finite volume corrections via 

5,w = i£^iM, (C5) 

53{x) = 2Y,Ko{rx), (C6) 

55W = -Y.-Ki{rx), (C7) 

where the argument is typically x = ^/XxL. For the numerical implementation, we made use of 
the multiplicities m{n) and rewrite the sum as 

£/(r) = £m(n)/(v/^), (C8) 

where sum is evaluated until the relative change is less than the required precision e or a maximum 
number A^ for n is reached. (Typically, we use £ = 5 ■ 10^^ and A^^ = 100, but checked that going 
to £ = 5 ■ 10^^, A^ = 1000 does not change the result.) The multiplicities for n <20 can, e.g., be 



found in [62]. Here we list m{n) for n = 1, . . . , 100: 
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6,12,8,6,24,24,0,12,30,24,24,8,24,48,0,6,48,36,24,24, 
48, 24, 0, 24, 30, 72, 32, 0, 72, 48, 0, 12, 48, 48, 48, 30, 24, 72, 0, 24, 
96,48,24,24,72,48,0,8,54,84,48,24,72,96,0,48,48,24,72,0, 
72,96,0,6,96,96,24,48,96,48,0,36,48,120,56,24,96,48,0,24, 
102,48,72,48,48, 120,0,24, 144, 120,48,0,48,96,0,24,48, 108,72,30. 



(C9) 



Note that we have the following relations involving the df. 

d 



dX 



xSiiVxL) = -HVxL). 



dZ 



&i{^L) = lH,{^L) 



Furthermore, by doing the substitutions in the finite volume correction part A^ 



'xy 



5i(v^L) ^ logfl, 
HVX^L) -. -(l+log^) 



L^SsiVX^L) 



1 

Xx 



(CIO) 
(Cll) 



(C12) 
(C13) 

(C14) 



one obtains at NLO the non-analytic part of the corresponding X^^xy (i-^- without the analytic 
terms multiplying the LECs). 

1. SU(3)xSU(3) 
a. FV Correction for the Squared Pseudoscalar Mass 
a. non-degenerate (m^ 7^ niy) 



m^^^nin ^ niy 



xy 



\Xr\ — Xl)\Xr\ — Xh) 

24K^fi[{Xv-Xx){Xv-Xy] 

[Xx ~ Xl) [Xx ~ Xh 



Xn^iiVx^L) 



+ 



[Xx Xy)[Xx Xrj, 



\Xy XxJKXy Xi]) 



(C15) 
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m^^nir^ (niy j^ nin) 



^xy 



24;r2/n' 



\Xx Xl)\Xx Xh) 



Xx X' 



HVXxL) 



'J 



, [■^Xx Xl Xh)[Xx Xy) [Xx Xl)\Xx Xh) _ o / nr:-T\ 

H -, ■■ — XxOii^xL) 



[Xx Xy, 



iXy-Xl){Xy-Xh) 
[Xx Xy) 



XyMVX^L) 



(C16) 



b. degenerate 
• nixT^mn 



vLm 



Uii^fl 



1/77 ILx) 



, [Xx Xr\)\^Xx Xh~Xl)~[Xx Xl)[Xx Xh) ^ c / nrrT\ 

H : XxOi [^/XxL) 



[Xx Xr\ , 



[Xx Xl)[Xx Xh) 



Xx Xrj 



SiiVTxL) 



nix^mn 



vLm 



XxSi iVx^L) - {2Xx -Xl- Xh)HVX^L) 



24K^fi [ 

[Xx — Xl) [Xx — Xh) j2^ ( /t;-t\ 
n. ^ ^5[y Xxi^) 



(C17) 



(C18) 



b. FV Correction for the Pseudoscalar Decay Constant 



a. non- degenerate m^ 7^ niy 

• nix^nin ^niy 






1 



Xx + Xi 



+ 



Xx + Xh 



5i 



Xx + Xi 



L + 



Xy'Xlc' I Xy'Xlr 



8 
1 



81 



96;r2/o' 



Xx + Xh T \ I X>' "T Xh c 
[Xl Xx) [Xh Xx 



Xy ' Xh 



Xx Xrj 



HVXxL) - (^^ 5>')(^^^ ^'h ,{^L) 



Xy Xr\ 



I I [Xl Xri ) [Xh Xrj ) [Xl Xl] ) [Xh Xrj ) _ ^ [Xl Xi] ) [Xh Xrj / \ y s ( /y— /■ ^ 



+ 



[Xx — Xf] ) [Xy — Xf] ) [Xx — Xl] ) [Xy — Xf] ) 

Xh ~T~ Xl — ^Xx [Xl — Xx) [Xh — Xx) ^ [Xl — Xx) [Xh — Xx 



Xrj Xx 



[Xrj Xx) 



[Xrj Xx) [Xy Xx) 



Xx5ii^/x^L) 



^Xh~\~Xl ^Xy [Xl Xy)[Xh Xy) r. [Xl Xy)[Xh Xy) \ ^ ;; / /v~T\ 



Xt] Xy 



[Xri Xy, 



[Xrj Xy)[Xx Xy 



(C19) 
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m^—^m^,my^ m^ 



'-^xy 



STt^f^ 



Xx + Xl. ( IXx + Xlj\ , Xy + Xl. ( 



Xy ' Xl 



, Xx'^Xhc- I Xx~\~Xhj\ , Xy^Xhc- I Xy^ Xh J 

-\{Xi-Xx){Xh-Xx)L^5s{^xL) - ^^'~^J^^^y^'' 5,{^yL) 



Xy Xx 



^(^_^ 2Xx-Xi-Xh ^^ {Xi-Xx){Xh-Xx) ^^^^^^^^^-^ 



Xx Xy [Xx Xy) 

-[2Xx-Xi-Xh-2^^'~^''^^^''^'' 
+ 



Xx Xy 

Xh + Xl-2Xy {Xl-Xy){Xh-Xy) ^ iXl - Xy) iXh - Xy] 



Xx Xy 



[Xx Xy) 



HVTxL 

2 



[Xx Xy 



XyMVX^'L) 



b. degenerate, nix = t^y 



(C20) 



A^/ 



8^%' 



Xx + Xi 



5i 



Xx~T~ Xl J \ . Xx~T~ Xh o I Xx^Xhj 



(C21) 



2. SU(2)xSU(2) 



a. FV Correction for the Squared Pseudoscalar Mass 



a. non-degenerate (nix 7^ ^yj 



vLm 



1 



b. degenerate (nix = %j 



Xx Xl 






vLm^ 



1 



16;r2/2 



Xx^i ( v^L) - {xx-Xi)H VTxL) 



(C22) 



(C23) 
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b. FV Correction for the Pseudoscalar Decay Constant 



a. non-degenerate (m^ 7^ niy) 
1 



xy 



32;r2/2 
1 



[xx + Xi)Si 



Xx + Xi 



l] + [Xy + Xl)Sl 



+ 



6An^f 



- {xx-Xi)HVx'xL) - {xy-Xi)^{VTyL) 



Xy Xx Xx Xy 

b. degenerate (nix = tny) 



Xy + Xl 



(C24) 



(C25) 
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TABLE I: Improvements to the RHMC algorithm were implemented as the ensemble generation progressed. 

The table details the algorithms used for each ensemble and the range, in molecular dynamics (MD) time, 

where measurements were made. 

mi RHMC MD range Measurement range 



0.005 II 0-4435 900-4460 

0-549 

0.01 I 550-1454 

II 1455-5020 1460-5020 

0-1519 

0.02 I 1520-1884 

II 1885-3610 1800-3560 

0-492 

0.03 I 493-929 

II 930-3060 1260-3040 
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TABLE II: Calculation parameters for the three separate datasets FPQ, DEG and UNI. The range and mea- 
surement separation A are specified in molecular dynamics time units. A'meas is the number of measurements 
for each source position 4^. The total number of measurements is therefore A'^meas x A'^src. where Ngrc is the 
number of different values for t^rr- 



mi 


Dataset 


Range 


A 


"meas 


'src locations 




FPQ 


900-4460 


40 


90 


5,59 


0.005 


DEG 


900-4460 


40 


90 


0,32 




UNI 


900-4480 


20 


180 


0,32,16 




FPQ 


1460-5020 


40 


90 


5,59 


0.01 


DEG 


1460-5020 


40 


90 


0,32 




UNI 


800-3940 


10 


315 


0,32 


0.02 


DEG 

UNI 


1800-3560 
1800-3580 


40 
20 


45 
90 


0,32 
0,32 


0.03 


DEG 

UNI 


1260-3020 
1260-3040 


40 
20 


45 
90 


0,32 
0,32 
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TABLE IIL Combinations of quark source and sink smearings used for correlation functions for each of 
the datasets FPQ, DEG and UNI. A is the time component of the axial current density operator and P is 
the pseudoscalar density operator. XY-XY denotes the contraction of two quark propagators with X-type 
smearing at source and a Y-type smearing at sink. W=wall source, B=box source, H=hydrogen S-wave and 
L=local. 



Operator Type 


Dataset 


{^snk^src) 


FPQ 


DEG 


UNI 


(A,A) 


WL-WL 


BL-BL 


HL-HL 
LL-LL 


{PA) 






HL-HL 
LL-LL 


{A,P) 


WL-WL 
WW-WW 


BL-BL 
BB-BB 


HL-HL 
LL-LL 


{P,P) 


WL-WL 
WW-WW 


BL-BL 


HL-HL 

LL-LL 



TABLE IV: The different source and sink contraction of the quark propagators for the UNI dataset. V = Yi 
(/ = 1 , 2, 3) is one of the spatial components of the vector current and T = 041 are the calculated components 
of the tensor current. XY-XY denotes contraction of two quark propagators with a X-type smearing at 
source and a Y-type smearing at sink. H=hydrogen S-wave, G=Gaussian wave function (cf. [16]), and 
L=point. 



Contraction 


UNI Dataset 


i^mk^src) 


m=0.005, 0.02, 0.03 


m=0.01 


{v,v) 


HL-HL 
LL-HL 
HL-LL 


GL-GL 
LL-GL 
GL-LL 


{T,V) 


HL-HL 
LL-HL 
HL-LL 


GL-GL 
LL-GL 
GL-LL 
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TABLE V: Lattice values for the pseudoscalar meson mass and ;^^/d.o.f. calculated from an uncorrelated 
simultaneous fit to five correlators of the FPQ dataset on m/ = 0.005 and m/ = 0.01 ensembles. 







mi = 0.005 


nil = 


3.01 


nix 


niy 


mxy, 


X^ld.oi. 


mxy, 


X^IA.oi. 


0.001 


0.001 


0.1402(9) 


0.36(7) 


0.1432(10) 


0.18(6) 


0.001 


0.005 


0.1680(8) 


0.38(8) 


0.1706(9) 


0.21(6) 


0.005 


0.005 


0.1915(8) 


0.39(9) 


0.1938(8) 


0.25(6) 


0.001 


0.010 


0.1971(8) 


0.38(8) 


0.1994(9) 


0.26(7) 


0.005 


0.010 


0.2172(8) 


0.38(9) 


0.2194(8) 


0.25(7) 


0.010 


0.010 


0.2400(7) 


0.36(9) 


0.2421(7) 


0.24(7) 


0.001 


0.020 


0.2449(9) 


0.36(9) 


0.2473(9) 


0.30(8) 


0.005 


0.020 


0.2611(8) 


0.34(9) 


0.2635(8) 


0.24(8) 


0.010 


0.020 


0.2802(7) 


0.31(8) 


0.2825(7) 


0.22(7) 


0.020 


0.020 


0.3155(7) 


0.26(8) 


0.3178(6) 


0.23(8) 


0.001 


0.030 


0.2850(10) 


0.33(10) 


0.2877(10) 


0.33(10) 


0.005 


0.030 


0.2990(8) 


0.32(9) 


0.3016(8) 


0.26(9) 


0.010 


0.030 


0.3158(7) 


0.28(8) 


0.3183(7) 


0.24(8) 


0.020 


0.030 


0.3477(6) 


0.23(7) 


0.3500(6) 


0.26(10) 


0.030 


0.030 


0.3775(6) 


0.21(8) 


0.3797(6) 


0.29(11) 


0.001 


0.040 


0.3204(11) 


0.34(13) 


0.3234(11) 


0.36(12) 


0.005 


0.040 


0.3329(8) 


0.31(10) 


0.3358(8) 


0.31(11) 


0.010 


0.040 


0.3482(7) 


0.26(8) 


0.3508(7) 


0.29(11) 


0.020 


0.040 


0.3776(6) 


0.22(8) 


0.3800(6) 


0.30(12) 


0.030 


0.040 


0.4055(6) 


0.21(8) 


0.4077(5) 


0.32(13) 


0.040 


0.040 


0.4320(6) 


0.21(9) 


0.4341(5) 


0.34(13) 
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TABLE VI: Lattice values for the pseudoscalar decay constant calculated using method described in 
Eq. (92). 



m„ 


niy 


/^.(m,= 0.005) /^.(m,=0.01) 


0.001 


0.001 


0.0829(7) 


0.0860(7) 


0.001 


0.005 


0.0846(7) 


0.0878(6) 


0.005 


0.005 


0.0860(6) 


0.0891(6) 


0.001 


0.010 


0.0867(7) 


0.0900(6) 


0.005 


0.010 


0.0880(6) 


0.0910(6) 


0.010 


0.010 


0.0898(6) 


0.0928(6) 


0.001 


0.020 


0.0904(7) 


0.0928(7) 


0.005 


0.020 


0.0915(6) 


0.0946(6) 


0.010 


0.020 


0.0934(6) 


0.0962(6) 


0.020 


0.020 


0.0969(6) 


0.0995(5) 


0.001 


0.030 


0.0936(8) 


0.0972(8) 


0.005 


0.030 


0.0947(6) 


0.0979(6) 


0.010 


0.030 


0.0965(6) 


0.0994(6) 


0.020 


0.030 


0.1001(6) 


0.1025(5) 


0.030 


0.030 


0.1033(6) 


0.1056(5) 


0.001 


0.040 


0.0965(8) 


0.1003(8) 


0.005 


0.040 


0.0975(7) 


0.1008(7) 


0.010 


0.040 


0.0993(6) 


0.1022(6) 


0.020 


0.040 


0.1029(6) 


0.1053(5) 


0.030 


0.040 


0.1062(6) 


0.1083(5) 


0.040 


0.040 


0.1091(7) 


0.1111(5) 
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TABLE VII: Lattice values for the pseudoscalar meson mass and ;^^/d.o.f. calculated from an uncorrelated 
simultaneous fit to four correlators of the DEG dataset on m/ = 0.02 and m/ = 0.03 ensembles. 





mi = 


0.02 


mi = 


0.03 


nix 


ntxx 


;^2/d.o.f. 


ntxx 


X^/d.o.i. 


0.001 


0.1493(11) 


0.7(4) 


- 


- 


0.005 


0.1993(9) 


0.5(3) 


0.2016(11) 


0.4(1) 


0.010 


0.2471(8) 


0.4(2) 


0.2500(10) 


0.3(1) 


0.020 


0.3227(7) 


0.3(2) 


0.3261(9) 


0.3(1) 


0.030 


0.3849(6) 


0.4(2) 


0.3884(9) 


0.2(1) 


0.040 


0.4396(6) 


0.4(2) 


0.4429(8) 


0.2(1) 



TABLE VIII: Lattice values for the pseudoscalar decay constant calculated using method described in 
Eq. (93). 



Mx 


fxx{mi= 0.02) 


fxximi = 0.03) 


0.001 


0.0926(18) 


- 


0.005 


0.0947(15) 


0.0961(14) 


0.010 


0.0969(13) 


0.0985(13) 


0.020 


0.1021(11) 


0.1038(12) 


0.030 


0.1077(10) 


0.1092(11) 


0.040 


0.1133(9) 


0.1145(10) 
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TABLE IX: Decuplet baryon masses for degenerate valence quarks, "traj." refers to the Monte Carlo tra- 
jectories used for measurements. Numbers in ()'s give the total number of measurements, separated by 40 
trajectories in each case. To reduce effects of auto-correlations, the measurements were block-averaged into 
bins of size 80 trajectories before fitting. 



mi nij traj. (#meas) ;|^^/d.o.f. mass 

0.005 0.03 900-4460 (90), 920-4480 (90) 0.5(7) 0.9647(65) 
0.005 0.04 0.6(8) 1.016(51) 

0.01 0.03 1460-5020 (90), 1480-4040 (65) 1.1(1.1) 0.9809(88) 
0.01 0.04 1.2(1.2) 1.033(60) 

0.02 0.03 1600-3600 (50), 1900-3580 (43) 1.7(1.6) 1.030(11) 
0.02 0.04 2.0(1.8) 1.074(10) 

0.03 0.03 1020-3060 (52), 1320-3040 (44) 1.4(1.7) 1.040(10) 
0.03 0.04 1.8(2.2) 1.082(6) 
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TABLE X: Low energy constants obtained from SU(2) x SU(2) fits witli a valence mass cut mavg < 0.01. 
For convenience the LECs L- are quoted in units of 10^^ at two different chiral scales, A^ = 770 MeV and 
1 GeV. The definitions of the low energy constants are given in the Appendix. (Only statistical errors are 
quoted.) 



B f 


/3 U 


2.414(61) 0.0665(21) 


3.13(.33) 4.43(.14) 




(2Lf-Lf)(2LP)-Lf)) 


770 MeV 3.3(1.3) 9.30(.73) 
IGeV 1.3(1.3) 5.16(.73) 


0.32(.62) 0.50(.43) 
-0.71(.62) 4.64(.43) 



TABLE XI: Lattice scale and unrenormalized quark masses in lattice units. Note rhx = nix +'«res- Only the 
statistical errors are given here. 

a^^ [GeV] a [fm] niud fnud nig fhs rhud '■ fng 

1.729(28) 0.1141(18) -0.001847(58) 0.001300(58) 0.0343(16) 0.0375(16) 1:28.8(4) 
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TABLE XII: Central value and statistical error (stat) for the fit parameters from SU(2) ChPT and systematic 
error estimates from four sources. These are: finite-volume effects (FV), lattice artefacts (a^), the chiral 
extrapolation in the light quark masses (ChPT) and the unphysical value of the dynamical strange quark 
mass (rris). For details we refer to Sec. VIC. The LECs are quoted here only at a scale A), = IGeV, but 
the absolute value of the error is independent of A^. The total systematic error is obtained by adding the 
separate errors in quadrature, neglecting any possible correlations between them. 











systematic errors 






value 


Stat 


FV 


a^ 


ChPT 


m. 


total 


B 


2.414 


(.061) 


(.049) (.097) (.034) (.017) 


(.115) 


f 


0.0665 


(21) 


(13) 


(27) 


(36) 


(01) 


(47) 


(2Lf)-Lf).104 


-0.7 


(0.6) 


(0.1) 


(-) 


(1.8) 


(-) 


(1.8) 


(2L?)-Lf).104 


4.64 


(.43) 


(.95) 


(-) 


(3.20) 


(-) 


(3.34) 


Lf . 10^ 


1.3 


(1.3) 


(0.3) 


(-) 


(1.8) 


(-) 


(1.8) 


Lf -104 


5.16 


(.73) 


(.16) 


(-) 


(8.20) 


(-) 


(8.20) 


h 


3.13 


(.33) 


(.20) 


(.08) 


(.10) 


(.02) 


(.24) 


'h 


4.43 


(.14) 


(.05) 


(.08) 


(.76) 


(.04) 


(.77) 


mud 


0.001300 


(58) 


(23) 


(52) 


(22) 


(09) 


(62) 


fhs 


0.0375 


(16) 


(00) 


(15) 


(04) 


(08) 


(17) 


mud ■ m. 


1:28.8 


(0.4) 


(0.5) 


(1.2) 


(0.6) 


(0.6) 


(1.6) 


AiMeV] 


124.1 


(3.6) 


(1.9) 


(5.0) 


(4.4) 


(0.2) 


(6.9) 


MMcY] 


149.6 


(3.6) 


(0.4) 


(6.0) 


(1.0) 


(1.5) 


(6.3) 


h/fn 


1.205 


(18) 


(14) 


(48) 


(34) 


(12) 


(62) 



TABLE XIII: Finite volume correction factors obtained from the SU(2) and SU(3) PQChPT fits compared 
to the prediction by Colangelo, Diirr, and Haefeli (CDH) [62] at the unitary pion masses of 331 MeV 

{nix = my = mi = 0.005) and 419 MeV (m^ = niy = mi = 0.01). 



m;i:[MeV] 




331 
419 



0.09(.01) 0.14(.04) 0.13(.04) 
0.03(.00) 0.04(.01) 0.04(.01) 



-Rf[%] 
SU(2) SU(3) CDH 



0.36(.03) 0.56(.16) 0.32(.00) 
O.lO(.Ol) 0.17(.05) 0.09(.00) 
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TABLE XIV: The measured pion mass and decay constant at nix = niy = nii = 0.01, m/, = 0.04 at two 
different volumes (for 16^ data see [34]). We compare directly the measured ratio to the predictions from 
SU(2) and SU(3) ChPT and the value interpolated from Colangelo, Durr, and HaefeU (CDH) [62]. For 
explanation see Eq. (109) and text. 



fn 



ir 



measured 
16^ 



ratio 



0.24211(75) 0.247(3) 0.980(12) 
0.0928(6) 0.0895(7) 1.037(11) 



SU(2) 



1-R(240 

SU(3) 



CDH 



0.9961(04) 0.9944(15) 0.9942(22) 
1.0155(14) 1.0252(70) 1.0118(03) 



TABLE XV: Comparison of the SU(2) NLO-LEC h and k (defined at a scale of 139 MeV, cf. Eq. (B44)) 
with results from other Nf = 2 + 1 and Nf = 2 lattice simulations [66-68] and a phenomenological estimate 
[64]. 







Nf 


type 


h h 


this work, direct SU(2) fit 
this work, conv. from SU(3) 


2+1 
2+1 


DWF 
DWF 


3.13(.33)(.24) 4.43(.14)(.77) 
2.87(.28)(-) 4.10(.05)(-) 


MILC 
MILC'' 




2+1 
2+1 


stagg 
stagg 


2.85(.07)(-) 
0.6(1.2) 3.9(0.5) 


ETMC [68] 




2 


TM-Wilson 


3.44(.08)(.35) 4.61(.04)(.ll) 


CERN [67] 




2 


impr. Wilson 


3.0(0.5)(0.1) 4.1(0.1)(-) 
3.3(0.8)(-)'^ 


phenom.[64] 


2.9(2.4) 4.4(0.2) 



"direct SU(2) fit, from [66] 
''converted from SU(3), see [75] 

'^First number obtained without additional NNLO-term, second number from fit including NNLO-term, 
see [67] for details. 
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TABLE XVI: Comparison of light physical quark masses and their ratio obtained from recent lattice 
simulations using various types of fermion discretizations and different renormalization techniques (non- 
perturbative and perturbative). (In the cases where no error is given for the quark mass ratio, no error was 
quoted in the cited work. Since we do not know the correlations, we could not provide the error estimate 
here.) 

Nf type m^(2GeV)[MeV] mp(2GeV)[MeV] m, : m„rf 

using non-perturbative renormalization 

this work 2+1 DWF 3.72(0. 16)stat(0.33),en(0.18)syst 107.3(4.4),tat(9.7),en(4.9)syst 28.8(0.4),tat(1.6)syst 

RBC[69] 2 DWF 4.25(0.23),tat(0.26)ren II9.5(5.6)stat(7.4),en 28.10(0.38)stat 

ETMC[70] 2 TM-Wilson 3.85(0.12)stat(0.40)syst 105(3)stat(9)syst 27.3(0.3)stat(1.2)syst 

QCDSF[71] 2 impr. Wilson 4.08(0.23).stat(0.19)syst(0.23).scaie lll(6)stat(4)sy.st(6).scaie 27.2(3.2) 

using perturbative renormalization 

MILC[66] 2+1 Stagg. 3.2(0)stat(0.1)ren(0.2)EM(0)cont 88(0)stat(3)ren(4)EM(0)cont 27.2(0.1)stat (0.3), 

PACS-CS [72] 2+1 impr. Wilson 2.3(1.1) 69.1(2.5) 30( 

JLQCD[73] 2+1 impr. Wilson 3.54(1 j^)^^,,! 91.1(1^4/),otai 25.7(? 



sy St 
■71 
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TABLE XVII: Fitted parameters from SU(3) x SU(3) fits with a valence mass cut mavg < 0.01. For conve- 
nience the LECs l\ are quoted in units of lO^'* at two different chiral scales, A^^ = 770 MeV and 1 GeV. 
The definitions of the low energy constants are given in the Appendix. (Only statistical errors are quoted 
here.) 



B() = 2.35 (.16) /o = 0.0541(40) 
A. Z.f Lf (2Lf-4^))(2Lr)-Lr)) 

770 MeV 1.39(.80) 8.72(.99) -0.01(.42) 2.43(.45) 
1 GeV -0.67(.80) 2.51(.99) -0.47(.42) 5.19(.45) 
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TABLE XVIII: SU(2) x SU(2) low energy constants obtained directly from the SU(2) x SU(2) fits com- 
pared to those converted from the SU(3) x SU(3) chiral fits with a mass cut of mavg < 0.01. 



B 



f 



h 



h 



direct SU(2) x SU(2) 2.414(61) 0.0665(21) 3.13(.33) 4.43(.14) 
converted SU(3) x SU(3) 2.457(78) 0.0661(18) 2.87(.28) 4.10(.05) 



TABLE XIX: Comparison of fitted SU(3) NLO-LECs hf^ at K^ = 770MeV in units of 10""* to phenomeno- 
logically obtained values by Bijnens (see Table 2 in [78]) and dynamical Nf = 2 + \ staggered lattice sim- 
ulation (MILC), [65] and [66] (there, LECs were quoted at K^ =m^, conversion was done according to 
Eq. (B20).) 





^4 


-^5 


4^) Lf) (2Lr)-Lf)(24^)-Lf)) 


this work" 


1.4(0.8)(-) 


8.7(1. 0)(-) 


0.7(0.6)(-) 5.6(0.4)(-) 2.4(0.4)(-) 0.0(0.4)(-) 


Bijnens, NLO 
Bijnens, NNLO 


= 
= 


14.6 
9.7(1.1) 


= 10.0 5.4 =0 
= 6.0(1.8) 2.3^ =0 


MILC, 2007 


i.3(3.o)(;]:;j) 


13.9(2.0)(+2-|J) 


2.4(2.0)(+2;0) 7.8(1.0)(1.0) 2.6(1.0)(1.0) ?,A{\.G){+J^^) 



"Vox reasons mentioned in Sec. VII B, we do not quote any systematic error for parameters obtained from 

the SU(3) fits. 
''This value was derived from the quoted single values for Lj ' and Lg ' ; since we do not know the 

correlation between those two, we cannot provide the error estimate. 
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TABLE XX: Pseudoscalar B-parameter B^ for pseudoscalar mesons with valence masses nix, niy for two 
different values of the light sea quark mass mi (mi, = 0.04) on 24'' x 64 lattices, fit range 12-52 (uncorrelated 
fit) 









Bxy 


ntx 


niy 


m, =0.005 m, =0.01 


0.001 0.001 


0.4691(83; 


0.4698(53) 


0.005 0.001 


0.4910(66; 


0.4953(41) 




0.005 


0.5079(46; 


0.5119(32) 


0.01 


0.001 


0.5136(57; 


0.5209(39) 




0.005 


0.5267(37; 


0.5312(28) 




0.01 


0.5421(29; 


0.5459(25) 


0.02 


0.001 


0.5494(52; 


0.5593(42) 




0.005 


0.5588(3i; 


0.5629(26) 




0.01 


0.5697(25; 


0.5727(23) 




0.02 


0.5911(23; 


0.5935(20) 


0.03 


0.001 


0.5777(55; 


0.5886(54) 




0.005 


0.5850(32; 


0.5884(28) 




0.01 


0.5933(26; 


0.5957(22) 




0.02 


0.6105(23; 


0.6126(18) 




0.03 


0.6267(22; 


0.6287(16) 


0.04 


0.001 


0.6012(64; 


0.6131(73) 




0.005 


0.6070(38; 


0.6105(33) 




0.01 


0.6137(30; 


0.6160(24) 




0.02 


0.6280(23; 


0.6299(18) 




0.03 


0.6418(2i; 


0.6435(15) 




0.04 


0.6550(19; 


0.6565(14) 
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TABLE XXI: Pseudoscalar B-parameter B^y for pseudoscalar mesons with valence masses m^-, my for three 
different values of the light sea quark mass mi (mi, = 0.04) on 16'' x 32 lattices, fit range 12-20 (uncorrelated 
fit) 



mx m 



0.01 0.01 

0.02 0.01 

0.02 

0.03 0.01 
0.02 
0.03 

0.04 0.01 
0.02 
0.03 
0.04 

0.05 0.01 
0.02 
0.03 
0.04 
0.05 



Bxy 

mi = 0.01 mi = 0.02 m/ = 0.03 



0.546(8 
0.577(6 
0.598(5 
0.600(6 
0.617(4 
0.633(4 

0.620(5 
0.633(4 
0.647(4 
0.659(3 
0.636(5 
0.648(4 
0.660(3 
0.671(3 
0.682(3 



) 0.539(8: 


0.527(7) 


) 0.569(6: 


0.556(6) 


) 0.589(5: 


0.580(5) 


) 0.594(6: 


0.579(5) 


) 0.609(5: 


0.600(4) 


) 0.626(4: 


0.618(3) 


) 0.616(6: 


0.599(5) 


) 0.627(4: 


0.618(4) 


) 0.641(4: 


0.634(3) 


) 0.655(3: 


0.648(3) 


) 0.634(6: 


0.616(4) 


) 0.643(4: 


0.634(3) 


) 0.655(3: 


0.648(3) 


) 0.667(3: 


0.661(3) 


) 0.679(3: 


0.673(3) 
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TABLE XXII: Parameters and ;^^/d.o.f. for fits of Bxy to Eq. (B48) as shown in Fig. 19. The parameters 
B = 2.414(61) and / = 0.0665(21) were fixed from combined SU(2) x SU(2) meson masses and decay 
constants fits, Sec. VI; A;^. = 1 GeV. 



Bg)(m,) b,{my)/f b2{my)/f x^/d.ol B, 



rriv 



udy 



0.02 0.513(10) 
0.03 0.544(10) 
0.04 0.568(11) 



-5.2(1.6) 8.18(.65) 0.27 0.5249(86) 
-5.4(1.5) 6.85(.65) 0.20 0.5550(87) 
-5.4(1.7) 6.05(.71) 0.14 0.5789(97) 
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TABLE XXIII: Fit parameters for fits of Bxy to Eq. (B22) as sliown in Fig. 21 . Values are quoted for different 
cuts nicut in the average valence mass value. The paraemters Bq = 2.35(.16) and /o = 0.0541(40) were fixed 
from combined SU(3) x SU(3) meson masses and decay constants fits, Sec. VII; A), = IGeV. 



mcut 



B 



PS 



;^Vd.o.f. Bk = Buds 



0.01 0.266(26) -4.5(1.6) -0.20(.18) 0.68(.60) 1.8 0.484(15) 

0.0125 0.249(24) -2.5(1.5) 0.87(.13) 0.96(.61) 2.7 0.556(12) 

0.015 0.232(22) 0.2(1.5) 0.43(.13) 1.25(.59) 5.4 0.546(11) 

0.0175 0.223(22) 1.0(1.5) 0.94(.12) 1.49(.61) 6.1 0.570(12) 

0.02 0.205(19) 4.2(1.4) 0.40(.10) 1.89(.58) 10.2 0.553(10) 
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TABLE XXIV: Finite volume comparison for the pseudoscalar bag parameter measured on 16^ (? 

(1.8fm)^) and 24-' (« (2.7fm)-') for dynamical quark masses mi = 0.01, m/j = 0.04. 



nix 


niy 


r(163) 


r(243) 


(163) (243) 
tixy — tixy 


/r(243) 
jBxy 


0.01 


0.01 


0.546(8; 


0.5459(25) 


0.0(1.5) % 


0.02 


0.01 


0.577(6; 


0.5727(23) 


0.8(1.1)% 




0.02 


0.598(5; 


0.5935(20) 


0.8(0.9)% 


0.03 


0.01 


0.600(6; 


0.5957(22) 


0.7(1.1)% 




0.02 


0.617(4; 


0.6126(18) 


0.7(0.7)% 




0.03 


0.633(4; 


0.6287(16) 


0.7(0.7)% 


0.04 


0.01 


0.620(5; 


0.6160(24) 


0.6(0.9)% 




0.02 


0.633(4; 


0.6299(18) 


0.5(0.7)% 




0.03 


0.647(4; 


0.6435(15) 


0.5(0.7)% 




0.04 


0.659(3; 


0.6565(14) 


0.4(0.5)% 
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TABLE XXV: Results for the measured vector meson masses (DEG and UNI data sets). 



DEG 



UNI 



0.005 0.005 0.005 0.5053(58) 0.507(19) 
0.005 0.04 0.5591(75) 

0.04 0.04 0.6227(19) 0.6183(48) 

0.01 0.01 0.01 0.5288(45) 0.527(17) 
0.01 0.04 0.5887(89) 

0.04 0.04 0.6295(18) 0.6319(48) 

0.02 0.02 0.02 0.5789(55) 0.584(16) 
0.02 0.04 0.612(13) 

0.04 0.04 0.6453(33) 0.6453(84) 

0.03 0.03 0.03 0.6317(33) 0.6239(76) 
0.03 0.04 0.6447(83) 

0.04 0.04 0.6622(27) 0.6609(69) 





TABLE XXVI: Results for the measured of couplings fy/fv 




aniq 


0.03 0.02 0.01 0.005 


physical 


zVd.o.f 


fp/fp 


0.6806(55) 0.6646(73) 0.645(13) 0.624(21) 


0.619(15) 


0.17 


fl/h- 


0.6893(55) 0.6781(56) 0.6667(43) 0.6570(68) 


0.6498(62) 


0.11 


/0 //<!> 


0.6947(36) 0.6933(38) 0.6881(22) 0.6866(32) 


0.6838(33) 


0.10 



TABLE XXVII: Previous quenched lattice results for the ratio of tensor and vector decay constants. 



Reference 


/J(2GeV) 
fp 


/J.(2GeV) 

fK' 


/;'(2GeV) 


Becirevic et al. [86] 
Braun et al. [100] 


0.720(24) (+16) 
0.742(14) 


0.739(17)(+3) 


0.759(9)(0) 
0.780(8) 
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(a) {b) (c) 



FIG. 1: (a) One loop correction to fK', (b) one-loop wavefunction-like contribution proportional to g^; (c) 
diagram contributing to the one-loop wavefunction (and mass) renormalization. The grey square denotes 
a flavor-changing axial- vector current operator, the black bullets in (b) represent the PP*7i vertex and the 
grey circle in (c) denotes a KKkk vertex from L^^ in Eq. (22) . 



^0 ^ ^u 

FIG. 2: One-loop contribution to Bk ■ The grey box denotes the insertion of the KKnn operator ^22 as 
defined in the text. 
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Molecular dynamics time Q histogram (normalized) 



FIG. 3: The left graphs show topological charge, as measured on each ensemble every 5th unit of molecular 
dynamics time. The right panels show normalised histograms of toplogical charge. Notice the width of 
the histograms decreases as the light quark mass is reduced. For both the left and right graphs, the light 
dynamical quark mass is m/ = 0.005, 0.01, 0.02, 0.03, from top to bottom. 
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FIG. 4: Integrated autocorrelation time for 12th timeslice of tiie unitary degenerate pseudoscalar correlator 
with local sources and sinks, {P^^,P^^), measured on the m/ = 0.005 ensemble. Measurements were done 
every fifth trajectory in the range 900-1530 and a bin factor of two was used in the analysis. 
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FIG. 5: Integrated autocorrelation time for 12th timeslice of tiie unitary degenerate pseudoscalar correlator 
with a local sink and a Gaussian source, {P^^,P^^), measured on the m/ = 0.005 ensemble. Measurements 
were done every twentieth trajectory in the range 500-4500 and a bin factor of four was used in the analysis. 
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FIG. 7: The residual mass as a function of the valence quark mass. The solid line is a linear fit to the four 
unitary points with mi = 0.005,0.01,0.02 and 0.03. The star is the extrapolation to nii = 0, which we quote 
as the residual mass in our simulations. 
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FIG. 9: Sea quark mass dependence of the decuplet baryon masses for degenerate valence quarks rux = 0.03 
(black circles) and 0.04 (red squares). Linear fits are shown by sohd Unes. Blue asterisks denote masses 
extrapolated to the physical value of the bare, light sea quark mass ruud- 
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FIG. 10: Combined SU(2) x SU(2) fits for tlie meson decay constants (left panels) and masses (right panels) 
at two different values for the ligiit sea quark mass, and with a cut on the valence masses m^vg < 0.01. Points 
marked hy filled symbols were included in the fit, while those with open symbols were excluded. 
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FIG. 11: SU(2) X SU(2) fits for the kaon sector. Left panel: Kaon decay constant, right panel: Kaon mass 
squared for nig = 0.04. Points with filled symbols were included in the fit, while those with open symbols 
were excluded. 
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FIG. 12: Interpolation to the physical strange quark mass {red diamond) for mj^ and Jk- Open symbols are 
not included in the interpolations. Left panel: Kaon decay constant, right panel: Squared kaon mass. 
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FIG. 13: Finite volume correction factor obtained from the SU(2) fits for the decay constants (left panels) 
and squared meson masses (right panels) at sea quark masses mi = 0.005 {(top panels) and 0.01 (bottom 
panels). The vertical dashed line denotes the lightest (valence) mass used in this analysis {my = 0.001). 
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FIG. 14: Comparison of the size of the LO, NLO, and NNLO contributions in the SU(2) x SU(2) fits for 
unitary degenerate meson masses squared mfi (upper left panel) and decay constants /// (upper right panel) 
as a function of the the light sea quark mass parameter Xt (for details see Sec. VI C 1). The LO-contribution 
is normalized to one. Vertical dashed lines indicate quark masses of m/ = m„^/, 0.005, 0.01, and 0.02 (from 
left to right). The lower panels give a comparison of the size of the LO and NLO contributions in kaon 
ChPT for thekmnmas3sqUfaedm^f^(lower left panel) and kaon decay constant fth (lower right panel) as a 



function of the the light sea quark mass. 
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FIG. 15: Combined SU(3) x SU(3) fits for tlie meson decay constants (left panels) and masses (right panels) 
at two different values for the light sea quark mass, with valence mass cut m^vg < 0.01. Points marked by 
filled symbols were included in the fit, while those with open symbols were excluded. 
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FIG. 16: Combined SU(3) x SU(3) fits for tlie meson decay constants (left panels) and masses (right panels) 
at two different values for the light sea quark mass, with valence mass cut mavg < 0.03. Points marked by 
filled symbols were included in the fit, while those with open symbols were excluded. 
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FIG. 17: Leading order (LO) and Next-to-Leading order (NLO) contributions in the unitary meson decay 
constant fit in SU(3) PQCiiPT as function of the light quark mass parameter (left panels) and the heavy 
quark mass parameter (right panels). In the plots in the top panels the other quark mass parameter was 
set to our highest value used in the fits (m/ = 0.01 or m/, = 0.04), while in those in the bottom panels the 
remaining quark mass parameter was set to zero (m; or m/, = — mres)- NLO/ denotes the NLO contribution 
proportional to logXi s^d NLO//, the NLO contribution proportional to log(;^/ + Xh)/'^- The LO contribution 
is always normalized to one. Vertical dashed lines indicate quark masses of m/ = niud, 0.005, 0.01, 0.02, 
0.03, and 0.04 (from left to right, left panels only show quark masses up to 0.02). 
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FIG. 18: A plot of our measured values for the decay constant, converted to physical units, versus the 
degenerate valence pseudoscalar mass squared. The data for degenerate quarks is denoted by filled symbols 
and for non-degenerate quarks, open symbols are used. The graph shows that we see small effects for non- 
degenerate quarks. The results of SU(2) and SU(3) partially quenched ChPT fits to our data are shown, 
and both fits agree well with the data. The unitary SU(2) chiral extrapolation is also given, along with our 
value of /, and two of our data points lie on this curve, as expected. The value of / differs by ss 30% from 
the decay constant measured at mji = 420 MeV. We also plot the SU(3) chiral limit curve, for which the 
horizontal axis is the unitary meson mass for three degenerate mass quarks. None of our measured values 
must lie on this line. The large difference between /o and our measurements is apparent, showing the poor 
convergence of NLO SU(3) ChPT, with LECs as determined from our data. 
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FIG. 19: Extrapolation of the pseudoscalar bag parameter B^v in the light quark mass rux to the physical 
point niud for niy = 0.03 (left panel) and my = 0.04 (right panel). Points with filled symbols are included in 
the fit while those marked by open symbols are excluded.For fit parameters see Table XXII. 
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FIG. 20: Linear interpolation in niy of the pseudoscalar bag parameter Bu^y to the physical strange quark 
mass m.v (red diamond). The open symbol at my = 0.02 was excluded from the interpolation. 
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FIG. 21: Fitting the pseudoscalar bag parameter B^^, to the SU(3) x SU(3) x?T formulae with a cut in the 
average valence mass value of 0.01 (top panels) and 0.02 (bottom panels). Points marked by open symbols 
are excluded from the fit. The left panels are for light sea quark masses m/ = 0.005, the right panels for 
m/ = 0.01. 
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FIG. 22: Dependence of the fit parameter Bp^ (red bars) and the extrapolated physical value of Bk (blue 
crosses) on the average valence mass cut, when fitting to the SU(3) x SU(3) fit formula. The blue lines 
indicate the result (with statistical error) for Bk obtained from the SU(2) x SU(2) fit. 
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FIG. 23: Comparison of the pseudoscalar bag parameter Bx) for valence quark masses nix, rriy measured 
on 16"^ X 32 (red symbols) and 24^ x 64 (blue symbols) lattices at m/ = 0.01, m/j = 0.04. (Red symbols are 
slightly shifted to the right on this plot.) 
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FIG. 24: The diagram illustrating how T decays can be used to deduce fp and /k*. 
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a linear fit to the mass behaviour and the solid grey lines a quadratic fit. 



